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PREFACE 


This  document  comprises  the  final  report  for  the  program 
’'Automated  Estimation  of  Gas  Phase  Thermokinetic  Parameters 
for  High  Ignition  Temperature  Solid  Gun  Propellants."  The  work 
was  performed  for  the  Ballistic  Research  Laboratories,  Aberdeen 
Proving  Grounds,  Maryland,  by  McDonnell  Douglas  Astronautics 
Company,  under  Contract  DAAD05-  72-C-  0098  .  This  final  report: 
covers  the  period  from  15  November  1971  to  15  January  1973. 

The  report  is  submitted  in  two  volumes: 

Volume  I  —  A  FORTRAN  Program  for  Computing  the 
Thermochemical  Properties  of  Complex 
Gas  Molecules  by  the  Method  of  Group 
Additivity 
(Books  1  and  2) 

Volume  II  —  Definition  of  Requirements  for  Estimation 
of  Kinetic  Parameters 


ABSTRACT 


A  FORTRAN  program  for  calculating  the  thermodynamic 
properties  of  complex  gas  molecules  by  the  method  of  group 
additivity  is  documented  herein.  Included  are  descriptions  of 
the  computational  method,  char acteristic s,  functions,  input/ 
output  formats,  and  logic  structure  of  the  program.  Seventeen 
case  calculations  for  diversely  structured  molecules  are 
presented.  Also  depicted  are  a  listing  of  the  program,  a 
complete  glossary  of  the  program  variables  and  the  tables  of 
data  comprising  the  data  library  input. 
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Section  1 
INTRODUCTION 

The  increased  emphasis  in  recent  years  on  high  level  material  performance 
and  controlled  chemical  reactions  has  intensified  the  need  for  versatile  ana¬ 
lytical  programs  capable  of  estimating  the  thermodynamic  and  kinetic  prop¬ 
erties  of  compounds.  Moreover,  there  are  more  specific  reasons  why  these 
capabilities  are  important.  First,  the  size,  cost,  and  purity  of  the  sample 
and  the  expense  and  time  required  for  measurement  of  the  data  are  frequently 
excessive.  Secondly,  an  accurate  thermokinetic  computational  method  can 
serve  as  a  material  screening  device  and  predictive  tool  for  estimating  the 
properties  of  promising  new  compounds  and  postulated  chemical  structures. 

The  objectives  of  this  investigation  are  related  to  the  above  requirements. 

They  comprised  the  (1)  development  of  an  automated  program  and  (2)  prelim¬ 
inary  tasks  prerequisite  to  future  development  of  an  automated  program  for 
estimating  the  respective  thermodynamic  and  kinetic  properties  of  gas 
molecules.  The  results  achieved  in  the  present  study  are  twofold:  (1)  develop¬ 
ment  of  a  large  computer  program  for  calculating  the  thermodynamic  prop¬ 
erties  of  complex  gas  molecules  and  (2)  definition  and  extension  of  methods 
and  requirements  for  estimating  the  kinetic  parameters  of  several  classes  of 
gas-phase  chemical  reactions.  The  latter  work  is  described  in  Volume  II  of 
this  report.  The  future  goals  of  the  program,  which  were  instrumental  to 
the  initiation  of  the  present  study,  can  be  summarized  as  follows:  (1)  char¬ 
acterization  of  the  thermodynamic  profiles  of  high- temperature  solid  gun 
propellants  and  (2)  estimation  of  the  burning  rates  of  gaseous  -  solid  propel¬ 
lant  mixtures. 

The  computation  of  the  thermodynamic  properties  in  our  computer  program 
is  based  on  the  principle  of  additivity  of  group  properties  developed  by  Benson 
and  coworkers  (1  through  3)o  This  method  was  selected  for  its  applicability  to 
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both  chemical  thermodynamics  and  kinetics  and  its  accuracy  and  versatility 
in  handling  gas  molecules  with  complex  functional  groups  using  only  very 
simple  molecular  structural  data  as  input* 

A  complete  documentation  and  description  of  the  program  are  presented  herein. 
Sections  2  and  3  contain  descriptions  of  the  computational  method  and  other 
basic  properties  of  the  program  as  well  as  the  input/output  procedures. 

Section  4  contains  examples  of  case  calculations  representing  a  broad  spec¬ 
trum  of  organic  compounds.  The  logic  structure  and  computational  criteria 
of  the  program  are  discussed  in  Section  5  for  those  interested  in  probing 
further  the  logic  operations  of  the  program.  The  appendices  contain  a  com¬ 
plete  listing  of  the  program,  a  glossary  of  the  program  variables  and  a  listing 
of  the  program  data  library. 
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Section  2 

PRINCIPAL  PROPERTIES  AND  CHARACTERISTICS 


The  Thermochemical  Group  Additivity  Computer  Program  (TGAP)  is  a 
highly  automatic,  self-consistent  FORTRAN  IV  computer  program  for 
calculating  the  thermochemical  properties  of  complex  gas  molecules  as  a 
function  of  temperature  and  one  atmosphere  of  pressure.  The  program  is 
capable  of  handling  neutral,  non-f ragmented  molecules  composed  by 
hydrogen,  carbon,  oxygen  and/or  nitrogen  atoms. 

The  input  to  the  program  has  been  minimized  as  much  as  possible  for  facile 
and  rapid  utilization  not  only  by  physical  scientists  and  engineers  but  also 
by  persons  with  a  more  limited  knowledge  of  chemistry. 

2.1  COMPUTATIONAL  METHOD 

In  TGAP,  the  thermochemical  properties  of  the  molecule  are  calculated  by 
the  Group  Additivity  Method  of  Benson  (1  through  3)0  Unlike  the  more  simple 
atomic  additivity  techniques  for  the  evaluation  of  various  types  of  molecular 
properties,  the  group  additivity  method  is  primarily  an  empirical  second- 
order  approximation  technique  that  generally  yields  results  well  within  the 
limits  of  experimental  data. 

The  group  additivity  method  treats  a  molecular  thermodynamic  property  as 
a  composite  of  the  thermochemical  contributions  resulting  from  the  various 
atomic  groups  in  the  molecule  and  from  the  higher  order  structural  inter¬ 
actions  and  symmetry  of  the  molecule.  The  latter  include  those  contributions 
arising  from  nonbonded  next- near e s t- neighbor s  (cis,  ortho,  gauche) 
interactions,  ring  (strain)  corrections,  rotational  symmetry  and  optical  iso¬ 
merism.  In  the  assignment  of  the  group  property  values  themselves,  account 
is  also  taken  whenever  possible  of  differences  in  the  bond  environments  of 
next-near est-neighbor s  and  next-next-nearest-neighbor s.  For  example,  the 
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thermodynamic  properties  of  the  group  C- (H^CKC*),  are  different 
depending  whether  the  ligand  C*  is  bonded  to  another  carbon  atom  by  means 
of  a  single,  double  or  triple  bond. 

In  the  terminology  of  group  additivity,  a  group  is  defined  as  any  polyvalent 
atom  (ligancy>2)  in  a  molecule  together  with  all  its  ligands.  Any  molecule 
composed  of  two  or  more  groups,  such  as  CH^CH^  and  CH^CHOHCH^,  is 
amenable  to  treatment  by  the  group  additivity  method.  In  this  report,  the 
central  atom  in  the  group  is  referred  to  as  the  "core"'  atom  and  the  remain¬ 
ing  components  of  the  group  that  are  bonded  to  the  core  atom  are  frequently 
described  as  ’’ligands.  ” 


The  equation  for  the  calculation  of  the  thermochemical  property  F  by  the 
method  of  group  additivity  is 
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where  the  number  and  type  of  each  contribution  are  specified  explicitly. 

In  the  Thermochemical  Group  Additivity  Program,  the  data  for  calculating 
the  first  five  terms  of  the  above  equation  are  derived  from  the  tables  of 
group  properties,  ring  corrections,  ^nd  second-order  interactions  of 
Benson  (1  and  2)  for  hydrocarbons  and  oxygen  and  nit  rogen- containing 
compounds.  These  data  are  listed  in  Appendix  A. 
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2.  2  PROGRAM  FUNCTIONS 

TGAP  is  a  structurally  integrated  program  designed  for  minimum  data  input 
by  the  user.  An  important  feature  of  the  program  is  its  capability  of  per¬ 
forming  a  detailed  structural  delineation  of  a  molecule. 

For  the  purposes  of  the  calculation,  the  molecule  can  be  regarded  as  a 
topological  network  of  interconnected  atoms  which  can  form  linear  arrays, 
branch  junctions,  closed  loops  or  rings,  as  well  as  complex  combinations  of 
these.  The  computational  objectives  are  achieved  by  a  series  of  scanning, 
pattern  recognition,  comparison  tests,  and  other  subtle  logical,  topology- 
oriented  operations  which  enable  the  computer  to  identify  the  molecule  and 
recognize  such  structural  complexities  as  ring  closure,  branch  splitting, 
and  symmetry.  The  principal  logical  determinations  and  calculations 
performed  by  the  program  are  listed  below. 

A.  Identification  of  core  ligand  and  constituents. 

B.  Dete rmination  of  bond  order  and  direction. 

C.  Identification  of  single  and  fused  ring  structures. 

D.  Deter mination  of  longest  chain  in  the  molecule. 

E.  Detection  of  gauche  interactions  in  non-cyclic  groups  and  between 
cycLic  and  non-cyclic  structures. 

F.  Detection  of  cis  interactions  in  non-cyclic  groups  and 
non-aromatic  ring  structures. 

G.  Identification  of  ortho  interactions  in  aromatic  -  type  structures 

H.  Detection  of  dite rtiary- ethe r  structures, 

I.  Identif  ication  of  pyridine  -  type  structures  and  of  ortho  and  par  a 
substitutions  therein. 

J.  Determination  of  internal  rotational  symmetry. 

K.  Determination  of  external  rotational  symmetry  and  optical 
isomer  properties-  (i.  e.  ,  identification  of  enantiomers, 
meso  structures,  pseudoasymmetric  atom. 

L.  Calculation  of  the  individual  thermochemical  contributions  and 
the  gross  thermochemical  properties  of  the  molecule. 
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Operation  K  is  somewhat  limited  in  its  applicability.  At  present,  the 

external  symmetry  number  is  computable  for  molecules  whose  central 

atom  is  not  part  of  a  ring  structure  or  else  is  part  of  a  monocyclic 

benzene-type  ring.  If  the  molecule  contains  one  or  more  ring  structures 

but  the  central  atom  is  not  a  ring  atom,  the  symmetry  number  can  be 

calculated  also.  The  reason  other  centralized  ring  structures  cannot  be 

handled  is  that  virtually  every  ring  constitutes  a  unique  symmetrical 

situation  that  must  he  treated  independently  in  the  absence  of  molecular 

coordinate  data.  In  the  ca.se  of  optically  active  molecules,  the  logic  and 
•  .  *  ' 
mathematical  constraints  (1)  require,  the  presence  of  asymmetric  carbon 

atoms  and  (2)  take  into  account  only  those  asymmetric  atoms  which  form 

part  of  the  principal  chain.  Fortunately,  these  conditions  are  satisfied  by 

the  majority  of  such  molecules.  . 

The  program  executes  the  calculations  by  performing  a  (1)  topological 
structural  analysis  of  the  molecule  for  the  identification  of  structural 
groups,  second-order  interactions,  and  ring  structures;  (2)  a  symmetry 
anaLysis  of  each  atom  in  the  molecule  and  of  the  molecule  as  a  whole; 

(3)  a  specification  of  the  secondary,  ring,  and  symmetry  contributions  to 
the  thermodynamic  properties;  and  (4)  a  cross  match  of  each  group  with 
the  cor  re  sponding  group  values  contained  in  the  Data  Library  and  summa¬ 
tion  of  these  results. 

If  aLl  the  thermochemical  data  for  a  particular  group  value  are  missing  from 
the  Data  Library,  the  program  identifies  the  group  or  groups  involved,  and 
prints  a  warning  message. 

2.3  COMPUTATIONAL  ACCURACY  • 

The  thermochemical  properties  calculated  by  the  program  are  usually  in 
good  agreement  with  corre sponding  experimental,  values  up  to  the  upper 
temperature  limits  of  availability  of  the  empirical  data  (about.  1000°  to 
1500°K).  The  original  heat  capacity  data  of  Benson  are  recorded  at 
regular  temperature  intervals  up  to  1500°K  for  those  compounds  of 
interest  here.  These  data  have  been  curve  fitted  to  a  regular  third  order 
polynomial  for  utilization  by  the  program.  The  temperature  range  of  the 
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thermochemical  properties  ..calculated  by  the  program  extends  from  298°K 
to  5000° K.  Obviously,  the  validity  of  the  data  at  the  upper  temperature 
range  is  uncertain  particularly  since  complex  organic  gas  molecules  are 
unstable  at  these  elevated  temperatures.  However,  it  does  serve  to  provide 
extrapolated  data  limits  for  the  unstable  superheated  gas  state. 

A  number  of  calculations  executed  by  the  program  are  presented  in 
Section  4  that  illustrate  the  versatility  of  the  program  and  its  accuracy 
at  398°K  and  1000°K.  The  molecules  selected  for  this  purpose  represent 
a  broad  spectrum  of  structural  and  topological  properties. 

2.  4  DIMENSION  LIMITS  OF  PRINCIPAL  VARIABLES 

The  program  is  capable  of  handling  essentially  any  size  molecule  the  user 
is  likely  to  input.  The  following  is  a  list  of  the  prevailing  dimensional 
limits  on  the  size  of  the  molecule  and  its  structural  properties.  If  any  of 
these  dimensions  are  exceeded,  the  program  will  print  an  explanatory 


diagnostic  and  terminate  the  case  calculation. 

Core  atoms  in  molecule  100 

Number  of  unique  rings  40 

Core  atoms  in  ring  structure  30 

Number  of  chain  segments  60 

Core  atoms  per  chain  segment  50 


2.  5  COMPUTER  CORE  REQUIREMENTS 

The  program  was  developed  and  tested  on  a  Control  Data  6  500  computer 
with  the  RUN  compiler  and  CP  loader  under  the  scope  3.  3  operating 
system.  In  this  system,  the  program  requires  34595  or  103443g  core 
locations  exclusive  of  the  auxiliary  system  routines.  The  only  auxiliary 
mathematical  system  routine  required  is  the  standard  natural  logarithm 
subprogram  which  occupies  55  or  67g  core  locations. 

2.  6  RUN  TIME  REQUIREMENTS 

The  computation  time  requirements  are  a  function  of  the  size  and  complexity 
of  the  molecule.  On  the  CDC  6500,  the  thermochemical  properties  of  a 
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large  molecule  composed  of  27  carbon  atoms  anil  46  hydrogen  atoms  with 

i 

one  ring  structure  and  13  second-order  corrections  can  be  calculated  in 
approximately  3  seconds.  Molecules  of  more  conventional  size,  i.  e.  ,  15  or 

less  core  atoms  require  on  the  average  1.5  seconds  or  less.  Compilation 
of  the  program  takes  8  seconds  and  loading  3  seconds. 

Many  data  processing  systems  are  capable  of  accepting  a  compiled  deck 
or  a  compiled  program  file  directly.  This  type  of  input  requires  approxi¬ 
mately  3  seconds  for  loading. 
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Section  3 


INPUT/OUTPUT  PROCEDURES  .  ; 

This  section  describes  in  detail  the  input/output  data  formats  of  the 
program.  Basic  operational  instructions  are  also  presented. 

3.1  INPUT  FORMATS 

The  input  to  the  program  consists  of  (1)  punched  cards  comprising  data 
which  are  characteristic  of  each' molecule,  and  (Z)  the  TGAP  GROUP 
ADDITIVITY  DATA  LIBRARY  which,  is  a  fixed  data  package  permanently 
stored  in  the  program.  - 

3.  1.  1  Punched  Card  Input  .  ’ 

The  definitions  and  formats  for  the  punched  card  data  input  are  given  in 
Table  3-1.  As  shown  therein,  each  case  deck  is  followed  by  a  ***  card 
whereas  the  data  deck  itself  is  followed  by  an  END  card. 

The  Tables  of  .Thermochemical  Group  values  are  printed  at  the  very 
•beginning  of  the  first  case  calculation,  providing,  of  course,  the  output 
flag  is  appropriately  set.  Accordingly,  it  is  necessary  to  set  this  flag 
only  once  in  the  data  deck  of  case  one.  This  printout  operation  is  never 
repeated  when  subsequent  cases  of  a  run  are  executed  regardless:  of  the 
value  of  the  output  flag  in  these  cases.  The  user  also  has  the  option  to  print 
an  identification  label  and  to  eliminate  the  printout  of  the  program  title  and 
case  number.  The  program  title  printout  flag  can  be  individually  set  for 
each  case  calculation.  The  printout  of  the  case  number  in  all  case  calcula¬ 
tions  of  a  given  run  is  controlled  by  the  case  number  printout  flag  of  .case 
number  onee 

The'.program  always  computes  the  external  rotational  symmetry  number 
and  optical  isomer  contributions  unless  the  molecular  structure  is  of  the 

type  described  in  Section  2c  2,  in  .  which  case  a  warning  message  is 
printed.  In  any  event,  the  user  has  the  option  to  input  the  symmetry 
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PUNCHED  CARD  DA 

Column 

Number  Number  Type/Symbol. 

•  1  4-6  .  Numeric 


1 

7-9  . 

Numeric 

1 

10-12 

Numeric 

1 

13-15 

Nume  ric 

1 

31-80 

2  to  N 

1-80 

H,  C,  O,  N. 

(N  £  51) 

$  . 

N+l 

1-3 

Last  Card 

1-3 

END 

*13  data  are  right  adjusted  in  the  data  field 


Table  3-1 

t  FORMAT  FOR  PROGRAM  TGAP 


Format  Description 

13*  Printout  options  (OJK): 

K=0  no  printout  of  Data  Library 
K^O  data  are  printed. 

J=0  title  and  case  number  printed. 

1  title  and  case  number  not  printed. 
J=2  title  -  yes;  case  number  -  no. 

J=3  title  -  no;  case  number  -  yes. 

13  -External  rotation  symmetry  number 

(optional). 

13  Number  of  enantiomers  (optional). 

13  Number  of.  me  so  structures  (optional). 

50A1  Molecule  identification  label  (optional). 

80A  1  Graphic  formula  of  the  molecule. 

A3  Denotes  end  of  case. 


A3 


Denotes  end  of  data  deck. 


number  and/or  the  asymmetric  properties  (number  of  enantiomers  and 
meso  structures)  directly.  These  are  utilized  by  the  program  to  compute 
the  corresponding  contributions  to  the  entropy.  The  program  then  prints 
both  the  input  and  computed  preliminary  parameters  and  entropy  contributions. 
When  the  symmetry  number  and/or  optical  parameters  are  input,  the  entropy 
contributions  derived  from  these  data  are  used  to  compute  the  total  entropy 
and  free  energy  of  the  molecule.  In  the  more  general  case  when  these 
parameters  are  not  input,  the  parameters  calculated  by  the  program  are 
used  instead  of  the  input  parameters  to  calculated  the  corresponding  molecular 
properties. 

Table  3-2  depicts  the  FORTRAN  symbols  available  for  inputting  the  graphic 
formula  of  the  molecule.  The  codes  for  the  chemical  atoms  and  radicals  are 
identical  to  their  conventional  chemical  symbols.  The  graphic  formula  is 
constructed  essentially  the  same  as  the  handwritten  version.  The  following 
rules  describe  the  procedures  for  punching  the  graphic  formula  onto  the 
input  cards. 

A.  Always  use  the  designated  bond  symbols  to  represent  the  bonds 
between  atoms,  but  do  not  insert  bonds  between  the  atoms  of  a 
radical.  For  example,  input  -CN  but  not  -C-N  lor  the  radical  CN. 

B.  Chemical  symbols  must  be  placed  at  the  end  of  the  bond  symbol, 

but  not  on  the  side.  For  example  C - H  not  G - 

C.  The  number  and  types  of  bonds  assigned  to  ah  atom  or  radical 
must  be  sufficient  to  satisfy  it’s  valence. 

D.  There  are  a  maximum 'of  eight  bond  positions  available  about  a 
particular  atom,  as  shown  in  Figure  3-1.  If  the  atom  is  located 
near  a  card  boundary,  the  number  of  locations  that  are  available 
for  bonding,  is,  of  course,  reduced. 

E.  The  bond  between  two  chemical  entities  can  be  any  length  but  its 
direction  must  be  linear.  Absolutely  no  deviation  from  linearity  is 
allowed  when  constructing  a  bond. 

F.  There  is  considerable  freedorfi  in  the  way  the  bonds  connecting  a 
radical  can  be  positioned.  Table  3-3  shows  various  permissible 
positions  for  the  different  radical  types  listed  in  Table  3-2.  The 
actual  order  of  the  atoms  in  the  radical  does  not  matter,  for 
instance  -N02,  -N20,  .-2NO,  etc.,,  are  all  acceptable. 
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Table  3-2 

SYMBOLS  FOR  SPECIFICATION  OF  GRAPHIC  FORMULA 


Atoms 

H 

C 

O 

N' 

Radicals 
.  CO 
CN 
NO 
NO  2 

Bonds 
Single  - 
Double  = 
Triple  $ 


column  - ^ 

CARD 

1 


7  ■■  - .  .  M  - ^  3 


CR44 


Figure  3-1.  Maximum  Number  of  Directions  Available  for  Chemical  Bonding 
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Table  3-3 

PERMISSIBLE  ATOMIC  POSITIONS  FOR  RADICALS 
Type  1.  Components  =  2,  Connectivity  =  1 


N  a 

C- 

N 

C- 

NC- 

C- 

N 

C- 

N 

45°  <  a  <  315° 

(a) 

(b) 

(c) 

(d) 

Type  2.  Components  =  2,  Connectivity  =  2 


O 

o  . 

V* 

-CO 

-C- 

-C- 

-C- 

c- 

o 

o 

(a) 

(b) 

(c) 

(d) 

(eH 

Type  3. 

Components 

=  3,  Connectivity  =  1 

2 

,  2 

.  -  N  . 

O 

O 

O 

-N 

-N 

-N02 

2 

(a) 

(b) 

(c) 

(d)  . 

Also 


N 

-O 

2 


45°  <  a  <  3  1  5C 


*Do  not  use  if  radical  is  in  the  top  leftmost  location  of  the  molecule. 


Atoms  and  radicals  cannot  be  adjacent  to  bonds  to  which  they  are 
not  connected.  Thus  the  linkage 
=C-H1 

"H2 

is  not  allowed  because  is  also  adjacent  to  the  bond  between 
C  and  H^.  r  This  is  easily  corrected  by  adding  an  additional  bond 
symbol  between  C  and  H^,  namely,  • 

=G- -H. 


1 


H- 
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This  rule  is  perhaps  the  easiest  to  overlook  and,  consequently, 
should  he  given  special  attention  when  formulating  the  structural 
input.  In  the  event  any  of  these  rules  are  violated,  a  diagnostic 
message  is  printed,  the  current  case  calculation  is  aborted,  and 
the  next  case  calculation  (if  any)  is  begun. 

3.  1.  2  Data  Library  Input 

The  TGAP  GROUP  ADDITIVITY  DATA  LIBRARY  is  composed  of  Benson’s 
thermochemical  group  additivity  values  and  interaction  corrections  for 
hydrocarbons  and  compounds  containing  oxygen  and  nitrogen.  These  data 
are  listed  in  Appendix  A.  Except  for  the  contributions  of  rotational  sym¬ 
metry  and  entropy  of  mixing  due  to  optical  isomerism,  which  are  handled 
separately,  the  DATA  LIBRARY  constitutes  the  basic  data  bank  for  the 
calculation  of  the  thermochemical  properties  of  the  molecule.  Once  the 
specific  molecular  groups  and  interactions  are  identified  by  the  various 
logic  routines  of  the  program,  the  group  contributions  and  corrections 
specified  by  the  program  are  abstracted  from  the  DATA  LIBRARY  and 
summed  to  determine  the  gross  thermochemical  properties  of  the  molecule. 

The  principal  components  of  the  DATA  LIBRARY  are  the  Tables  of  Thermo- 
chemical.  Group  Values  which  the  program  will  output  when  so  directed  by 
the  user.  A  listing  of  this  output  is  given  in  Table  A-l.  The  remaining  data 
consist  principally  of  the  thermochemical  corrections  for  ring  structures 
(Table  A-2)  and  the  thermochemical  contributions  for  various  types  of 
gauche,  cis,  and  ortho  interactions  (Table  A -3). 

The  data  in  the  LIBRARY  are  permanently  stored  in  the  program  by  means 
of  data  statements.  These  are  positioned  in  the  block  data  subprogram 
DATA1  and  Subroutine  CORCIG. 

3.2  OUTPUT  FORMAT  . 

All  output  occurs  via  the  printed  mode  and  comprises  both  input  and  calcu¬ 
lated  data.  In  addition  to  the  thermodynamic  properties,  the  output  of  the 
program  provides  information  on  the  structure  and  symmetry  properties  of 
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the  molecule.  These  data  are,  in  essence,  by-products  of  the  - 
thermochemical  computations.  The  individual  output  arrays  are  described 
below  in  the  order  in  which  they  are  printed  by  the  program. 

3.  2.  1  Molecular  Structure 

This  is  a  printout  of  the  graphic  formula  of  the  molecule  in  exactly  the 
same  form  it  is  input  into  the  program.  The  structure  is  automatically 
centered  on  the  printout  sheet  regardless  of  where  it  appears  on  the  input 
card. 

3.  2.  2  Atomic  Composition 

The  types  of  atoms  in  the  molecule  and  the  number  of  each  are  printed  next. 
This  includes,  of  course,  the  atomic  constituents  of  each  radical. 

3.  2.  3  Group  Structure  and  Composition 

This  is  a  printout  of  a  multi-group  property  array;  it  comprises  the 
following  items. 

A.  Group  number  —  the  ordinal  number  assigned  to  each  group. 

This  number  identifies  the  group  as  well  as  the  core  atom  thereof. 

B.  Group  weight  —  an  integer  which  is  unique  to  each  chemical  group 

and  which  correlates  the  group  with  its  thermochemical  group 
values.  ' 

C.  Chemical  symbol  of  the  individual  components  of  the  group. 

D.  Row  and  column  coordinates  of  group  components.  For  single 
atoms,  the  row  and  column  coordinates  represent  the  particular 
card  and  column  in  the  structural  input  deck  containing  the  species. 
For  radicals,  the  coordinates  may  be  shifted  one  or  two  locations 
away  from  the  input  coordinates  since  the  program  ultimately  packs 
all  the  symbols  of  a  radical  into  one  word. 

E.  Bond  vector  —  the  bond  direction  between  a  ligand  and  the  core  atom 
in  a  group.  The  numbers  represent  the  bond  directions  depicted 

in  Figure  3-  1 . 

F.  Bond  type  —  the  multiplicity  of  each  chemical  bond  in  the  group. 

G.  Group  number  of  core  atom  —  this  number  identifies  those  ligands 
in  a  group  which  are  also  core  atoms. 
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3.2.4  Ring  Specifications  ,y 

If  ring  structures  are  present,  the  program  prints  the'  total  number  of 
members  in  each  unique  ring  in  the  molecule  together  with  the  group 
number  of  each  member. 

3.  2.  5  Ligand  Identification 

This  printout  identifies  each  unique  pair  of  complex  ligands  in  the  molecule 
which  have  identical  structures.  The  corresponding  core  atom  is  also  speci¬ 
fied,  Here  the  term  ligand  refers  toall  the  atoms  in  a  particular  chain  or 
branch  which  is  bonded  to  a  specific  core  atom.  The  ligand  is  considered 
complex  if  it  has  one  or  more  core  atoms,  and  not  merely  atoms  with 
connectivities  of  one  like  hydrogen,  CN,  etc. 

3.  2.  6  Atomic  Composition  of  Group  Components  . 

In  this  printout,  the  atomic  composition  of  each  ligand  and  core  atom  in  a 
group  is  specified.  Here  the  term  ligand  has  the  same  significance  as  in 
3.2.5.  If  the  core  atom  happens  to  be  part  of  a  ring,  the  atomic  dis tribution  of 
the  ligands  which  form  part  of  the  ring  structure  is  quite  arbitrary.  In  either 
case,  the  total  number  of  atoms  of  a  particular  type  which  are  contained  in  a 
group  should  correspond  to  the  total  number  of  such  atoms  in  the  molecule. 

3.  2.  7  Properties  of  the  Longest  Chain 

This  printout  gives  the  number  of  core  atoms  present  in  the  longest  chain 
in  the  molecule  and  the  group  number  of  each  atom. 

3,  2.  8  Nongroup  Interaction  and  Symmetry  Contributions 

The  nongroup  interaction  correction  for  the  heat  of  formation  of  ring  struc¬ 
tures  includes  a  correction  for  the  strain  energy.  The  majority  of  the  ring 
corrections  are  for  rrionocyclic  rings  whose  ring  number  is  printed  in  column 
two.  A  few  are  for  dicyclic  and  tricyclic  fused  ring  structures. 

Other  nongroup  interaction  corrections  include  ortho,  para,  cis,  gauche 
alkene,  and  ditertiary  ethe  r  corrections.  In  this  case,  the  group  numbers 
of  the  core  atoms  involved  in  the  interaction  are  printed  in  column  two  and 
the  number  of  interactions  of  each  type  is  displayed  in  column  three. 
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The  symmetry  corrections  refer  to  internal  rotations  of  twofold  and 
threefold  symmetry. 

3.  2.  9  External  Rotational  Symmetry  Contribution 

The  calculated  symmetry  number  and  the  corre sponding  entropy  contribu¬ 
tions  are  printed  next.  If  the  symmetry  number  is  not  calculable,  it  is 
indicated  in  the  printout.  In  this  case  the  calculated  symmetry  number  is 
set  to  one  and  the  entropy  to  zero.  If  the  symmetry  number  is  input,  it  is 
also  output  together  with  the  entropy  contribution  computed  therefrom. 

3.  2.  10  Optical  Isomer  Contribution 

This  contribution  accounts  for  the  entropy  of  mixing  associated  with  the 
presence  of  optical  isomers  and  asymmetric  carbon  atoms.  The  program 
computes  and  prints  the  number  of  asymmetric  and  pseudoasymmetric 
atoms,  enantiomers,  meso  structures  and  the  corresponding  entropy 
(Section  5.  3.  7.  5).  If  these  properties  are  not  calculable,  it  is  indicated 
in  the  printout  and  all  pertinent  values  are  set  to  zero.  If  the  number  of 
enantiomers  and  meso  compounds  are  input,  they  are  also  output  together 
with  their  entropy  value. 

The  program  also  outputs  the  group  numbers  of  the  asymmetric  atoms  and 
of  the  pseudoasymmetric  atom  (if  any). 

3.2.  11  Thermochemical  Properties  of  the  Molecule 

The  coefficients  for  the  heat  capacity  of  the  molecule  are  those  for  the 
regular  third-order  polynomial 

Cp  =  cQ  +  CjT  +  c2T2  +  c^T3 

°  0 

The  data  in  the  DATA  LIBRARY  are  expressed  in  terms  of  AHf?  298>  S298’ 
and  the  four  heat  capacity  coefficients.  The  remaining  thermochemical 
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properties  of  the  molecule  are  computed  by  melts  of  the  standard 
relationships  below: 

HT  -  H0  =  mo  (coT  +  1/2  =1t2  +  1/3  c2T3  +  1/4  c3T4) 

scons  =  S298  '[coT  +  1/2  C1  (298)  +  1/2  c2  (298>2  +  1/3  c3  <298> 

S°  =  c  InT  +  1/2  C?T2  +  1/3  C,T3  +  s 
1  °  ^  3  cons 

-  (G^,  -  H-98)/T)  =  -  [lOOO  (H^  -  H"98)  -  TS^,Jy/T 
Hi  (To, al)  =  4Hfi  2,8  +  H^.  -  H°98 
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Section  4 

EXAMPLES  OF  COMPUTATIONAL  OUTPUT 


This  section  contains  complete  listings  of  17  thermochemical  calculations 
performed  by  the  program  for  a  variety  of  different  molecular  structures. 

As  shown  in  Table  4-1,  the  examples  include  hydrocarbons,  saturated  and 
unsaturated,  symmetric  and  asymmetric  molecules,  as  well  as  fused  and 
monocyclic,  aromatic  and  nonaromatic  ring  structures.  Nearly  all  types 
of  nongroup  interactions  and  symmetry  contributions  are  represented. 

Table  4-2  presents  a  comparison  of  the  calculated  and  published  thermo¬ 
dynamic  properties  of  12  of  these  molecules.  The  agreement  in  the  majority 
of  cases  is  very  good. 

Case  17  requires  some  special  comments.  The  molecule  in  question  is  a 
highly  complex  hypothetical  structure  formulated  with  the  objective  of  demon¬ 
strating  the  program’s  symmetry  detection  capabilities.  It  has  the  basic 
structure  CX^YZ,  where  carbon  atom  C  is  defined  by  group  number  12.  An 
analyses  of  the  printout  reveals,  that  the  two  X  ligands  (group  numbers  9  and 
16)  have  been  identified  correctly  as  well  as  the  six  subordinate  pairs  of 
complex  ligands.  These;data  were  used  to  compute  the  listed,  symmetry  and 
optical  isomer  properties  of  the  molecule.  The  computation  time  for  case  17 
was  approximately  3  seconds  on  the  CDC-6500  computer. 

Table  4-3  lists  the  input  data  for  all  17  cases  in  the  same  order  as  the 
computational  output. 
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LIST 

Table  4-  1 

OF  SAMPLE  MOLECULES 

Case  Number 

Name 

i 

Ethane 

2 

Octyl  alcohol 

3 

Valeraldehyde 

4 

tert-Butyl  ether 

5 

Butyronitr  ile 

6 

1  -  Cyclohexylheptane 

7 

Phenol 

8 

2,  2,  3  -  Tr imethylhexane 

9  ■ 

Styrene 

10 

Aniline 

11 

Naphthalene 

12 

cis  -1  -Methyl  —  2  -  ethyl  cyclopentane 

13 

T  rimethylmethane 

14 

Benzenehexacarbonit  rile 

15 

bis  (2-nitro  -5  —  hyd  roxy  cyclohexyl)  - 

malononitrile 

16 

F used  nitrogen  tricyclic  ring 

17 

Hypothetical  molecule 

I 
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Table  4-2 


COMPARISON  OF  CALCULATED  AND  PUBLISHED  THERMOCHEMICAL  PROPERTIES 


Case 

Number 

- - - - 1 

1 

Temperature 

°K 

AH°  (calc) 

AH" 

S° (calc) 

S° 

C  ° (calc) 

P 

C" 

P 

Ref 

Kcal/  mole 

cal /  0K-mole 

i 

298 

-20.  16 

-20.  24 

55.  08 

54.  85 

12.  27 

12.  58 

4 

i 

1000 

79.  75 

79.  39 

29.  55 

29.  33 

4 

2 

298 

-86.  16 

-85.  34 

123. 85 

124. 14 

48.  18 

48.  17 

4 

2 

1000 

217. 11 

217. 08 

107. 72 

107. 58 

4 

3 

298 

-54.  58 

-54.  45 

91. 60 

91.  53 

30.  28 

29.  96 

4 

3 

1000 

147. 70 

147. 69 

63.  86 

64.  00 

4 

4 

298 

-87.  78 

-87. 20 

105.  18 

102.  12 

48.  92 

48.  76 

4 

4 

1000 

201.  02 

195.  26 

109.  77 

107. 86 

4 

5 

298 

7.  47 

8.  14 

77.  85 

77.  78 

22.  63 

23.  19 

4 

5 

1000 

120. 56 

120.  40 

48.  40 

48.  22 

4 

6 

298 

-65.  63 

-65.  73 

138. 58 

137.  51 

64.  27 

65.  89 

4 

6 

1000 

274.  63 

274. 79 

162.  61 

162.  70 

4 

7 

298 

-23.  20 

-23.  03 

74.  17 

75.  43 

24.  26 

24.  75 

4 

7 

1000 

124. 1 1 

124.  65 

55.  50 

55.  49 

4 

8 

298 

-57.  70 

-57. 65 

112.  00 

111.  34 

50.  49 

50.  43 

4 

8 

1000 

213.  00 

212.  42 

117.  29 

117.  38 

4 

9. 

298 

35.  22 

35.22 

82.  46 

82.  48 

29.  20 

29.  18 

4 

9 

1000 

141.  52 

141.  50 

67.  83 

67.  92 

4 

10 

298 

20.  80 

20.  76 

76.  29 

76.  28 

25.  70 

25.  91 

4 

10 

1000 

128. 17 

128.  38 

59.  05 

59.  18 

4 

11 

298 

36.08 

36.08 

80.  22 

80.  22 

31.  67 

31.  68 

4 

1 1 

1000 

147. 41 

147.  36 

78.  30 

78.  38 

4 

12  1 

298  j 

-35.  66 

-35.  9 

_ 1 

1 

GROUP  ADDITIVITY  THtRHOCWEHU^V,  FRORtRfjjfe'  RPt>D*AM 
- CASE — NUMBER _ 1 - - 

. .  . MOLECULAR  .STRUCTURAL  --■  — . — 

- - - - - H - ~ - - - 

m  *  y 

H  H,  ! 


atomic  composition  op  molecule 

. ATOM . . NUMBER  _ 

H  6 

- c - 3 - — 


GROUP’  STRUCTUHE  AND  COMPOSITION 


-ORQU?. _ QKQLP CHEKICAUS.YM3DL ^OR  til  RQH  .^fijtLD-COLLMN - BCNIL BONO — OR  CUE -NUMBER  — 

Number  weight  of  component  coordinate  coordinate  vector  type  of  core  atom 


c 

3 

’  3 

1 

. C  -  ----- 

-  3  ,  .  ... 

6 

.  3  .. 

1 

-  -  -  -2  - 

y 

3 

i 

7 

1 

0 

y 

1 

3 

— . -  3- 

1 

1 

- 1 - 

t 

1 

! 

T*  ° 

i 

i 

; 

2  296 


ATOMIC  COMPOSITION  OF  CORE  ATOM  AND  LlGANCS  IN  EACH'  OROUP 


GROUP  NUMBER  SUBGROUP 

1  CQHt 

. -  •  LI  GANG 

LIGAND 

LIGAND 

- - - ’- - :L1GAND— 

2  CCHt 
LIGAND 
LIGAND 
LIGAND 
LIGAND 


H 

0 

J 

1 

1 

1  - 
0 
i 

1  . 
1 
1 


C  C  M 

1  0  0 

-  1  -  ...  0  0  • 

C  0  c 

C  0  0 

- o - 0 - 0-  ’ - 

1  0  0 

1  0  o 

-  0  0  0 

0  0  0 

C  0  0 


. NUMBER  OF  CQR£  ATOMS  IN  LONGEST  CHAIN-  --  —  2  --  -  •  - 

GROUP  MIJNBER  OF  CORE  CONSTITUENTS  CF  LONGEST  CHAIN  •  1«  2, 

.  -  -  _  NOUbKCLP  INTERACTION  AND  INTERNAL  ROTATIONAL  SYMMETRY  CQNTRIBUT IONS 

TYPE  RING  NUMBER  CR  AMOUNT  HEAT  np  FORMATION  bNT«QPY  HEAT  CAPACITY  COEFFICIENTS 

" . . . - — OR QUP--Num.fi . . .  -  -  KCAL  -  - - •  - -CAL/Dfcfi -W -  CAL/CEG  K  CaL'DEG  *  *  *2  CAL/CEG  *•  *3 -CAL/C6G--M** 

INTEHNAL  RUTATION  2  1  -2*183 


EXTERNAL  ROTATIONAL  SYMMETRY  CONTRIBUTION 

S OUNCE  SYMMETRY  NUMBER  ENTROPY 

CCMPL'TEC  6  -3,56069 


SUl-'PCE 

COMPUTED 


- - XPT1CAL-  ISOMER-  CCNtHlBUT  TON - -  -  -  - 

ARYMWTRIC  PSCUDCASYMM.ETHJC  ENANTIOMERS  HESC  TOTAL  OPTICAL  isomers 
0  0  0  0 


ENTROPY 
0  ,  00000 


TMEPMCCFibHlCAL  PROPERTIES  cr  MOLECULE 

heat  capacity  CCEfFlClfcUS 

CAL/DfOMOLE  CAL/CEG**2-MOLE  CAL/0fcfi**3.MCU  C  AL/OEn***-vrLE 

-T,3tJEA?6eeE-01  4,eM56852E-H?  -  2 , 334?6  7*  ?  b •  (,  5  « , 36p70?70b* 09 


T 

CtP) 

S 

H(T).HfG) 

• (Q( T )»H<29e ) )/! 

M  Y ) »w ( ?90  ) 

DELHF  f296  >*M  (  T  NM298  I 

GEU  K 

-  CAL/DEG-MOLE 

-  CAL/DEC-MOL6  • 

KCAL/HOLE 

-  CAL/DEC-MCLE 

kcal/mole 

KCAL/MCLF 

#i) 

I2,?.665b6 

55,076076 

1,672337 

55,076070 

0,000000 

*20,160000 

300,0 

'  12,340730. 

55.15P384 

1jB96V46 

95,0763?? 

,024609 

•20.135391 

_ A  OC-.O-- 

_ 15 ,73465a _ 

_ 5 9  rl E  Si 16 _ 

_ "2 ,304843. 

.  55r6ft3P5l _ 

_ 1-^4325  06 _ 

—  •$0,727494.  -  - 

500,0 

IE, 606365’ 

63,034981 

5*035803 

56,7Ceit4fl 

3,163466 

•16,996534 

600,0 

21.522042 

66,70  F.  734 

7,054915 

56,071104 

5,182578 

•14,977422 

.703,0 

23,927046  .  . 

7u,2Jnfi4G 

9,32930? 

59,557211 

7,457546 

•  12.702454  .  ... 

*00,0 

26,049943 

73.347506 

11,031027 

61,099145 

9, 95669 H 

•10.201310 

930,0 

27,914496 

76,725799 

14,531206 

62,6<5C300 

12,650949 

•7,501051 

on  ka-t/laq  . 

70  _ 

_ 17 ,406214- 

_ 44 ,£59490  - _ 

4-5^  533^7  7- 

_ _ •4-.626123— .  - 

i2:c'o 

— 9  }W*FOO  » — ■ — 

32,224534 

, — 1  - 1  r*  *  *  O  ^  t  ~ 

65,387997 

23,599371 

07,245469 

71.723334 

1,563034 

14T0  ,0 

34  ,2*9852 

9H, 517041 

30,255257 

70,243526 

28,352921 

6 , 22292 1 

160  u  r  U 

35,952937  - 

.  95,207694 

37,28449ft 

73,075094 

35,41219? 

15,252159 

1030, 0 

37,4^3101 

99 , 32T091 

44,623571 

75,777465 

42,751234 

22,591234 

2300,0 

30,909660 

103,546886 

52*254029 

78,3556*0 

50,382493 

30,222493 
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GROUP  ADC  IT  IV  ITT  THERMOCHEMICAL  PROPERTIES  PROGRAM 


-  _  .  .CASE  ■  NUMBER - »- 

— .  -  -  --  HOUUULAR— STRUCTURE- 

- ■ - M-H„N  M  -M-M--*-  H - 

*••«*••* 

H*c*t;-c*C-C»C*C"C-0*H 

HHHHHHHH 


atomic  composition  of  molecule 
atom  "number 

h  te 

- . — ~c - ' - o - 

u  i 


GROUP  STRUCTURE  and  COMPOSITION 


GROUP 

NUMBER 

—  1-  — 


GROUP 
HEIGHT 
-  -296  - 


.  359 


399 


399 


399 


chemical  symbol 
or  component 
-  c 
c 

h 


URID  ROM  GRID  CCLLKN  BON’D 
COORDINATE  COORt IN* TE  VECTOR 

3  "  13  "~3~ 

3  9  7 

- 1 - - 11 - 1 - 

5  11  5 


BONO 

TYPE 


GROUP  NUMBER 

or  core  atom 

- 1 - 

l  2 

1  0 

1 - 0 - 

1  0 


13 

11 

15 

13 

13 

.15 

13 

17 

-15- 

15 

17 

15 

19 

-17- 

17 

—19 

17 

21 

-49- 

19 

21 
19 
23 
21 
21  . 

23 
21 
25 
23  - 
23 

25 

23 

27 

-25- 

25 

27 

25 

29 


7 

3 

..-1  . 
5 


1 

1 

-4— 

1 
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ATOMIC  COKPCS I T I  ON  OF  CORE  ATOM  AND  LlfiANDS  '  IN  EACH  QRQUP 


CROLP-MJy.E£fv- 

1 


..5  _ 


..SUBGROUP.  . 
CCHE 
LIGAND 
LIGAND 
LIGAND 
ligand 

_ COKE™ 

LIGAND 
LIGAND 
LIGAND  — 
LIGAND 
CONE 

-  -  ligand — 
ligand 
ligand 

LIGAND 

COME 

LIGAND 

— ligand-- 

LIGAND 

LIGAND 

COKE 

LIGAND 

ligand 

--LIGAND-- 

LIGAND 

cone 
LIGAND 
UGAND 
LIGAND 
—ligand- 
CORE 
LIGAND 
LIGAND 
LIGAND 
LIGAND 
— -  GORE  - 


_  h _ 

o 

15 

1 

1 

1 

__0— 

13 

3 

-  1 
1 
0 

-11 

9 

1 

1 

0 

9 

1 

1 

0 

7 

V 


l 

0 

5 

11 

1 

. 1~ 

0 

3 

13 

1 

1 

—  0  - 


....  c 

1 

7 

0 

0 

0 

— 1 
6 
1 
c 
0 
1 

...s... 

2 

D 

0 

1 

4 

— 3- 
0 
Q 

1 

3 

4 

— D-- 
0 
1 

-  2 

5 
D 

1 

1 

6 

P 

0 

—  1- 


1 

0 

0 

0 

0 

-.1 — 

0 

0 

O' 

0 

1 

—ft — 


Q 

0 

-  0 
1 

•o 

- o  - 

0 

0 

■  1 

0 
0 

- 0— 

0 

1 

..  o  . 

0 
0 

- o .. 


_ N  . 

0 

O' 

c 

0 

0 

—  0.. 

o’ 

0 

.  0 
0 
0 

0, 

0 

0 

O' 

0 

0 

-0- 


LIGAND 

1 

0 

1 

0 

ligand 

is 

7 

0 

0 

ligand 

1 

0 

c 

0 

LIGAND 

1 

0 

0 

0 

COHfc 

0 

0 

1 

0 

ligand 

17 

6 

0 

0 

LIGAND 

1 

C 

0 

0 

NUMREI»  DF  core  atoms  in  longest  CHAIN  «  9 

GROUT  l.luBtfi  Of  CODE  CONSTITUENTS  CF  LONGEST  CHAIN  »  1, 


2. 


3, 


TYPE 

jntehnal  rotation 


NONGRCLP  I  f  TFRACT  ION  AND  INTERNAL  HCTaIIONAL  SYMMETRY  CCN  TK  It-UT  1 CNS 


Pjf»G  NI.MPfch  CR 
GROUP  NUMI-EA 

1 


AHCLNT 

1 


HLaT  OF  FORMATION  tsTROPY 
KCAL  CAL/UtG  K 

*2,1*3 


h£  *  T  CAPACHv  COEFFICIENTS 
CAL/DtC  K  CAL/DEG  K * • 2  CAL/CEG  K**3  CAL/CEG 


EXTERNAL  ROTATIONAL  SYMMETRY  CCMP  !  feL  T  I  CD 

SOURCE  SYMMETRY  NUMBER  ENTRCFY 

CCMPUTtr  1  o.ooocc 


OPTICAL  ISOMfcR  C0NTHJ8UT10N 

SOURCE 

C  0  P  L‘  T  E  C 

ASYMKE  TR  IC 

0 

PSEITCASYKMFTRIC  ENANTIOMERS  r  ESC 

C  0  0 

TOTAL  OPTICAL  ISOMERS 

0 

f  N  T  R  CP  T 
0  ,OQCOfl 

TrERSCCHE-N  I  CAL  PROPERTIES  CF  HClECULE 


Tic  at  capacity  coefficients 


CAL /DEG 

-^CLF  CAL/EEG**2-me 

CAL/DfcG**3*M0lE  - - 

•€AL/DtG**4-MCLE 

•  l ,  0  2  05  P  5  2  F.  •  u  1  1 .  P 19C7042E-- 1 

•1.07097057L-C4 

<  ,3011  9994E-O0 

T 

C{P) 

H  (  T  )  -  H I  O 

•(GIT)-hl29P))/T 

H  T  I  -  n  {  2  9  6  > 

DEIHF ( 298 ) *H  <  7 ) »h 1 

DEL  K 

CAL/Cfi-MOLE 

cal/cfc-hdle 

kcal/mplc 

C  A  L / r  F  r.MCLE 

XCAL/MCLE 

KCAl/PCLE 

,d 

•  «ft,ii:34?6  - 

■-  l?3.M*77'j  -- 

7,59.194ft 

- 123 ,04677  0  - 

—  OrOOOnon- . 

-  •06,140000 

3?C,0 

48 , 4  ’1496 

1?4 , J  6  9  <5  ft  4 

7 , 60  7  5E 1 

123  ,  B47B4B 

.096635 

-06,063365 

-loo,  0 

A 0 , 006P&U 

139,6AA.1<36. 

13,173*21 

125,909209 

5, 502*75 

•00,577125 

rj?c  ,  0 

7 1 , 9^2552 

154,682614 

19,033963 

130,196560 

12,243017 

•73,916983 

4CC.0 

*1,456672 

lAfi, 662194 

27,515945 

135,453062 

19,924999 

•66,235001 

7  o.: .  c; 

''•9.64  7  2R4 

lPl,fl49fl53- 

36,l'fll5i3 

141,149044 

20,490566 

•57,669432 

coc ,  0 

-  ■  9  6 , 6  m  2  4  £  fi 

-  194,29.;  042 . 

....  45,4(16220- 

. — -  14?,t)?l74a  - 

-  37.31 *275 - 

-  -48,344725  - 

9  r,  g  ,  r, 

102,640267 

2  0  6 , 0  3  G  6  7  5 

55,379426 

152,932364 

47,706480 

*36,371520 

moo, 3 

in? 1 716703  . 

217,115036 

65,9u4297 

156,601665 

58,31 3352  ■ 

•27,046646 

1230,0 

115,730224 

237, 50 '1522 

rt6,29;:742 

170.250692 

00,699796 

•5,460204 

1 4  0  o  ,  o 

121.801.1*7 

255 , 01 41 fi  A 

112,067033 

161,108400 

104,476087 

18,316087 

1*30,0 

127,0367^7 

272,423437 

•  136 , 955569 

191,570546 

179,364623 

43,204023 

ifloo.u 

132,541001 

—*87,695969 

-162,899659 

- 201,413350 - 

-155.308713 - 

-  -  *9,146713 

2000,0 

139, 4 JP665 

301,909628 

190, <’63531 

210,763336 

182,472585 

90,312585 
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GROUP  ACmiVlTT  7MEBMOCHEHICAL  properties  PROGRAM 

- - -CASE — HUMBER _ 3 _ — _ 

■  MOLECULAR  S1RUCTUR6 -  - - - 1  - 


- K-.H-  M..JH - 

•  •  -  •  0 
M-C-C-C-C-C-N 
m  ••  «..• _  .. 

H  M  H  H 


ATOMIC  COMPOSITION  of  MOLECULE 

ATOM  '  “  NUMBER 

H  10 

- c - a - 

o  i 


group 

NUMBER 

_ 


i 


3 


4 


3 


GRCUP  STRUCTURE  AND  COMPOSITION 

“group"  chFmTcal Ftp pci  grid ~r  oiTYFT  F"c c  lTmn  bcnd  bond  Fpoup  "number 

WEIGHT  OF  COMPONENT  COORDINATE  COORDINATE  VECTOR  TYPE  OF  CORE  ATOM 

.  C  *  3  ”  12  3  1  2 

H  3  8  7  10 


ATOMIC  COMPOSITION  OF  CURE  ATOM  AND  LIGANDS  IN  EACH  GROUP 
GROUP  nukrer  slbqhoup  h  c  c  n 

1  COMt  0  1  0  o 

LIGAND  7  4  1  0 

- - . - - -  ...  . .  LIGAND - 1 -  o  o  o  —  -  - - - - 

LIGAND  1  o  o  o 

Lir.AND  10  0  0 

LIGAND  3310 

LIGAND  3  1  0  C 

-  .  . .  ■  LIGAND  •  !-•-  -  0  -  -  0  -  0  -  . .  . . 

LIGAND  1  0  0  C 

3  COME  0  1  0  C 

LIGAND  3200 

ligand  1  CO  0 

. — . . .  -  -  LIGAND  1 .  0  0  0  •  - - . 

4  CORE  0  1  0  0 

LIGAND  1  1  1  0 

- LIGAND- - 7- - 3 - 0 - 0 - 

LIGAND  1  0  0  0 

LIGAND  1000 

. -  .  -  _ -  5  -  CORE  0  -1  1 -  0 -  — . -  . . — . 

LIGAND  9  4  0  0 

LIGAND  1000 


--  -■  NUMUER-OF -CORE  ATOMS  IN  LONGEST  CHAIN  -■  5  -  - - - - 

GROUP  NUMBER  CF  CORE  CONSTITUENTS  CT  LONGEST  CHAIN  •  1,  2,  3,  4*  3# 
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. - .  NONGHCUP- . INTERACTION  AND  ...  INTERNAL  ROTATIONAL  SYMMETRY  CONTRIBUTIONS 

TYPE  RJNG  NUMBER  cr  amount  heat  of  FORMATION  entropy  neat  capacity  coefficients 

. .  - . . GROUP --NUW.ER - - - XCAU - CAL/Dbr,  K . GAl ^DE6— K — WUBFM***  GAfc/C6S-**“5 -«**4 

internal  rotation  1  1  *2,103 


EXTERNAL  rotational  symmethy  ccntr 1 90  t ion 
source  symmetry  dumber  entrofy 

CCMPLTEC  1  0,00000 


. . .  ' . .  -  -  --  OPTICAL-  ISOMER-- •  CONTRIBUTION -  - — . - . . . 

SOURCE  ARYMMLTRIC  PSEUDCASYMM6TKIC  ENANTIOMERS  PESO  TOTAL  OPTICAL  ISOMERS  ENTROPY 

COMPUTED  0  .  s|  0  0  0  0.00000 


TNERMCHHEMICAL  PROPERTIES  of  kglecule 
.HEAT  CAPACITY  COEFFICIENTS 

CAL/DEG-MOLR  CAL/CE'i*»Z-NULE  C  AL/OEG*  *3-POLfc*  C  AL/DbG**4-rCLfi 
4,313923728*00  1 ,  C1945774ET 1  ®5 , 28l2u52E-C!3  1 ,  O4Q7425OE-C0 


T 

C(p) 

s 

H  (  t  >«m  0  > 

•  tQIT ) pH 1 29S  > ) /T 

NITI-MI 2901 

pELNF<298)*MT).H<29ai 

TJl-6  K  - 

-  CAi./nfcfi.Mou  - 

-  CAL/DFfi-PHLE 

-  KCAL/MOLE 

— CAL/DLG.MCLE 

-  KCAUHOlfc 

MCAL/PCLF 

290,0 

30,279255 

91,5*1677(1 

5t36679fl 

91,556770 

0,000000 

•54,500000 

3oo ,  0 

30,425566 

91 ,  79  9797 

5,427393 

91,557447 

-,  060705 

■54.519295 

<  n  C  .  0 

37,303370 

101,515347 

8,621160 

92,079391 

3,454382 

•51,125610 

?  0  o ,  o 

43,3*4710 

U0,50?62l 

12,862295 

95,516627 

7 ,495497 

•47,084503 

4  o  c  #  c 

48,  717030 

118,095745 

17,473322 

96  ,7?2?C5 

12,106524 

.42.473470 

700 , 0 

53 , 367776 

126,767110 

22,532962 

102,172562 

17,216184 

•37,303610 

F  CO  » 0 

r>7 , 3993V  J 

134 ,103108 

23,126239 

105,713907 

22,759441 

•31,820559 

9  G  G  ,  0 

*0,674319 

141,159566 

34,04431-4 

109,265671 

28,677906 

•25,902494 

1000  ,  C 

63,855007  . 

147,701647  .  -- 

40,294628  ' 

—  112,78381 7 

34,917830 

■  ■19,602170 

1200,0 

*6,582437 

199,704193 

53,553096 

119,628945 

40,1*6298 

•6,393702 

1400,  C 

72,0*4235 

170,631040 

67 , 636161 

126,1*3724 

62,269303 

7 ,6  09  363 

1600,0 

74,656959 

100  ,  443684 

82,338252 

132,336325 

76,971454 

22.391454 

lflOO  ,0 

77,401166 

109,407257 

97,563710 

138,106750 

92,196912 

37,616912 

2000,0 

00,216410 

197,703424 

113,316788 

143,720429 

107,949909 

53,369989 

ORCUF  ADDIT  I V  t  TY  TMERHOCHEM ICAi.  PROPERTY*  PROGRAM 
- - - CASE  _  .NUMBER _ 4 _ 1 - 


MOIBCULAR-  STRUCTURE  - 


- -  H  M- 


-  c  --  - C  -• 

•  *  •  • 

H..c-C-«*-0— C-C—H 

- ..  .  -« — ~m~  — 

C  C 

H  .  H 


H  H  H  H  H  H 


ATOMIC  COMPOSITION 
ATOM 


— C  - 

0 


or  molecule 

nuhbgr  " 
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ORCUP  STRUCTURE  AND  COMPOSITION 


CROUP 

NUMBER 

-  -  1 


GROUP 

CHEMICAL  STMROl 

uft{p  ROM 

GRID  CCLLMH 

BCtvD 

BOND 

CROUP  NUMBER 

HEIGHT 

Cf  COMPONENT 

coordinate 

COORDINATE 

VECTOR 

T  TPE 

or  CORE  ATOM 

-  26* 

c  -  • 

.  6 

.  -  A  _ 

.  . 

-  1  —  _ _ _ _  .... 

c 

6 

6 

3 

1 

2 

H 

3 

1 

8 

1 

0 

9  - 

4  _ 

_  g  -  -  . 

b 

6 

1 

7 

1 

0 

640 

C 

6 

6 

-  2  - . -  .  -  ■  ... 

C 

6 

4 

7 

1 

1 

C 

4 

_  _  (,  . 

6 

1 

_ J._ 

1 

3 

C 

8 

"  6 

5 

1 

9 

296 

C 

4 

6 

■ . -  3 

c 

6 

6 

5 

1 

2 

H 

1 

3 

e 

1 

0 

b 

1 

9 

2 

I 

0 

412 

c 

6 

10  - 

. .  4  ...  •  1 . . 

c 

6 

6 

7 

1 

2 

c 

6 

14 

3 

1 

6 

296 

c 

8 

6 

9 

c 

6 

4 

1 

1 

2 

b 

11 

9  ■ 

4  ...  - 

1 

o  . - 

h 

11 

6 

5 

1 

0 

H 

11 

3 

6 

1 

0 

6'C 

C 

6 

14 

6 

D 

6 

10 

7 

1 

4 

C 

4 

14 

1 

1 

7  . .  . 

c 

6 

16 

3 

1 

e 

c 

8 

14 

5 

1 

9 

~  29  6 

.  "  C . 

4 

14 

7 

c 

6 

14 

5 

1 

6 

h 

1 

11 

0 

1 

0  .  ■  »■ 

b 

1 

14 

1 

1 

0 

H 

1 

1? 

2 

1 

0 

296 

c 

6 

"  16 

. 

6  . 

c 

A 

14 

7 

1 

6 

h 

3 

19 

2 

1 

0 

b 

6 

19 

3 

1 

0 

b 

9 

19 

4 

1 

0 

~?9* 

'“c” 

8 

14 

9 

c 

6 

14 

1 

1 

6 

h 

U 

IT  ■ 

4  . 

1 

0 

H 

U 

14 

5 

1 

0 

b 

11 

11 

6 

1 

0 
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COMPLEX  LIGANDS  7  AND  6  OF  COHE  ATCH 


HL  .  ■ 

ARE  EQUAL | 


COMPLEX  LIGANDS  7  AND  9  OF  COHE  I  TCP.  6  ARE  EQUAL  * 


COMPLEX-  UCANCS 2 — AHO - QF-COKE^T-t* - 4 — AfcS--fiQLAW- 


CONp.SX  LIGANDS  1  AND  3  OF  COKE  ATOM  2  “RE  EQUAL t 


COMPLEX  LIGANDS  1 

AND  5  OF 

COME 

ATCH  2 

ARE  EQUAL, 

-AT-CXIC  GOMPGSITfOK  -~4f- 

-CORE  ATOM- 

-A  NO  — 

trKANDR- 

IN 

GROUP  NUMBER 

SURGHOUP 

H 

c 

c 

K 

1 

CQwfc 

0 

1 

c 

0 

LIGAND 

ID 

7 

1 

0 

LIGAND 

1 

0 

G 

0 

L  i  G *NP 

1 

0 

0 

2 

COHfc 

0 

1 

0 

0 

LIGAND 

s 

1 

G 

0 

LIGAND 

3 

1 

0 

0 

LIGAND 

9 

A 

1 

0 

3 

CORE. 

0 

1 

0 

0 

LIGAND 

ID 

7 

1 

c 

ligand 

1 

0 

c 

0 

LIGAND 

1 

0 

0 

0 

LIGAND 

1 

0 

0 

0 

4 

CORE 

0 

0 

i 

0 

LIGAND 

V 

4 

c 

c 

LIGAND 

9 

4 

0 

0 

5 

COKfc 

0 

1 

£ 

0 

LIGAND 

ID 

7 

1 

0 

LIGAND 

1 

0 

0 

0 

. I...  . 

-  -  LIGAND  - 

-1 

P 

0 

-0  - 

LIGAND 

1 

0 

0 

0 

6 

CORE 

0 

1 

s 

•  0 

LIGAND 

9 

4 

1 

0 

L  1  GAM> 

3 

1 

0 

0 

LIGAND 

3 

1 

G 

0 

... 

- LIGAND 

3 

-  -  -  1  -  ■ 

-  Q 

0  .  .  - 

7 

COHfc 

0 

i 

o 

0 

LIGAND 

ID 

7 

l 

0 

ligand  . 

1 

0 

L 

0 

LIGAND 

1 

0 

0 

0 

LIGAND 

1 

0 

o 

Q 

.  -  In  -  -  — 

CG«c  - 

■  ■  0 

i 

0 

G  •  -  •>  - 

ligand 

ID 

7 

1 

o 

LIGAND 

1 

0 

C 

0 

LIGAND 

l 

c 

0 

0 

ItGANfi 

1 

0 

0 

'  c 

9 

CURE 

D 

1 

0 

0 

-  -  LIGAND - 

AD 

7  — 

1 

...  -0  ■  • 

LIGAND 

1 

0 

Q 

c 

LIGAND 

1 

0 

Q 

•  0 

LIGAND 

1 

0 

0 

0 

NUMNfcR  OF  LORE  ATOMS  IN  LONGEST  CHAIN  *  *> 

HOUR  NUMBfcH  OF  CORE  CONSTITUENTS  CF  LCNGEST  CHAIN  «  1,  2,  4>  6,  7, 
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TYPE. 


NONORCUP  INTERACTION  AND  INTERNAL  ROTATIONAL  lYRMETRY  OONTRliUT ION* 


QAUCH6  EMER 
gauche  etuer 
DIURTIAPY  cTHER 
INTERNAL  ROTATION 
internal  POTATION 
INTERNAL  ROTATION 
INTERNAL  ROTATION 
INTEHNAL  ROTATION- 
INTERNAL  ROTATION 
INTEHNAL  ROTATION 
INTEHNAL  ROTATION 


.  — RING  NUMBER  .CH. 
GROUP  NUMBER 


1 

2 

3 

5 

— A - 

7 

0 

9 


-..amount. ~.heat  or  formation..  entropi! _  _ heat  capacity  epirriciCNTi . .  . 

KCAL  CAL/DEQ  K  CAL/OEG  %  CAl/OtG  K*iB  CAL/DEO  MU  ClL/069  **«4 

2  ,600 

- a . — taoo-* - 

1  0,400 

1  -2,1*3 

1  -2 1 

1  -2,103 

- 1 - - - -Brl*3 - - 

1  -2,183 

1  -2,163 

1  -2,163  -  --  - 


EXTERNAL  ROTATIONAL  SYMMETRY  CCNTRIBUTIOA 

SOURCE  -  SYMMETRY  NUMBER _ ENTROPY - - 

CCHPUTEC  2  -1,37746 


CP  T I  CAL  JSQMER  CONTRIBUTION 


SOURCE 

ASYMMETRIC 

PSELDCASyn METRIC 

enantiomers 

RESO  TOTAL  OPTICAL  ISOMER* 

ENTROPY 

COHruTEE 

0 

0 

0 

0  0  ■ 

0,00000 

THfcPHCCHEMICAL 

PROPERTIES  OP 

MOLECULE...  .  .  - 

HEAT  CAPACITY  COEFFICIENTS 


. 

. .  CAL/Dtq*MCLF  -  CAL/CtG*  *2 -MOLE-  •  - 

-  a , 7  0  221 192E  *  DC  2 ,  470ft24l7£-Cl 

-CAL/DtG*O«K0Lt-  ■ 
*1 , 06496026  E  *04 

— 4ALVCEG*M.*ClrS 
3  ,69O5793lE«O0 

T 

C(P) 

S 

M(T)-HO) 

-(0(T)«H(29f  )  )/T  H(T>-M129BI 

DilNF(295>*MTM 

PE<’  K 

CAL/PER-MOLE 

CAL/LEG-fOLE 

kcal/holf 

CAL/CEG-HCLE 

KCAL/MCLE 

KCAL/KCLE 

290,0  . 

.  ....  43,919654  _ _ 

...  105,176696  ... 

_  6,549824 

-105.176696 — . 

-  0,000000  -  - 

.  -17,780000 

330 , 0 

49,222327 

105  ,5C493i 

6,647^6^ 

105,177791 

♦098142 

*67 ,681858 

400.0 

53,053.161 

121 , 1 1  9 1 64 

12,202524 

107.2P7412 

5,7.12701 

•62,047299 

'j  3  0  ,  o 

74,566214 

136,96565u  . 

19,181431 

.  -  111,702633 

12,631606 

•75,148392 

430,0 

P4,l?0320 

151,432697 

27, 13 ''742 

117,131166 

20,560919 

-67,199061 

70C.0 

«2, C6311* 

1A5.0H773 

35,95185? 

123,011073 

29,402029 

-36,377971 

.-.acc,o  - 

-  03 , 7SC031  --- 

—177,756622 - 

- 45,5015.-2- 

- -129,067024  — 

- 3**951670 - 

- *46, *28321 

?oo,o 

104,534505 

169,726798 

55,671777 

135,148850 

49, 171953 

-36,658047 

1300,c 

109 , 769972 

201,016644 

66,39M:3 

141,176564 

59,B402B0 

•27,939720 

1200,0 

120,007621 

221.925792 

89,357338 

152,919489 

82,807313 

•4,972489 

1*01, i 

132,2«0457 

241,302094 

114,529497 

164,174756 

107,979674 

20  < 199674 

1  *  0  C  ,  0 

149,445957 
_ 124  ,4di,599 _ 

260,016097 
_ 37E^  965995 _ - 

142,398365 
_  ^74 ,021223 

174,987759 

- 185*4018*4— 

195,021500 

136,046542 
lit  ,271490 

46 i 266342 
•0,491409 

2  JOC  ,0 

209,664662 

299,057011 

213,020846 

206,471023 

110,691023 
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GROUP  ADDtTJVm  TmERMOCHEMICAL  PROPERTIES  /TOPOGRAM 

- CASE — NUMBER- - 9 - 


hOLtCULAR  -STRUCTURE 


- - - - - - - H— H~H - 

•  -  m 

H-C-C-C-CK 
H  H  H 


ATOMIC  COMPOSITION  Of  MOlECILE 

ATOM  "  NUMBER 
M  7 

N  l 


ORCUP  STRUCTURE  anc  composition 

.  GRriJp*'  GROUP  ~e  H  E  M  c7l~S~V  K IJOL  UR  I  fl  ROW  GRID  COLUMN^  QoTiB_”TonO 

NUMBER  HEIGHT  of  COMPONENT  COORDINATE  COORDINATE  VECTOR  TVPF 

—  •  i ...  -  2*6  -  ■  Q  .  ..  .  3  ....  U 

C  3  13  3  1 

¥  3  9  7  1 

H  5  11  5  1 

C  3  13  2 

C  3  11  7  1  .  1 

C  3  15  3  13 

¥  *  13  5  i  0 

3  679  C  .  .  ..  3  .15 

C  3  13 

¥  1-15 

CN-  — . 3  - .  17 

H  5  15 


ATOMIC--  COMPOSITION  l)F  CURE  .  ATOM  -  -AND— t  KaNM  I  fcACH-  -  GROUP - 

GROUP  NUMBER  SUBGROUP  H  c 

1  CORc  0  -  1 

LIGAND  5  3 

LIGAND  1  0 

•  LIGAivD  1 - 0 

LIGAND  1  d 

2  CORE  0  l 

LIGAND  2  2 

LIGAND  3  1 

LIGANP  1  o 

ligand  1  0 

3  CORE  0  1 

LIGAND  .  5  2 

LIGAND  1  0 

LIGAND  0  1 

LIGAND  1  D 


NUMEER  OF  CORE  ATOMS  IN. LONGEST  CHAIN  ■  3 

GHOUP  NUMBER  OF  CORE  CONSTITUENTS  CP  LONGEST  CHAIN  ■  1,  2,  3, 


C  N 

-  0  C  • 

0  1 

0  0 

0  0 

0  0 

0  1 

C  0 

0  0 

G  O 

0  G 

0  0 

0  0 

0  1 

0  0 


GROUP  NUMBER 
OF  COPE  ATOM 
1  - 
2 
0 

0 
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■  ■*  '> 

GROUP  ADDITIVITY  THERMOCMEHlC*t  PROPER^fls  PROGRAM 
- i - CASE — NUMfi£A - 1 - 1 - : - 


- -MQUfcCUUAR  - -STRUCTURE--- 


- k - N - - - 7 - 

H  .  C*  H 

H--C  C*-H 

KkHHHHH 

Hp-C— C— C— C--C—  C--C—  C  C--H 


- k  —H — H--H  • -H-  —  H.  Q  M 

*  • 

H  M 


ATOMIC  -  -COMPOSITION  -  OF  .  MOLECULE 
ATOM  KUHHER 

.  c  u 


ORCUP 

STRUCTURE  AND  COMPOSITION 

GROUP 

GROUP 

CHE  H 1  PAL  ST  POOL 

GRID  ROW 

GRID  CCLUHN 

AC  NO 

BOND 

GROUP  NUMBER 

NUMBER 

WEIGHT 

OF  COMPONENT 

COORDINATE 

coordinate  vector 

TYPE 

OF  COPE  ATOM 

i 

296 

c 

.  a 

. 37 

1 

c 

e 

40 

3 

1 

2 

h 

a 

34 

7 

a  ...  . 

1 

D 

.  h 

10 

39 

6 

1 

0 

2 

399  . 

C 

a 

40 

2 

C 

a 

37 

7 

1 

1 

C 

a 

43 

3 

1 

3 

.  —  - 

.  .  -  -H  - - 

- -  6 - 

-  -3P  -•  -  • 

.* - 

-  1  - 

. - . 0  •  -  - 

V 

10 

3  a 

6 

1 

0 

3 

399 

C 

a 

43 

3 

c 

a 

40 

7 

1 

2 

c 

H 

46 

4i . . 

3 

a-- 

1 

4 

K 

10 

.  41 

6 

1 

0 

4 

399 

C 

0 

46 

4 

C 

0 

43 

7 

1 

3 

c 

a 

49 

3 

1 

5 

_  ....  _ 

-  -  ►* . . 

6 _ 

-  — -  -44 - 

.  .5  . 

..  .1-.. 

_ _ _ _.,Q  _ 

V- 

10 

4  4 

ft 

1 

0 

5 

359 

c 

a 

49 

3 

c 

a 

46 

7 

•  1 

4 

c 

a 

5? 

3 

1 

6 

H 

6 

47 

b 

1 

0 

b 

1U 

47 

ft 

1 

0 

6 

399 

C 

a 

52 

6 

C 

a 

49 

7 

1 

5 

c 

a 

55 

3 

1 

7 

.  H . 

.  .  6 .... 

50  - 

0  . 

1 

...  0 

I- 

10 

-  50 

ft 

1 

0 

7 

399 

C 

« 

59 

7 

C 

a 

52 

7 

1 

b 

C 

a 

59 

3 

1 

a 

..  . 

-  .  .  h  . 

..  6  .  .. 

.  .  ...  53 

ti 

■  1 

0 

b 

10 

53 

ft 

1 

0 

3 

5Q2 

C 

0 

59  . 

6 

C 

10 

61 

4 

1 

10 

c 

b 

59 

1 

1 

9 

.... 

_ 

_ _  c  _ 

..  ...  a _ 

-  .  55  . - 

7 

—  -  1 

.  7 

H 

10 

57 

6 

.1 

0 

-  9 

399  .. 

C 

.  ...  5  . 

59 

9 

C 

a 

59 

'  5 

1 

a 

c 

j 

. . . .  5  _ 

61 

....  ...bis 

2 

..  ,_,7. _ 

1 

■  ~  !A- 

u 

_  _  {j  _ 

h 

3 

57 

a 

1 

Q 

'  13 

399 

C 

10 

61 

10 

C 

a 

63 

2 

1 

12 

c 

a 

59 

1 

0 

H  . 

.  . . . 12  -  - 

-  63 

4  . 

-----  1 

.  0 

b 

12 

59 

'  6 

1 

0 

.11 

.  399  . 

.  c .  . 

3  . 

61 

11 

c 

5 

59 

ft 

1 

9 

c 

a 

63 

4 

1 

13 

.  •  k . 

.  ..  ..j  . . 

. . 59.-  _ 

...  a  . 

_ Q  .  . 

b 

1 

■63 

. 

2 

*  — * 

.1 

0 

.  12 

399  - 

C 

a, 

.  63 

12 

C  . 

5 

63 

1 

1 

13 

C 

_  y. 

10 

. . a _ 

61 

_ _  6 _ _ _ 

6 

_ 1 _ 

1 

10 

h 

10 

65 

4 

1 

0 

13 

399  - 

c 

5.  ■ 

63 

13 

c 

3 

61 

8 

> 

11 

C 

j.  - 

B 

-J 

63 

.  ~4  5 

5 

(1 

_ 

12 

- Q- - 

D  1  fur 

h 

5 

66 

3 

1 

0 

— . -M  J  NU  ’ 

CROUP 

■  4U4iHER~ — 

MUMMERS 

OF  RING  CONSTITUENTS  ■ 

e»  .  9,  n, 

13, 

12 

.  10i 
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)  . 

NONbKCLP 

INTERACTION 

AND  INTERNAL 

ROTATIONAL 

Sv^'try 

contributions 

TYPE . . . 

— RltlG  fcUKBEH  CR.. 
GPOUR  NUMbER 

^JLfcCLNV.  hEA^ 

cp  pdrha-T.ion . 

►  UL  . 

.  ENTROPY - 

CAL/OfcG  K 

CAL/DEC 

-HEAT,  CaP AC I TY_cO£FP IC I  ENTS 

K  CAL/CEG  K**2  CAL/CEG  M*3 

CAL/CEG  H • • 4 

RING  CORRECTION 

1 

1 

,000 

18,000 

•1 *3120E« 

Oi  2.5547E-02  •1,133?E»05 

1.0331E*P5 

GAUCHE  ALKANE  . 

INTERNAL  ROTATION 

-  8  7 

1 

-  -1 

1 

,boo 

-2,103 

external  rotational  symmetry  ccntr  irut  ion 

Sm;MCE  SYMMETRY  NUMBER  ENTROPY 

ceMPUTfr  1  oiooooo 


.  OPTICAL  ISOMER  CCNTHIBUTICN 

Sfil'NPE  ASYMMETRIC  pSELDCASVPhETRlC  ENANTIOMERS  Mf;SC  TCTal  OPTICAL  ISOmBRS  FNTRCPY 

•CCMfU-Tfit . . -~C- - . ~-ii - ‘ — -q -  -0 - -  — 0 - 4tOOOOO 


TPERpCCHEMJCAL  PROPERTIES  OP  hCLfOlLE. 

. .  .  .  -■  _ HUAT-  CAPAC1  TYL.COEPP  10 IEN7S  . . . .  J  . . 

CAL/DIG-MOLE  CAL/t£n**2“MULE  CAL/CfcO*  .  3-MOLfc  C  AL/CEG*  M-PCLE 
-1 , ?a?103P5E*01  3,28ia64JlF-01  ,  •! ♦  8923574Pfc»C  4  4 , 14826526E-0B 


T 

ctni 

c 

h<t  >-««;■> 

-{GIT)-Hf29fc))/T 

PIT  >-H(?9e > 

DELHP(298  } *M  T )*H(  298  ) 

_ CtC  k - 

— CAL/LLt»MCuE — 

- CAU/LF-G-Mntt  — 

- KCAL/HOLP  ... 

- CAU/ni R-MCLfc - 

—  KCAL/HOLE - 

■  —  - . kcal/mcie  . 

29»’(  o 

64, 271413 

138,576770. 

7 , fc  7  JV65 

130,576770 

0  ,000000 

-65.63C0CP 

3  30,0 

ft* ,7/3711 

I39,00?i92 

7,002960 

130,578200 

,120995 

•65,501005 

«.K,o . 

f*  ,83071(1 

160,5Pt:412 

15,354961  . 

141.3P59/2 

7 ,60099ft 

•57,949004 

b  «.  {•  t  j 

li'4,i4057N 

API.  75‘270 

24,876131 

147,353946 

1  7  ,  2C2166 

•  4  A ,427634 

* ;  t- ,  g 

11 5, 9/621 

202.17-1 -? 

36  t  III  “O'!  P 

154,001376 

~8 , 4?6r)37 

•37,203963 

.....  7  Z  u  ,  0  -  - 

..  -l33,UiS03-  - 

- 221 ,7003vg  - 

-46,704995  — 

—--162,975355 - 

-  -4i, n nun — 

. *34, 5\69  7  n - 

nv.D 

144,8‘'N3b? 

24C  ,2R:*5o 

#>2,714421 

171 ,407979 

55,040456 

•10,589544 

<00,0 

J.T-4  ,  VifcAM 

257,923717 

n  |  6  9  6  4^1 

180,12092? 

.70,022516 

4,392516 

.10  00,  u 

162,612304 

274,634590 

93,564264 

188 ,744?C9 

85,690299 

20,260299 

1?3G,  0 

305,457310 

127  ,<13K  .*8 

205,674116 

119,739843 

54,109043 

1 4  r,  c ,  o 

1o4,*>A13C4 

333,194977 

163,425566 

221,943034 

155,751601 

90,131601 

UflL.U  ... 

A'72, 74660b  — 

358,376721 - 

-  201,16'  ?84  - 

.  237,447771  - 

193,486319 

. — 127 ,056319  -  ■ 

1POO,  o 

201,717551 

301,579466 

240,576940 

252,108924 

232,902975 

167,272975 

2010,0 

213,470059 

4 Q 3 ,407665 

202,^32747 

266 , 2P6274 

<74,350782 

208,728702 
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GROUP  ACC  1  T I V I TY  THERMOCHEMICAL  PROPERTIES  PROGRAM 


- CASE — NUMBER - 7 - 

- - -HOUCUUAR- -STRUCTURE— 


ATOMIC  COMPOSITION  OF  MOLECULE 


ATOM 

H 

C 

0 


NUMBER 

6 

6 

1 


GPcir  structure  anc  composition 


CROUP 

GROUP 

CHEMICAL  SVMPCU 

CAR  ID  ROM 

GRID  COLUMN 

NUMBER 

MU  GMT 

CF  COMPONENT 

COORDINATE 

coordinate 

■  1  — 

-m* 

*  c 
c 

- 7  ~  . 

3 

10 

-  3  - . 

3 

3 

M 

.  7 

1 

8C*'C  BCNO  group  number 
VECTOR  TYPE  Of  CORE  ATOM 

2  1  2 

P  2  3 

7  .  i  .  0 


A  Ufcft  C 

c 

h 


i 

5 

1 


3  3  016 


6  1016 


7 


1016  C  10  7 

. . C -  -  12.  - 9 - .6 

C  7  7  1 

H  10  9  3 


- Z 

.  1 
1 


RING  NUMBER  1  IS  6  MEM.FKEU 

_r,PCUP.  NUMBERS  OF  RING  CONSTITUENTS  •  2*  Si-.  7,..  .  6,  3,  it 
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NOMiftCUP — INTERACTION-  -AM) — INTERNAL- • RGTAt  I  ONAl.— 6  YMMfeTR  v  •-  CONTRIBUTIONS  -  .  . 


TYPE 

JKUHNAW  HQTaT  I  DN 


R  | ri G  NUMOL^  CR  AMOUNT  MEAT  i)F  FORMAT  ION  ENTROPY 
UKQUP  MJMHfcR  KCAL  CAl/ObG  K 

2  1  -1.477 


MEAT  CAPACITY  COEFFICIENTS 
CAL/r;£G  K  CAL/DEG  K**2  CAWCEG  K •  * 3  CAL/CEG  M*« 


EYTERNAL  ROTATIONAL  SYMMETRY  CCNTFI6UT10A 

SCLPCE  SYMMETRY  NUMBER  EMRQFY 

CCKFLTEC  1  0,00000 


ECt'l-.CE 
-CorruTtc ■ 


optical,  isomer  contribution 

AFYMMI  TMC  PSELT.CASYMKETMIC  ENANTIOMERS  f'ESC  TOTAL  OF'TIGAI  ISOMERS  FNTRCPY 

c  «  _  O  .  .  f.  .  . .  0  --  ■  - .  0,00008 


ThERMCCMEM  ICAL  PHOPENTIES  CF  MCLFCULE 

.  .  ...HEAT  CAPACITY  .COEFFICIENTS  -  - - 

CAL/DEG-MOLE  CAL/CEO* ♦2-POLE  CAL/DEG* *3-hCl E  C AL/Cb^ ••4-kCLE 

-p  ,  09  ac  1391E*0L»  1. J«374651E*C1  *•!  ,C2Q6i739E-04  2  , *3403300  4E-O0 


T 

f.  i  P  > 

c 

H  (  T  )  •  H  (  0  ) 

r.M-  k  ••• 

*  CAt  />: EG-hCLE 

- CAL/CFG-KOL-P 

-  KCAL/Mt)LE 

<  9  *- ,  0 

24,2558ot 

75 ,172540 

?,  RO  >581 

Z  c  c  ,  o 

?o , 4?2?14 

75,335343 

2,h49?5V 

4  1:0,0 

31, 998180 

83.432052 

5,682164 

SuU.O 

36,225790 

91.246346 

9  1 2  U  4  ?  8  3 

t  oc  ,0 

43,310666 

98,702449 

13,29  ’5 7  1 

701., 0  — 

.47,399  2W 

-  10S  ,  69743;! 

1 7 ,833565 

F  OC ,  0 

SC  ,676464 

112.24971? 

22,743385 

900 , 0 

!»3, 319070 

HE  .375497 

27,947733 

icoc .  o 

55,502199 

1?4, 109472 

33,391891 

l  z  r.  «, ,  o 

*•9,193967 

134,561641 

*4,868682 

1  4  0  fi ,  0 

63,158546 

143.971673 

57 | 0B7465 

It $0  ,0 

68,0O249L 

- • 152.747297  . 

.  7U| 244058 

IF 30, 0- 

77,53237? 

161,317942 

84,6l439i 

2000, 0 

90,734736 

170,125540 

101|5565!.16 

.(GU  )-K(29P  )  )/T  MU-k<?98)  OELRF  (  298 )  *M  T  ) -h  1 298  » 

C*WPf  f-MClit -  -  KCAL/MOLE  . .  MCaWMCLE  - 

75,1725*0  O.GCOOQO  *23,200000 

75,173043  ,  0*8*78  -2  3,15132? 

76*228*94  2.6P138  3  -20, 318417 

7tt,4^894l  6,403702  -16,796290 

01,219152  10,4^9990  -12,710010 

R4,2?l73ft -  li,0429«4-  -  -8,147016  - 

8  7 ,320206  17,942*105  -3,257195 

90,434217  25,147152  1,947152 

93 , 5l £  16  2  30,591310  7,391310 

99,504890  42,048101  18,860101 

105,195164  54,287084  31  ,087084 

U0.U95124 -  47,4,434/7 . —  - 44,243477 

115,754713  62,013511  58,813811 

120,747577  98.755926  75,539926 
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- AtCMC  COKPOSITJOV-  0F--CORE  ATOM  —  AND- --L  [CAM'S — EACH  GROUP 


2 

-  co«t - 

LIGAND 

LIGAND 

L 1 Q  and  - 
ligand 

CURB 

17 

1 

-l-  - 

1 

0 

4 _ 

—  1- — 

8 

0 

0 

0 

1 

0 

0 

-  o . 

0 

0 

— c  — - 

0 

0 

0 

c 

3 

LIGAND 

ligand 

LIGAND  - 
CORE 
LIGAND 

3 

11 

.  . .  3  .  _ 

0 

17 

1 

1 

5 

...  1 

1 

6 

0 

0 

0 

0 

0 

0 

0 

0 

D 

-4 

ligand 

ligand 

C0H6  - 

Ligand 

ligand 

ligand 

i 

l 

0  - 
V 

3 

, 

o  »  wr  ri  r 

: 

0 

0- 

-0  - 

0 

0 

0 

0 

0 

0 

LIGAND 

1 

0 

0 

0 

5 

CORE 

0 

1 

0 

O’ 

--  -  - 

LIGAND 

~V  - 

-8 

-  c 

G  .  ..  .  . . . . . -  . . 

ligand 

i 

0 

0 

0 

LtGAI.r 

i 

0 

0 

o' 

LIGAND 

i 

0 

c 

0 

6 

CORE 

0 

1 

0 

0 

ligand 

i/ 

h 

0 

0 

.  • _ ...  „ 

-LIGAND- ... 

. i 

0  ■ 

0 

c 

LIGAND 

i 

0 

0 

c 

ligand 

i 

0 

0 

0 

.  J. _ _ 

...  .COPE. _ 

_ 0 

. 1  .... 

_ D _ 

ligand 

13 

6 

0 

0 

LIGAND 

5 

2 

0 

0 

ligand 

1 

0 

0 

o  .  _ .  . .  .  . . 

LIGAND 

1 

D 

0  * 

0 

t 

CORE 

0 

1 

0 

c 

LIGAND  - 

.  15  — 

7  - 

.  -  .0 

.. .  c  •  .  ,1  .  — . . . :  .  . . 

LIGAND 

3 

1 

0 

0 

LIGAND 

1 

C 

0 

c 

LIGAND 

1 

0 

.  0 

0 

9 

CORE 

0 

1 

0 

0 

ligand 

17 

8 

0 

0 

.  . 

LIGAND' 

- 1-  . 

.  -0 — 

. 0 ... 

_ o  . . . . —  _ 

LIGAND 

1 

0 

0 

0 

UGNHB 

1 

D 

0 

0 

- - NU/IDEH-  Of  -COR6-4X OHS  -4  N  LONGS  ST--CHAJN _ a - 6 _ — 

GROUP  NUMBER  OF  CORE  CONST  I TUCN  TS  OF  LONGEST  CHAIN  ■  1,  2,  4,  7,  ft,  ,  9, 


- - -NONuRLUP —  INTERACT  I  ON — AND — INTERNAL — ROTATIONAL  -SYMMETRY.  CCMR.J  HUT  J  CNS 


T  Y  pfc 


U/.UCHE  SL*AM 


GAUCHE  A L I* A f . f 
INTERNAL  ROTATION 
IMlhnAL  rotation 
IKTfcHNAL  R  0  T  A  T  JON 
INTERNAL  ROTATION 
[ntfknai  r o t a t i n n 

INTI WWAL  ROTATION 


Mur,  NlJpnfe’P  CP 
GROUP  NUhliER 

£  4 

4  7 

2 

3 

5 

6 
9 


AMOUNT 

4 

1 

- 

1 

1 

1 

1 

1 


HF.AT  OF  FORMATION  EKTHOPY  heat  CAPACITY  COEFFICIENTS 

KCAL  CAL/OkG  K  CAWfcEG  K  CAL/PEG  K*»?  CAL/CEG  K®»3  CAL/C6C  M*4 

3,  £00 
.600 

-2,183 

-2,183 

-2,183 


EXTERNAL  ROTATIONAL  SYMMETRY  CONTRIBUTION 

SOURCE  SYMMETRY  NUMBER  ENTHCFY 

- COMPLTfet- -  - - 1 . . B?OGO*-& - 


CFTICAL  ISOMER  CCNTHIBUTICN 

.  SCUfcCE - AS-YHtftl&iC - fSiLD-C-ASTHKET-R'lC - -E.MANT  JOKERS - _Ji&£0 _ TOTAL  CRJ-UM.  1S0M685 - tURCP-Y- 

CUmFUTEC  1  U  202  1,37746 


GRQUP  M'HHgo  OF  ASYMMETRIC  CARBON  ATQMCST  ■ 


CAL/DEu-MCLE 
-7  ,lb?l5l33E*00 


ThERMCCHfcMlCAL  PROPERTIES  CF  MOLECULE 
HEAT  CAPACITY  COEFFICIENTS 

CAL/eF.O**<£-M«LE  CAL/DtG**3»MC’LH  C  AL /DEO*  •  4-mCLE 

2,3aS»8ll5t.-{'l  .  •l,41303A3Ch-04  3. 31543384E*08 


T 

rcr,  k 

2  JG.f.  .. 

3  0  c  ,  u 
4S0,rj 
5“JC,0 
ft?  0,0 
710,0 

f 00. C  - 
V  0  o ,  c 
isec.o 
iroc, e 
1430,0 
lftoc.o 

-iLOC.O 

soco.o 


C(PT 

CAL/r.EA«KOLE 

- !»C,4r.aia‘2.  . 

yo,3.)2l23 
*5,346391 
77 (960290 
fid  ,  692747 
^7,  7*52667 
1L  ■»,  4  59037- 
111»*°0723 
117,2066711 
xrs ,  7  ft?  Oi> * 

1 ‘5  2, 491577 
139,0^1707 
-  -147,038793 
138,044263 


CAL/CEG-HOLE 

111,99  »>  084  - 

112,336840 
129,001390 
-144,979719  - 
160,16*521 
174,543621 

lfit,U73lH - 

200.921083 
212.597437 
235,16772?  - 

255,074161 
273,167213 

-290r0045lO  - 

3  06 » 0  3  2  7  3  4 


MtT  >-HI .  ) 

kcal/molf 

-  7t,:i4V34..... 
7,116225 

12  |  9  4  3  H  V  9 
20,127024 
28,474909 
37,811956 

-  47,985660  _ 

58,8626(6 

70,329280 

-  94,677178 
120 j  519045 
147,662854 

176,23486* - 

206,679597 

38 


•  ton) 
CAL/r 
— Hi, 

m. 

114, 

-ue, 

124 
130, 
— -134, 
143, 
149 
-  162 
173, 
185 

- 196 

206, 


Ht298))/T 

ec-mcle 

99  808  4 - 

999213 
178978 
753940 
4  0 1 R  9  7 
547876 

003911- . - 

3} 2556 
683001 
117652 
999796 
282263 

*9*44* - 

2004Q3 


hm-w<  ?9a ) 
KCAL/MCLE 
—  -0,000000  — 
,101290 
5,928065 
13,112890 
21,459975 
30,797022 
-  4  0,97  07  26-— 
51,8*7674 
63,314346 
87,662244 
113,5G4111 
140,6*7920 
-169,319926— 
199,664663 


DELMF  (298)TKONha9ej 
xcal/mccf 

- *>97,700000 - 

-57,598710 
-51,771035 
•  44,587116  ■  ■■  - 

-36,240025 
-26 ,9029  7  P 

- - *16,739274 - 

•5,892326 
5,614346 
-  29,962244 

55,804111 
82,947920 

- -111,619926 - 

141,964663 


RING  NIH'HFR  1  IS  6  H^EFRfcU 

GROUP  NUMBERS. Or  RING  CONSTITUENTS.  «  .  2,  8,.  ft*..  6,  _-3#  i, 
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ATOMIC  COMPOS  I T I  OK  QF  CORE  ATOM  AND  UGAKCS  tyfr  EACH  ,QRQUP 

. . . QROL'P  NUKB6.R _ SUBGROUP _  H _ C  -  C  _ N  -L _ 

1  COHE  01  0  0 


SUBGROUP—- 

CORE 

ligand 

-LIGAND _ _ 

__ 

0 

7 

1 

—  -C 

1 

7 

_ g _ 

--  C  ... 

0 

0 

.0 _ 

..  _N  ...J.... 

0 

0 

0 

.  COHE 

0 

1 

0 

0 

LIGAND 

5 

.  5 

0 

0 

.  ligand  - 

-  3 

2 

0  •  , 

0 

NUMBER  HF  CORE  ATOI'S  IN  LONGEST  CHAIN  • 
GRQur  UUMfeER  CF  CORE  COMST ITUENTS-  CF  LONGEST  CHAIN  » 


IS  SYMMETRICAL  ABOUT  ATCM 


JNTEHNAL  ROTATION 


NONGHCtP  INTERACTION  AND  INTERNAL  ROTATIONAL  SYMMETRY  CONTRIBUTIONS 


R I Nli  NUMPfcR  CR  AMOUNT  HEAT  Of  FORMATION  ENTROPY 
G«OUR  NUMBER  KCAL  CAL/DEG  « 

2  1  -1,377 


MEAT  CAPACITY  CCEFflCUNYS 
CAL/DEG  K  CAL/DEG  K •* 2  CAL/CEG  K**3  CAL/CEG  K**4 


external  rotational  symmetry  contribution 


SOURCE 

eCPPl'TEC 


SYMMETRY  NUMWER 
1 


SOURCE 

C?»TL’T£t 


CFTJCAL  isomer  CONTRIBUTION 


ASYMMETRIC  PSEL’OC  ASYMMETRIC 

0  0 


ENANTIOMERS 

0 


total  OPTICAL  isomers 

0 


ENTROPY 

0,00000 


CAL/DECoMOlF 
-7 ,  2Bft997l4E*Q0 


TRERHCCHKMICAL  proper  t I k  s  CP  molecule 
HUAI  CAPACITY  CCEFF  I r if i;ts 

CAl  /CfcG»*i-HULE  CAL/nfcO**3.pnLF  ....  CAUCtG**4-MCLE 

l,SO0:HV7OE-Ol  •! ,  OOlSlRtiPE’-riA  2,52514903^*08 


T 

CIP) 

S 

M<  T  > - H  C -  ) 

H(T)-H(?90) 

DELNF<298>*M(1 >-H(29B) 

iipn  y 

C  AL/PE  ri»HQl.fc 

CAI./CPG-POLE 

KCAL/MOLE 

cal/deg*hcle 

kcal/molf 

MCAL/rCLE 

. —  -35,220000  -  -  _ _  — 

S)C,  n 

?9’,3791Cf> 

62,6?t»  4  75 

3,714612 

02,463193 

,058564 

35i278504 

4 ;  :■ .  o 

33 ,31552b 

«2.37u93B 

7,112121 

03,730707 

3,456093 

38,676093 

5110,0 

45, 044467 

m,7S4is? 

-  11,331135 

.  ■  ■  06,403930  • 

7,675107 

.  42,895107 

6C0.O 

52,1 27311 

110,685711 

16 , 239477 

09,713295 

12,5R3450 

47,803450 

7  *!  0  ,  U 

57,315449 

/,  1  kA'i'O  AQ 

119,123664 

21 , 721.1  )7  • 

27  A7l 1 3, 

93,310037 

_ .97  .045162 _ — 

18,064079 
_ g4  f  0  ISO 95 

53,284079 

59,235095-  - —  . . 

<?co!u 

-«>  i  f  DO\)  r  O  » —  - 

65  ,{113162 

■  — " —  i^/»oo^uoi — 

134 ,521177 

*  7  y  O  lAcO  — 

34  ,1(05763 

100  ,799243 

30,349736 

65,569736 

l  a  o  u. ,  o 

67,325515 

14i,52i5Ufi 

40,6524: 4  . 

104,525130 

36,996376 

72.216376 

1 2  0  C ,  0 

72,134159 

154.2B097& 

54,67.1800 

111,776506  - 

51,014861 

66,234861 

1400,0 

75,69731? 

165,677915 

69,456369 

118,677671 

65,800341 

101,020341 

1600,0 

79 , 72&09U 

176,036999 

84,901656 

125,206962 

il,  324020 

116,544828 

100C ,o 

...  05,4*1599 

-165,734424 - 

•101,458587- 

- 131,399669 - 

- . »7-r*02**9- 

- 133,0225*9 - 

2010,0 

94,054953 

195,149651 

119,346018 

137,304656 

115,689991 

150,909991 

40 


ATCMC  COMPOSITION  OF  CURE  ATOM  AND  LIGANDS  EACH  ,  GROUP 


_fiROU£_MlKB.EB _ 3U  BiiRQU  P _ M_ 

1  CDMfc  0 

LIGAND  i 

.  . . ...  .  ligand  - . 1 

2  CORE  0 

LIGAND  5 

- - ugand - <£— 


CORE 
LtGANp 
..LIGAND  . 

ligand 

COHE 

— LIGAND- 
LIGAND 

ligand 

CORE 


1 

4 

1  ..  _ 

0 

0 

-4 — 

0 

0 

1 


LIGAND 

CORE 

LIGAND 

LIGAND 

LIGAND 

CORE 

-LIGAND- 

LIGAND 

LIGAND 


..a 

o 

0 

-  C— ■ 

0 

0 

-« — 

0 

0 


_ _ NRHW-OF  -C<*S.-A*OH& 

CROUP  NLI'BE R  CF  CORF  CONSTITUENTS  Cf  LONGEST  CRAIN  « 


_ _ _ ring - 1 - l  *  --SVMM6  TR  hGAk-^ROGT-A-T-G^ - - - '  “  " _ 

NON-r.RCLP  INTERACTION  AND  INTERNAL  RCTATICNAL  SYMMETRY  CCNTR I  BUT  I ONS 

HEAT  CAPACITY  COEFFICIENTS 
CAL/DEG  V  CAL/DEG  K**2  CAl/CEG  M*3  CAl/CEG  M*4 


EXTERNAL  ROUT.IONAL  SYMMETRY  CCNTR  I8UTI0A  . - 

F.CIRCE  SYMMETRY  NUHbER  ENTROPY 

COMPUTE C  1  0.00000 


OPTICAL  ISOMER  CCNTHIBuTION 

50Pfc.cE  ASYMMETRIC.  RSEL'CCASYMMETHIC  ENANTJCMER5  MESC  ..TOTAL  CPT 1CAL  ISOMERS  ENTROPY 

CCm*’L'T£C  C  C  0  0.0  0,00000 


TYPE 

R  p.T,  NUKCEF  CR 

amount 

hFAt  Op  FORMATION 

ENTROPY 

INTERNAL  ROTATION 

GROUP  NUHHER 

2 

1 

RCAL 

cal/det,  k 
*1,377 

TRERKCChENICAL  PROPERTIES  CF  MCLECILE 
MEAT  CAPACITY  COEFFICIENTS 

C  AL/DEli«HCLE  CAL/CEH***-P°LE  CAL/DEG**3«KGLF  C*L/DEG**4-MCLE 
-o,  121^C2JIE*00  l,47970*m»'Jl  ,  I4b9l  51 OE-D  4 - 4, 509  M  «  A?  t-04 - - - 


T 

C{P) 

$ 

V4| 

-{GiT  >-M| 29A ) )/T 

DELRf  <298)*MT)-H<298) 

pi. r,  k 

CAL /r-Ef- MOLE 

cal/cfg-mole 

kcal/molP 

cal/deg-kclf 

KCAL/MCLE 

KCAL/PCLE 

?«r> ,  a 

25,6*596:3 

76,292540 

2 , 9  a  7  6  9  o 

76,292540 

o.ooonoo 

20,600000 

.-3r.t,o  - 

■  - 

-  76,445037  - 

...  2  i  959264  - 

-  74,293115 - 

- ,  C  5 1 5  7  7 - 

-  20,851577  - 

4  it  0 , 0 

33,933063 

P-5 ,04561? 

5j962*25 

77 ,4SZ2£1 

i, 354035 

23,654935 

br.ii ,  a 

40f|>r7D77 

93,354718 

9,b9«>759 

79,776579 

6,789070 

27,509070 

bOo  ,  0 

45,0.10059 

101,234194 

14,  P26639 

32,7:2612 

11,118949 

31,918949 

7-0,0 

c:  i  i^ficav 

10P . 6*^35 

ie, 

85 , 8P4837 

15,930599 

36,730599 

a  3 ,  g 

53,6 '12167 

11* ,5F 1946 

?4, *1367^9 

09 . 1AQC72 

21,131099 

41,931099 

-aJb,S7«*Al3 

-  122, 0772113 - 

29,556277 

. 92,467742 - 

- 26, *44887 - 

- 47 ,440567 - - - 

i  c  n  C ,  o 

59 , 0  *”•  0  8  6  i> 

120,10*696 

35x33994f> 

95,736441 

32, 4. 3  2256 

53,232256 

12'iO,  a 

63,640632 

139,339607 

47,607073 

102 , 0p94£4 

44,700104 

65 , 5001  84 

14CC.0 

69,1*5944 

149,54  .773 

60,857234 

108,146241 

57,949545 

7fl  ,749545 

lAO'G ,  q 

77,^1279 

159,264606 

75,439586 

113,932371 

72,531896 

93,331896 

1GC° ,  0 

<<9,641115 

169,026767 

92,043379 

119,508940 

09,135690 

109,935690 

203C,G  - 

..  -iC7r90993ti . 

•179,3*3(114 - 

111,693964- 

- 334  ,970477 - 

- 4JV8t7*4£7-4  — 

- 129T*M274 - 
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GROUP  ADDITIVITY  TMftRHOCMEHlCAU  MDHRTIEI  PROOPAH 
- CASE.  ..  fiUKBER - U - 


-  - -flOtBCULAfi  -  STSUCtU**—. - - 

- - - _H - R— - - - - 

•  • 

C  i  C 

-  _ _ - - 

H-C  C  C*M 

•  ■  • 

- Hmt~  — C -C-H- - 

i  -  •  * 

C  C 
H  H 


ATOMIC  CDHPOSJTION  OF  MOLECULE 

A 7  OH  MJHBER . 

H  B 

- C - 10 - 


OR  CUP  STRUCTURE  AND  COH POSITION 


uRCilF 

NUMBER 

1 

— UKQLE— 

height 

4016 

- CHEMICAL  -SYKOI 

OF  COMPONENT 

c 

.  . . c  __ 

It - QHXD-BflH — 

COORDINATE 

3 

_  7 

COLUMN 

COORDINATE 

3 

3  —  . 

_ BC^D _ 

VECTOR 

5 

ROND - OfiOLP— NUKEB* 

TYPE  OP  COPE  A  TOP 

1 

c 

H 

3 

3 

! 

1 

2 

7 

2 

1 

2 

0 

2 

ion 

C 

3 

3 

2 

C 

c 

3 

3 

3 

7 _ 

6 

.4 

2 

-1  - 

1  .  - 

b 

1 

5 

1 

,  I 

0  •  „ 

..  a 

mu 

c  - 

t 

1 

c 

9 

3 

4 

2 

3 

c 

. h- - 

3 

3 

. . -1 - 

1 

- 7  — 

1 

- 1 - 

1  _ 

4 

999 

c 

3 

7 

4 

c 

3 

9 

2 

1 

6 

c 

3 

e 

1 

2 

3 

1046 

c 

9 

3 

3 

c 

. .  Q  _ _ 

7 

7 

7 

3 

2 

o 

1 

7 

p 

11 

3 

3 

1 

0 

.  .  ^ . . . 6  . 

„  .1016 

.  c 

. 3  . 

. 9  .  . 

_ 

«  _  -  .  .  _  - . 

c 

3 

7 

6  ■ 

V  ' 

4 

c 

3 

U 

4 

2 

B 

7 

h 

1 

9 

1 

1 

0 

7 

999 

C 

_ _ c . . 

7 

. . . 9 . . 

7 

. . 9 

4 

1 

7 

g  .  . 

c 

3 

7 

1 

2 

4 

c 

9 

3 

6 

1 

3 

a 

~iou~ 

9 

11 

B 

c 

3 

9 

8 

2 

6 

_ — ...  .  _ _ _ 

. c 

...  .  7  ... . 

- U 

—  3  — 

-  -1  - 

- ID-—*.  -  - - -  -  — 

9 

ni6- 

H 

_ _  _  C _ 

9 

. . s _ 

13 

_ g _ 

3 

1 

0. 

C 

7 

11 

2 

2 

10 

c 

7 

7 

fl 

1 

7 

* - - - - - . 

.  - - 

. . -N 

11 

-  .  9 

. 3  - 

•  •  -4  — — 

-  0 .  . . . . . 

40 

1016 

c 

7 

11 

to  ’ 

. 

_ _ 

G  - 

_ 8 _ _ 

- 44 - 

- £ - 

- i - 

- e - : - 

C  9  9  6  2  9 

H  7  13  3  1  0 


RING  NUMBER  1  IS  6  memFEREU 

c-  - - -  - - JSROUP-mjMBERS-OrJUNG  . CONSTITUENTS - 7* - 4<  —  6* - D,_  -IQ-i- —  -9-«  - 

RING  NUMBER  2  IS  6  Kf^eFREU 

. . . GROUP.  NUMBERS  Of  .RING  CONST ITiifcNTS.  .  .  __1* _ 3,.  .S»- . 7, _ 4,...  _4’, _ 
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ATCM1C  CCKpeSlTlOK  OP  CORE  ATOM  ANO  LIGANDS 

. GROIF  NUMBER  - 

1 


EACH  GROUP 

n  ; . 

o 

0 

— 0- - 


.  ..  4  ... .. 

- -5 - 


LIGANP  .  . 
L1UANP 

CORE 

_  LIGAND _ 

7 

1 

0 

__5 _ 

9 

0 

1 

7 

0 

0 

0 

_ 9 _ 

0 

0 

0 

_ 

LIGAND 

.2  . 

2 

0 

a 

LIGANp 

1 

0 

0 

0 

CORE 

•0  .. 

1  ... 

0  - 

0 

LIGAND 

4 

5 

0 

0 

ligand 

4 

4 

0 

0 

_  CORE  _ 

_ 4 _ 

„ _ ^ _ 

_ $ _ 

LIGANP 

4 

6 

0 

0 

LIGAND 

3 

3 

0 

0 

--  LIGAND  - 

.  1... 

-  0 

--0  . 

-o  - 

CORE 

0 

1 

0 

0 

ligand 

4 

6 

0 

0 

. . LIGAND--- 


LIGAND 

1 

0 

0 

0 

7 

CORE 

0 

1 

0 

0 

..  ... 

.  LIGANP 

■  —  4  - 

...4. 

-0- 

— 0 . 

LIGAND 

4  . 

5 

0 

0 

e 

COKE 

C 

1 

0 

0 

LIGANP 

5 

7 

0 

0 ' 

LIGAND 

2 

2 

0 

0 

LIGAND 

1 

0 

0 

0 

9  . 

CORE 

0  - 

1 

-  -  c 

. .  .0 . 

ligand 

7 

9 

0 

0 

LIGAND 

1 

D 

0 

0 

10— 

_ _ CORE  — 

_ a _ 

.  4„- 

_ 0—  -- 

_ o„... .. 

ligand 

6 

6 

0 

0 

LIGAND 

1 

1 

0 

0 

.  -  ligand 

1 

•  0 

-  a 

0 

NUHPER  OF  CORE  ATOMS  IN  LONGEST  CHAIN 
GRCUP  PLM&ER  OF  CORE  CONSTITUENTS  OF  tONOEST  CHAIN. 


11 

2#  - 
2. 


7. 


10. 


EXTERNAL  rotational  symmetry  CONTRIBUTION 

SCONCE  SYMMETRY  numher  entpofy 

INPUT  4  "2.75453 

- - NOT  -CALCULABLE.  - - -  -  1 -  - - - Ot-OCOCO - 


CPTUAL  ISQHER  CONTRIBUTION 

.  .  . _ itQL'kCE _ A3YMM.tTF.tC. _ PSELCC  ASYMMLTHIC ...  ENANUOMEBS  .  _hfcSC _ TCUL  DPT  1CAU.  ISOMERS - EN.tRCPY. ...  - - 

Cf  I'PIjTEE  C  fj  0  0  0  0,00000 


TLERMCCHfcMICAL  PROPERTIES  CF  MCLECULE 
hem’  capacity  COEFFICIENTS 

CAL/DEG-MOUE  CAL/tER4**-MULE  CAL/D£G**3*MQLE  CAl./CEf!**4-MCLE 

-1 ,360*'5232E#01  1.S7339V39E-01  ■  l,2a79494:E«.i;4  3 . 3355Q42OE-O0 


C((?)  - 

CAl/WEG-MOl  F. 
31,6 69907 
3i,9ofl50i;  - 

42,«609b4 
!*2, 036091 
itf  ,639P5«;- 
65 ,66669? 

TC ,91844? 
74,995344 
7  fl , 297516 
*3,370199 
-87, 7615 AU— 
93,040579 
100.840390 
112,737921  - 


CAL/13Pfi-MOLf 
00,270073 
.  «0. 432710 
91.157058 
101*73fil£3 

-111  .9r*)l«2 - 

121 ,59770fc 
130,77471ft 

139,331675 
147.410053 
162,157277 
-±75c337982— 
1A7. 363129 
196,760141 
-  209.957459  — 


— »<<}<T)-L<29e  >  MI - MAT  Mw(298  ) — 


cal  /  rf  f-mcle 

8  Q  ,  2  2  0 1'  7  3 
80,22.0761 

fll,612«54 

B4.5M991 

80 , R98873 - 

92, 370972 
96 , 6  D 1C  1 2 
1OO.076H3 
105,129844 

113,429990 

12lr348AC0 - 

128  ,85  9625 
135,996245 
142,633724 - 


KCAL/M0L6 
0,000000 
,063579- 
3,017662 
8,576596 
1 4-,  172  791— 
20 ,4^0749 
27.3C6964 
34,609043 
42,280200 
58,472744 
- — 7  5^56  49  9ft— 
93,637607 
112,971413 
134,247469 


-.^)EWHf429d)*L(T)-H(39iF-- 

HCAL/MCLF 

36,000000 

36,143570 

39,897602 

44,656596 

. . 50,252791- - 

56,538749 

63,306964 

70,609943 

70,360200 

94,552744 

- Jrl  li  6  6  4  9  9  0 - 

129,71760? 

149,051413 

--  - 1 7 0  r33 7 469 
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GROUP  Actftmtrv  TMER*OCHE*ICAt  PROPERTIES  PROGRAM 


H  H 


-  •  ATO*U  — COMPOS  l  T-I  OK — OF  — HflfcSCtiWS - 

ATOM  NUMBER 

M-  16  -  • .  -  - -  - . - 

c  b 


GROUP  STRUCTURE  AND  COMPOSITION 

OROUP  QRQtP  “  CHEMICAL  STHROU  URID  ROW  GRID  C0LtKN~  BCND  BONO  GROUP  NUMBER 
NUMBER  height  of  COMPONENT  COORDINATE  COORDINATE  VECTOR  TYPB  OF  CORE  ATOM 

—4 - - -C - —  9 - 3 - 1 - 

C  7  3  4  1  2 

H  3  1  8  10 

h  .  3  9  2  l  ■  -  —  ■  0  - -  - . . 

h  7  1  6  1  0 

- 2.  .......  502 - C 

C 
C 
C 
H 

C  11  9  5  1  6 

C  b  U  2  1  5 

M  ~  9  7  6  1  0 


— RING  NUMBER  1  IS  5  MfeMBERtO  -  „  _  -  > 

GROUP  NUMBERS  OF  RING  CONSTITUENTS  •  2,  3,  it  7,  4i 


-  ATCM1C.  -  CCPPCS IT IOK  OF. ...CORE  ATOM  .  ANQ._.UGANBS- IK- •  EACH. ..GROUP  - 


CROUP  HUMBER 

- 4 - 


SU3GHOUP 
-CORE- - 


LIGAND 

LIGAND 

. LIGAND 

LIGAND 

COKt 

- UQAMJ- 

LIGAND 

LIGAND 


13 

1 

.. 

1 

0 


iz 

l 


c 

— 4— 

7 

0 

...  0  ■  - 

0 

1 

- 1— 

6 

0 


3  - - 

...  CORE  - 

LtGAND 

LIGAND 

_ ligand— - _ 

0  .....  1  - 0 

4  2  0 

5  2  0 

0 

0 

4 

LIGAND 

CORE 

LIGAND  - 

1  0  O' 

0  1  0 

-14  -7  -  -  0 

0 

0  , 

•  -  0  - . .  -  . 

•c. 

LIGAND 

LIGAND 

- COHfc-  - 

LIGAM1 

LIGAND 

-  ligand- 
ligand 

CORE 

- UGAKD- 

LIGAND 

LluAND 

-  -  LIGAND 

CORE 

LIGAND 

LIGAND 

LIGAND 

LIGAND 

. .  -CORE- 

LIGAND 

LIGAND 

- L I OAND 

LIGAND 


X 

\ . 

11 

3 

—  1 . 

1 

0 

10 - 

4 

1 

...  l  ...  - 

0 

12 

2 

1 

1 

...  0  --- 

11 

1 

- 1 - 


NUKPE"  OF  CORF  aTOHS  in  LONGEST  CHAIN 
CROUD  M.  i  bYh  DF  ~CQRE  C 6 K ST  I TL' E N f$  CF  L 0 NOEST  CH AIN ” 


It 


NOhPRcur  interaction  and  internal 


TYPE 

T»inr,  NHM.HL^  C  R 
GROUP  NUl'.l'p.E 

AMCLNT 

HFaT  of  FORMATION 

kcal 

RING  CCHHbCT  ir.t.-  - 

1  - 

1 

6,3u0 

CIS 

2  3 

1 

1,000 

GAUCHE  ALKAf'L 

3  5 

1 

,aoo 

J  NT  LKtiAL  KOT-TlilN - 

. . 1  -  ■  — 

. 1 

. 

INTERNAL  ROTATION 

6 

1 

RCTATICNAL  SYMMETRY  CONTRIBUTIONS 

ENTROPY  HE«f  CAPACITY  CCEFfJClFNTS 

CAL/Ott)  K  CAL/D&G  K  CAL/DEG  K**2  CAL/CE&  ***3  CAL/CEG  K««4 
2  7 » 300  -l,2010t*Ql  1,73226-02  •e.&AJlE-Ofc  1,6*566*09 

0.000  -2,719S£*00  9  j 8531 E -03  -4.7848E-06  1,3*426*09 

-2,183 - - - - - - 

-2.183 


CCNTR  IBUT  I  ON 
EUhCFY 

0 , oooto 
0,00000 


EXTERNAL  ROTATIONAL  SYMMETRY 

SGLPCE  SYMMETRY  number 

ILFLT  1 

KQT  CALCULABLE  l 


CPTICAL  ISOMER  CCNTHIBUirOK 

S2U-XE  ASYMMETRIC  PSELCCASYMMETRIC  ENANTIOMERS  M£SC  '  TOTAL  CpTlCAL  ISOMERS  EMNCPY 

COmH'TEL  0  (j  0  0  0  0,0000c 


TMERMCCHEHKAL  PROPERTIES  CF  MOLECULE 
VIEaT  CAPACITY  CCEFFlClfm 

CAL/DEG-MCLF  CAL/rEG**2-MULfi  CAL/DkG** J-MCLF  CAL/CEG*  •  4-pCLE 

•1,91,3403776*01  2,1*4^739^6 *<<1  —  ~«l-,a0b3221'-h-G4 . -3-, 05* 1-6293E •■06-  - . 


T 

r.  tr ) 

c 

•  <G(TMi(?98)  )/T 

MT>-M?9fU 

DEIHF  <  290 ) *M ( T ) -M 

pt  .r.  k 

CAL/i:i-r.-r  oi  L 

CAL/CFG-MOLF 

KCAL/MOLE 

CAL/ff-G-MCLE 

kcal/mcle 

kcal/mcle 

.796,0 

34,213673 

9  ? , 2$  35  4  u 

2  t  7  04093 

97,2915*8 

o, cotton 

-35,660000 

1CL.U  -■ 

-  .  .U,bi'tD55.  - 

07,5?J375 - 

2t772fll3- 

. -•••  97,29  4306— - 

.  r0*672i — 

— - -35,591279 

*OC,0 

4  c,. 142360 

109,365060 

6,921736 

90,62095? 

4,217643 

-31,442357 

t  C  G  ,  o 

59,694259 

121.403*57. 

12|33(J5.-1 

102,134640 

9 , 634  4  Oe 

-2  6 ,02  558  2 

t  C.  G ,  y 

69, <545741 

1 33, 236369  - 

l8t6A39:if, 

106,336677-  .  . 

16.139015- 

-19,520185 

7  0  C  ',  (, 

7  f1  ,  4  (i  Q  0  6  2 

144,677750  ‘ 

26(2770/8 

111,002057 

23,572906 

-12,007014 

MJC.p 
-  -UOC.y. 

l'-5|  6F0524  . 

91.730321 

155,64117? 

1  A  A  .  09  2S  Af. 

34|49546l 
„  43j374a2A 

115,901961 

3 1 , 791360 
_ *0,670734 _ 

•3,868632 

f  01A736- 

....  -  ^  A  f  4  il  “4  4 

- 4,  C*  O  t  y 

U00,0 

^6, 612721 

176,023051 

52(609276 

125,922665 

50,105186 

14,4*5186 

U’CC,  J 

10 A , Sc 6331 

194,422393 

73(011443 

135 , 632934 

70,307351 

34,647351 

1*00,0 

111,133353 

211,065734 

94,62123* 

1*5,410*30 . 

81,9}7145 

56,25714? 

1 6  C  0  i  U 

117,25.37  61* 

226,299229 

117 1*51145 

154,542321 

11*. 7*7052  ■ 

79,007052 

If’CC,  0 

124,635625 

240,516393  . 

141(606046 

163,346197. 

.  130,902753 

103,242753 

24CC  ,0 

134,74487? 

254,140926 

167(48722? 

171,749361 

1**, 703130 

129,123130 
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CROUP  ACCmvm  THfcRHOCHEM ICAL .  PROPERTIES  pACORAH 
- CA8£..-JIL*BER _ 11 _ 


-MOLECULAR  STRUCTURE 


H  M-C-M  H 


- ATOMIC— MWWUOK.  -Ot— M0L6CLLE- 


ORCUP  STRUCTURE  and  composition 


CHEMICAL  STpnfIL 
Of  COMPONENT 


OR  ID  ROW  GRID  COLUMN  BCNC 
COORDINATE  COORCJNATE  VECTOR 

- 6 _ _ _ d - - - - 

6  6  3 

3  18 


CROUP  NUMBER 
OF  CORE  atoh 


— COMPLEX  LIGANDS 


of  CORS  ATOM  —a  ARt-COUAtt - 


COMPLEX  LIGANDS  1  AND  4  OF  CORE  ATOM  2  ARE  FQUAL , 


COMPOSITION  OF  CORE  ATOM  AND  LlCANDS  IN  EACH  GROUP 


QROUP  NUMBER  SUBGROUP 

1  CORE 
-  .  ligand 

ligand 
ligand 

- - LIGAND  — 

2  CORE 
LIGAND 

-  _  .  ligand  — 

LIGAND 
ligand 

-1 - CO«6 - 

ligand 
ligand 

. -  -  LIGAND--  - 

LIGAND 

4  CORE 

- LIGAND - 

ligand 
LIQANO 
LIQAND 


o  0 

0  o 

-0 - 0— 

1  0 

1  c 

-l — 0- 

1  0 

0  c 

-1 - 0-~ 

3  o 

0  0 

c  —  0  -  . 
0  o 

1  C 


NUMPER  OF  CORF  ATOMS  IK  LONGEST  CHAIN  • 
..CROUP  ..NUMBER..OF  CORE  CONSTITUENTS  Of  LONGEST  CHAIN  .■ 
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NONURCUP  INTERACTION  AND 


INTERNAL  RCTATtOKAL 


HE  TRy 


CONTRIBUTIONS 


TYPE  _ .  RING. NUMBER  OR.  ...AMOUNT.  HEAT-OF  FORMAT  ION  ...  .ENTROPY  - 


GROUP  NUMBER  KCAL  CAU/DfcO  K 

JNUHNAU  ROTATION  1  1  -2,183 

IKTfcNAiAl-NOTA-UnA, - S -  1 -  -2,183— 

INTERNAL  ROTATION  A  1  -2,183  . 


- . -HEAT -CAPACITY.  COEFFICIENTS . 

CAL/CEG  K  CAL/CEG  M*2  CAL/CEG  K*«5  CAL/CEG  «**4 


EXTERNAL  ROTATIONAL  SYMMETRY  CCNTPIBUTION 

SOURCE  SYMMETRY  NUMBER  ENTROPY 

COMPUTE!  3  *2,10323 


SOURCE 

COMPUTED 


! 

! 


. . . . OPT  I  CAL — ISOMER  — CCNTH1BUTI  GN - - - 

ASYMMETRIC  PSELT,CASYHHrTHlC  ENANTJGMEHS  MLSC  TOTAL  OPTICAL  ISOMBRS  ENTRcPY 

■  -  G  -  . 0 . 0  -  0  - - - ■■■<  ■  -0  - -  0,00000 


THEPHCCMEMICAL  PROPERTIES  CF  MOLECULE 


cal/deomole 

-2,2270fl940E*00 


MEAT  CAPACITY  .COEFFICIENTS  . 

CAL/CEG**2*H0LE  CAl/DkG*  •3-MOlR  C  AL /CEG  •  •  4 -POLE 

0,96394 1& IE* ^2  *5,943346B2E-C5  1 , 21 7 N5257E-0B 


■aii 


T 

G(r  j 

§■ 

H(T)-hC  ) 

* ( G ( T ) -H ( 296 )  >/T 

PIT  )-M(290 ) 

DELHF(290)+H(T)-H(29fl) 

CLU  K  .. 

.  „_CAL/CEG*MOLE - 

-CAL/ LEG- MOLE  - 

kcal/mglf. 

-  -  CAL/CfcC-MCLE 

-  KCAL/hCLE- 

...  -  -  XCAL/MCLE  - 

29 n ,  o 

22,9?4630 

7C,4?7u  HZ 

3,304418 

70,42708? 

0,000000 

-32,140000 

JOO,  u 

?3, 0*454 3 

70,5P5,P9? 

3 | 39 0407 

70,427594 

i  045909 

*32,094011 

ACD.tl-  - 

- 29,618745 - 

...  7B,i34l?y  - 

5,9916  79-... 

-  71,416017-.  — 

2,6872*1- 

. . . -59,452739  - - 

soo ,  n 

35,356567 

85,374275 

91246V44 

73,489223 

5,942526 

*26,197474 

aco.o 

40,351073 

92,272745 

13 i 030216 

76,049749 

9,733798 

•22 • 406202 

7C0.0 

-14,075337  ..  - 

96,025761 

17,294817 

.78,639477  —. 

13,990399 

-18,149601  . 

F'CC.O 

*B,  4:?242R 

105,04i!fi63 

21 1 953377 

01,729664 

18,640959 

•13,491041 

Per ,  o 

r-1 , 605421 

110,931560 

26,957833 

84,649988 

23 ,653415 

•6 ,4P6565 

■  lDCO.D  . 

-  .  .  54,357384 . 

116,514802 

. 32,259427 

■  ■  87,559793. 

-  re, 9 55 13 tio 

■3,184991 

12CU.U 

5c,  B  ;i  05  Ob 

126.836594 

43,595542 

93,260661 

40,291124 

8 , 151124 

1400,0 

62,316370 

136,172933 

55,717813 

98,734793 

52 , 413395 

20,273395 

IfctiC  ,0 

65 , 4G9539 

144,701845 

68,499244 

.103 , 955078 

65,194326 

33,054626 

1HOL.O 

*8,904505 

132,607342 

81,929754 

108,926600 

70 , 62533* 

46,485336 

20C0.O 

73,146076 

160,076075 

96,116179 

113,670957 

02,811757 

60.671757 
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OROUF  ADDITIVITY  THERMOCHEMUK  PROPERTIES  PROGRAM 

CA»E  NUMBER  14 

MOLECULAR  STRUCTURE 

N 
C 

C 


»  • 

NC»C  C*CN 

■  •* 

NC*C  C-CN 


f 

N 


ATOMIC  COMPOSITION  OP  MOLECULE 


ATOM 

- «L 


number 

u 

— e 


—GROUP  - 
height 
1399 

GROUP  STRUCTURE  AND  COMPOSITION 

-CNBM4CAL  STWMJfc — -ARID  RON  —GRIP  COLUMN  BOND  - 
OP  COMPONENT  COORDINATE  COORDINATE  VECTOR 

C  A  4 

BOND 

TTPB 

NUMBSR 

1 

-ONOWr  NIIMV0N . 

or  CORE  ATOM 

1 

L 

c 

CN 

4 

6 

- n 

6 

2 

6 

7 

2 

7 

6 

- 4 - 

1 

4 

1 

t 

1 

1 

1 

3' 

2 

0 

2 

1399 

c 

c 

4 

« 

' 

CM  w4 

1 

c 

CN 

c  - 

c 

c 

Z 

10 

6 

A 

6 

4 

2  - 
1 

1 

2 

0 

. 1399  * 

9 

1 

...  -  on  ’■ 

8 - 

* 

- f  — 

O  '  "  "  " 

4 

1399 

c 

6 

0 

'  c 
c 

0 

0 

0 

9 

~Z 

1 

\  * 

CN 

6 

10 

3 

1 

9 

1399 

c 

10 

6 

1 

c 

0 

e 

2 

2 

* 

CN 

12 

.....  # 

9 

1 

1 

0 

Y 

. A 

1399 

c 

* . 

c 

6 

0 

1 

1 

— 

c 

CN- 

10 

6 

—to  ■ 

6 

- 3 - 

2 

“1 - 

9 

RING  NUMBER  I  JS  A  MIMBEN6U 

GROUP  NUMBERS  OP  RING  CONSTITUENTS  ■  l\  J,  9,  A,  4,  2, 


ATOMIC  COMPOSITION  Of  CORE  ATOM  AND  LIGANDS  ]N  EaCH  QROUP 


croup  number 
i 


s 


•4 - 


9 


A 


SUBGROUP 

CORE 

- CIGAND— 

LI04ND 

COPE 

-  LIGAND 
LIGAND 

CU«fc 

— LIGAND — 
LIGAND 
LIGAND 

—  CORE 
LIGAND 

ligand 
— LfOAND- 
CORfe 
ligand 
LtGAhD 
LIGAND 
CORE 

— ligand- 
ligand 

LIGAND 


C 

1 

“10 

1 

1 

10 

1 

1 

-6  - 

2 

1 

1 

2 

6 

-1- 

1 

6 

4 

1 

1 


-O’ 

0 

0 

o  -  • 

0 

0 

"0 - 

0 

0 

o 

a 

0 


- NUMBER- or  CDRg  ATOMS  IN  LONQGST  'CMAf~N — ■ - 0 - 

GROUP  NUMBER  Of  CQ&E  CONSTITUENTS  tf  LONGEST  CHAIN  •  1,  J,  9.  A,  4,  g« 
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TYPE  MKQ  NUMBER  CR  AMOUNT  HE*T  Of  FORMATION  EhTRO*Y 


- 0*0UM -NUMBER-  - .  - . ~KCAL - :CAl/D£0  X 

ORTHO  CORRECTION  1  2  1  ,5*0  -1,610 

ORTHQ  CORRECTION  13  1  ,370  -1,610 

- 0*T»0  CORRECTION - 2 - 4 - J - ySTO - -  1,6 10 — 

ORTHO  CORRECTION  35  1  ,370  *1,610 

ORTHO  CORRECTION  A  6  1  ,970  -1,610 

— gRTHO-ccRRemcN - 9 — 6 - r . — i»ro - »i,6io  ■ 


HEAT  CARACITY  COEFFICIENTS 
CAL/DE6  K  CAL/3EC  K**2  CAL/DSG  K-«J  CAWC83  K*T 
3,67656.01  *»40A7B.03  -6,26666.06  2,0P8CEo09 

3,67696^01  4,40476-03  *6,26666-0*  2,09*06-09 

3,6765E.Cr — 41!  4  0476-03  -6  >26666*06 - 9  ,"06606909 - 

3 , 67696*01  4,40476.03  -6,26666-06  2,09806*09 

3 , 6769E.01  4^40478-03  -6,264*6-06  2,09*06*0* 

3,47456*01  Hiimim  1,0^9B»0f  " 


external  ROTATIONAL  symmetry  CONTRIBUTION 

source  . -  SYMMETRY-NUM8BR  ENTPORY 

CCHPUTEC  IS  *4,93816 


OPTICAL  ISOMER  CONTRIBUTION 


SOURCE  ASYMMETRIC  PSEUD0ASYMMETR1 C  ENANTIOMERS  MESC  TOTAL  OPTICAL  150MIRI  ENTROPY 

C0MPUT6C  0  00  0  0  0,00006 


,  THERMCCHGMKAL  pROPEBTiea  QF  MOteCULt  - 
HbAT  CAPACITY  COernctPNTS 

-tAt/&ttwftCLE  CALVCfft»*2-M,JLE  - C^L'/OEO.  .3-MOLE - CAL/OEO*  >4-M0lC -  " 

2,099246096*01  2 , 03406124E-Ol  •! ,  9961219H.04  6 , 82401839E*00 


T 

PEG  K 

C<p) 

CAL/DfcG.MOLE 

S 

cal/deg-mole 

H  ( ▼ > -M  C  0  > 
KC*L/h0lE 

•(O(T).Hf 296 > )/T 
CAL/DEG-MCLE 

NCT  l-H 1 290 1 
KCAL/MOLE 

08LMF<298i*H(TJ.H12t«> 

KCAL/MOIE 

- - 

- 69  ,'6 P 7 4  J3 - 

■~1 0  8,40104  4 - 

1 3 ^  6  6l233 

- 108,401844 

“  0,000000 — 

'  — . 216,220000  • 

320,0 

65,892093 

103,841911 

13i792813 

108,403312 

'  ,131580 

216,33x380 

400  ,0 

74,764739 

129,07699a 

2Di847961 

111,110128 

7,186748 

225,406748 

— 

- 5J0tG - 

—41,322904  - 

146,507039  - 

^  26 j67l27fl 

'  116,406943 

13,010043 

233,230043  -  “  -  ' 

fcOO  ,  0 

7  00 , 0 

65,916032 

161,746361 

37  A  04  7727 

122 ,7nBB9l 

23,386*94 

241,606494 

4d,97i5ol 

175,257  ICO 

45i8C329B 

129,339663 

32,142064 

250,362064 

■  *00,0 

9CT*04935~ 

'1877274202 - 

54} 804 9 01" 

- 133,844617 - 

'*1,1*3660  - 

239,363668 

900,0 

92,119994 

198,056421 

63 1 96  rj394 

142,146*64 

*0,29*161 

260,319161 

1000 , 0 

- i2oo,  cr  — 

93,026978 

99,557689 

207,610057 

224,963669 

73 j  218377 
92A042937 

1*8,232713 
.  159,6*5382 

59,5573*4 

78,341704 

277,777344 

296,601704 

1400,0 

101,772597 

240,092462 

UlA6a7266 

170,0738*7 

98,026033 

316,246033 

1600,0 

114,947231 

254,446011 

133121595b 

179,725108 

119,594724 
144,687276 — 

337,774724 

“  '  t-frpo  ,0 

— j-36l  35712  Q— 

4PT, 470*1* 

188,83899? 

362,607276  _ 

2000,0 

175,277793 

289,976127 

1B9a459925 

197,676981 

175,791291 

394,018291 

50 


- ATOM  ft— COMPOS  I HCN  OF  HOLBCUkE  - 

ATOM  NUMBER 


H  -  4 


51 


KING  NUMBER  i  [S  4  MEHBeRtU 
CROUP “NUWBERE  Of  ■RrNG  'TC^rsi  tTUEHTS 


r,  5t  7#  6.  3, 


TUNG  NUMBER  2  IS  6  HEKbFRED 
OROUP-  NUMBeRS-Cr  RING  CCNSUTUtRTS- 


*«  10. '  X7t  -  H»  13. 


COMPLEX  L 1 Q AND S  7  AMO  9  CP  CORE  ATCH  8  ARE  EQUAL, 


52 


ATOMIC  COMPOSITION  or  CORE  ATOP  AND  LIGANDS  IN  SAC*  GROUP 


GROUP  NUMBER 
1 


SUBGROUP 

CORE 

— LIGAND — 
LIGAND 
LIGAND 
—  CORE  “ 
LIGAND 
LIGAND 
— trtGAftO — 
CORE 
LIGAND 
“LIGAND  - 

ligand 

LIGAND 

- CORE - 

LIGAND 
LIGAND 
CORE ' — 

.ligand 

LIGAND 


etganu - ff- 


L I Q  AND - 16 - 13 - 

LIGAND  X  0 

LIGAND  2  1 

LIGAND  -  ■  1-  0 

CORE  0  1 

LIGAND  14  12 

LIGAND - 4 - 8— 

LIGAND  1  0 

LIGAND  1  0 

-CORE'  ~  0  0 

LIGAND  19  l» 

LIGAND  1  0 


- : - trtmafcR-ap  CORr^TDMS  “1* -LONGEST  C**IN — ■ - 13 

group  number  of  core  constituents  or  longest  chain  »  4* 


« I  NO  CORRECTION 
-^•INC-CCRWT-Km- 
QAUCHE  ALKANE 
QAUCHE  ALKANE 


nongrcup  interaction  and  internal  rotational  vmHWf  contribution! 


INTERNAL  ROTATION 
INTERNAL  ROTATION 


GROUP  NUMBER 
1 

— z~ 

7  9 

9  e 


1 

l- 

4 

A 


KCAL 

•  ooo 
-  ,000- 
3,200 
3,200 


-  NEAT  CAPAC ITT  COEFP ttltNT0 - - 

CAL/OEG  K  CAL/DEQ  K  CAL/OEG  K**2  CAl/DEO  K**3  CAL/W  K**4 
10,800  ■l,3im*01  2,55478-02  -1,1337E«09  l,033;E*o? 

18,000  ■-  -  -1,3120E*01  8f5547i«0a  -1,13376*0»  -1,0**16.09 


— iT*7f' 
-1.377 


10 


EXTERNAL  ROTATIONAL  SYMMETRY  CONTRIBUTION 


SOURCE  SYMMETRY  NUMBER  ENTROPY  ' 

COMPUTE!  2  *1,37746  ' 


- : - cmeRt  — tSOHfeft- — CONTR10UT  TON - - -  - . . 

SOURCE  ASYMMETRIC  RSEUDOASYKHETRIC  ENANTIOMERS  PESO  TOTAL'  OPTICAL  T30HERI 

-COMPUTES - -  - . j . - . . —  *  C —  ‘0  0 


ENTROPY 

0,00000 


thermccmemical  properties  or  molecule 
-•  -  •  — heat  capacity -eoerneiENTS 

CAL/DEG.MOLE  CAL/CEG**2-MULE  CAL/DEO--3-MOLE  C*L/nEG*-4-K0LE 


▼ 

•3,103317556*01  3,J36l7i0$E**0l  »2 

,14954274E*fl4 

5 , 537166846*00 

T 

C(P> 

S 

•  ID<T)-H<29aM/T 

Hm-w<298) 

DSLHf I299I*NIT)pH(29#I 

- 

CAt-/Dfcf!»HOL£  - 

"CAL^DEG-MOtP' 

-mcal/mole 

C  Al  /DTO-MC 

kcal/molp 

kcal/kolf 

298,0 

46,990056 

104,107616 

3*0*9*54 

184,107616 

0,000000 

•148,480000 

300,o 

A7,41HU7 

104 ,503342 

3|184350 

184,106668 

^  1  0944Q2 

•140,385990 

- -  ■  4  o  ore  “ 

- 66  't7°T77o 

2tf0 ,872353 

~ti*92r,T06 

100^29527  3" 

-  5,833532 

-  *142,029180  - 

500,0 

83,21-0263 

217,592487 

16*444305 

190,863765 

13,394331 

•135,103649 

ftQO.O 

96,9H383u 

234.Q16329 

25*474906 

196,710739 

22,304554 

*•126,075446 

~  ’00 iC 

- 100,483687 

249,897266 

35*764214 

-  203,179752 

32,674260 

■115,785740 

000,0 

117,980068 

264,979272 

47*099463 

209,987366 

44,009509 

■104,450491 

’00,  c 

125.&C5205 

279,33979c 

59*299919 

2l6,8M7i? 

90,209565 

*92,250433 

1 C  C  0  •  0 

1"3  2,35x325 

292,94  39^4 - 

72* 210870 

223, 816608 

6  9 , 1 2  6  9 1 6 

•■79, 333004  . .  ’ 

1200,0 

142,645440 

310.O4U073 

99*779517 

237,466237 

90,689563 

•91.770437 

nao.o 

152,090254 

340,793900 

129*271760 

250,024047 

126,181805 

*22,270195 

lftOO',0- 

- 162,743607 . 

-  361,727007  * . 

160*709921 

263,215578 

157,019567 

9,15950? 

1000, c 

177,463330 

301,693217 

194*040294 

275,276301 

191,550340 

43,090340 

2000,0 

198,907208 

401,430045 

232*143138 

200,903453 

229,033104 

00,593104 

GROUP  STRUCTURE  I  NO  C0HPQ8  fT  1 0 K 


GROUP  QROLF  CHEMICAL  SYMBOL  GRID  ROW  GRID  COLUMN  BOND  BOND 
number  weight  of  COMPONENT  COORDINATE  coordinate  vector  type 
- 1  -  bl5  C .  7  4 


C  10  451 

N  4  7  2  1 


GROUP  NUMBER 
Of  CORE  ATOM 
1 
3 
2 


RING  NUMBER  1  IS  6  MEMbEHU* 

RING  NUMfcfR  2  IS  6  MEMfaFRfcO 


tm tf  mim*- 


> . 

6 


CURE  ATOP 

and 

LIOANDS 

— SUB GROUP - 

M 

*  *c~ 

0 

OGRE 

0 

1 

0 

ligang 

10 

5 

c 

—  LiQANr  -■ 

-  1- 

0“  ■ 

c 

ligand 

1 

0 

0 

CORE 

0 

'  0 

Q 

L"1  G  AND - 

0 

4 

"O'" 

lIGAND 

4 

2 

0 

CQ«E 

0 

1 

: 

—  LIGAND  - 

0 

4 

0 

LIGAND 

2 

1 

0 

ligand 

1 

0 

0 

- LtGAND 

'1 

0 

"0 

CORE 

0 

1 

0 

LIGAND 

0 

4 

0 

LIGAND 

2 

“1 

0 

LIGAND 

1 

0 

0 

ligand 

1 

0 

0 

CORE 

0 

1 

0 

LIGAND 

0 

0 

0 

CORE 

0 

0 

D 

—  L1GANO--  ' 

-4 

2 

0 

LIGAND 

4 

2 

0 

CORE 

0 

1 

0 

- tlGAND 

to~ 

*0~ 

ligand 

1 

D 

c 

LIGAND 

1 

0 

c 

-  —  CORE  ” 

0 

-  %' 

0 

LIGAND 

0 

0 

0 

C 

2 

C 

0 

1 

V 

0 

0 

2 

0 

0 

“0“ 

0 

2 

0 

0 

0 

0 

0 

1 

0 

1 

0 

2- 

0 

0 

■  0 


NMPER  OF  COPE  ATOMS  IN  LONGEST  CHAIN  ■  a 


CROUP  MP86R  CF  C0«E  CONSTITUENTS  CF  LONGEST  CHAIN 


1 1  3 <  6i  6*  5 ,  2 «  4 1  7» 


nokgRclp  interaction  am  internal  rotational  symmetry  contributions 


TYPE  ring  number  CR  AHOUNT  HEaT  Of  FORMATION  ENTROPY 

GROUP  NUMBER  KCAL  C Al /DEG  H 

RING  CORRteTICN - j — 2 — a - 1 - 3rw - — '•  — 


HEAT  CAPACITY  COEFFICIENTS 
CAl/DEG  K  CAL/06G  K*«2  CAL/DEG  ***Z  CAL/CEG 


EXTERNAL  ROTATIONAL  SYMMETRY  CONTRIBUTION 

- SUli»rCk - SYMMfcTRV-NUHBPR - - entropy 

NOT  CALCULABLE  1  0,00000 


OPTICAL  ISOMER  CONTRIBUTION 


SOURCE  ASYHH6TF1C  FSEL'DOASYMKETR  1C  ENANTIOMERS  YtSO  TOTAL  OPTICAL  ISOMERS  ENTROPY 

COMPUTED  0  00  0  '  C  •  0,00000 


-RlYtROfY'-AhO-PPeE-tMtMtiY-  U*rA-Bf*ttM-MAY  BP  *fTP  A  FEM  PERCEnT  DUE  TO  POSSIBLE 
CONTH I EUTIOM 6 )  FROM  EXTERNAL  "OTATJONAL  SYMMETRY  OR/aND  OPTICAL  1  S OHE R  1  8H  t 


THERmCCHLMICAL  PROPERTIES  cr  MOLECULE 


heat  CAPACITY  COEFFICIENTS 

CAL/DFG-MOLF  CAL/Ctn**Z-MULE  CAL/DfcG**3-MOLE  CAl/DEG*-4 -MOLE 

*l,4973B645E*ol  2,27523l7fefc-Ul  -»i ,817 1 1793t "0  4  5 , 692l2R3lE»0B 


- T - 

- C(PJ  - 

- - J - 

- -Ml  f  M-MtcM" 

— «tfltT)pHC299)  m~ 

— m  t  >-M(?ver 

— 08LHFt29BMM<T  >.f 

DEG  K 

CAL/UEC-MOlE 

CaL/DEG-MOLE 

KCAL/MOLF 

CAL/DEOMCLE 

KCAL/HCIE 

kcal/kolf. 

296,0 

36,167004 
-  3 8 1 4P5l 02 

aitflfiuooG 

32,136329 

4  j 1*8791 

4  j  225433 

31,660000 
-  •  31,660054 

0,000000 

,07664? 

12,600000 

12,676642 

a  0  c ,  0 

50,576881 

44,916117 

Qi697526 

33,5*6276 

4,546736 

17,146736 

500.0 

AC  1 424936 

57,301349 

14*265269 

37,066360 

10,116494 

22,716494 
- 29|169103 

- POtrnj- 

■  69 ,  Q43939 - 

41,478707 

10 , 5691 03  ’ 

7  C  G  1  0 

7* ,6*0360 

60 ,069453 

27,878448 

46,169982 

23,729657 

36,329657 

Pc  0,0 

7V,  6 P, 8 C 2 6 
-  63,614457 

90,377977 

100,009498 

35,603950 
-  43,765346 

51,059018 
-  55,968061 

31,455167 

39,636555 

44,055167 

52,236555 

- 9(J0  i'C  • 

1000,0 

87,358602 

109,027297 

52*347447 

60,826641 

46,198656 

60,798656 

1200,0 

94,057743 

125,541111 

70,462690 

70,262861 

06,333699 

78,933699 

1400,0 

102,497669 

140,634323  ■ 

90,086623 

79,250157 

65,937632 

98,537632 

1*00,0 

115,392204 

155,09P19B 

111,778785 

87,629451 

107,629994 

120,229994 

1 e  oc  ,0 

135,453767 

169,764139 

136,721320 

96,112734 

132,572529 

145,172929 

2000,0 

f65  1  393562 

~  Xfl5»4666fi6 — 

- t6d,6t8R75~ 

104 ,291594" 

162  j  47  0 1 55 

I75y070t85' 

56 


GROUP  ADDITIVITY  TH£RM0CH6P | CAL  PROPEPTIE*  PRQORIH 


CASE  NUMBER  17 

MOLECULAR  STRUCTURE 

HYPOTHETICAL  MOLECULE 

MM  H  H  H  N 

•  ■  .  H  *  H  N  ■  M  H 

H 

- - - - - - c "H — •• - <nrf - n — »«Ti — f 

m  •  m  C-»H  H  N  C  *  C  C**** 

h  C  *  m  m  m  H""C  *  •  • 

*  W  ~W '  W  -  - - - -  -rn  wmm  ■■  p-  -  ■ 

K  C  H**C  C  H  H  C  H 

•  m  *  *  ■  *  * 

C--H 

« 

- : - v - c - n - c - pt-w~k' 

H*»c  -wmm  c  C  »  ••• 

*•  H  *  «  p  •-  h  c 

H  H  C  C— h  • 

pm  ••  ••  NC»*C»*CN 

rrt - V  w  ® “  n  ■ 

...  »  H  H  M  W  • 

H  H  •  *  C 

n 

Cr*H 

»  ■  »■  ■» — h"—  0- - C**«C*«M 

•  c  •  •  ■  *■* 

H— C  H  *•  OC  *•  OC  •  H  ■ 

H 

• 

—  ,  *«  ■  *  *  m  ~  •  -*■"'■■  m  * 

m  mpm  R  0  •  -  *  •  * 

H  J- 

* 

C  *  H-OH  H 

h*“C  *  * 

■*  ..... 

M 

h  M  •  *•  •  •  ■ 

H  p  *  0®*N  C  M 

- -  - - -  H  H  —  - 

*  •  •  H  • 

H  M  H 

ATOMIC  COMPOSITION  OF  MOLECULE 
ATUH  HUMBER 

c  4b 

U  4 

. *  ~  4 


GRCUP  STRUCTURE  and  composition 


- --group  - 

- -EROL'P 

-  CHEMICAL  SYMBOL 

MR  I D  ROW 

GRID  COLUMN 

BOND 

BOND 

GROUP  NUMBER  - . 

NUMBER 

height 

CF  COMPONENT 

COORDINATE 

COORDINATE 

VECTOR 

TYPE 

OF  CORE  ATOM 

1 

296 

C 

3 

3 

1 

c 

5 

5 

4 

1 

2 

h 

1 

1 

e 

1 

0 

h 

i 

5 

2 

1 

0 

2 

399 

H - 

C 

5 

5 

1* 

0 

2 

c  ■ 
c 

7 

7 

8 

4 

...  j. 

1 

3 

H 

3 

7 

2 

1 

0 

M 

7 

3 

6 

1 

0 

3 

609 

C 

7 

7 

3 

r 

5 

5 

-  6 

"  1 

2  .  . " 

c 

4 

10 

2 

1 

4 

c 

9 

9 

4 

1 

5 

C -  ' 

'  '  10 

4 

6 

1 

6 

4 

296 

c 

4 

10 

4 

c 

H 

7 

1 

10 

A 

1 

"  1  ~ 

1 

3 

0 

H 

1 

13 

2 

1 

0 

H.  .  .  - 

- - 4 

13 

3 

1 

0 

5 

399 

C 

9 

9 

_ _ 3  _ 

- c - 

c 

"7 - 

13 

- 7  -  • 

13 

4 

- j.  - 

1 

7 

k 

12 

9 

3 

1 

0 

. 

..  H  .  — 

- 11 

7  ■ 

6 

1 

-  c 

A 

£96 

c 

ID 

4 

6 

- c - 

H 

13 

4 

- 2r 

5 

- 1  ~~ 

1 

0 

M 

13 

1 

6 

1 

0 

- - 

-  M -  '  — 

- 10 

1 

7 

1 

0 

7 

116* 

c 

13 

13 

7 

t 

9 - 

- 1! — 

- v 

3  -  - 

c 

15 

13 

4 

1 

0 

LN 

10 

16 

2 

1 

0 

" 

. . “CN -  •- 

- 15 

11 

6 

1 

0 

e 

902 

C 

15 

13 

0 

c 

13 

13“ 

0  "  " 

? 

c 

17 

17 

4 

1 

9 

c 

19 

- -  »-  — 

11 

—  *-  i  m  — 

6 

...  a- 

1 

10 

m  - - 

- ^ - 16  - 15  -  -—  5  1  o 


57 


59 


CO*PU*  UfiANfiS  1’  A*D  21  CF  COHS  ATOM  '15  ARE  EQUAL. 


-■co»*rttmioANo«— a*d— ar — efL-coRB--*Tctr--« — are  mualv 


COMPLEX  LIGANDS  45  and  46  CF  CQ«E  ATOM  34  ARE  EQUAL, 


COMPLEX  LlOANDS  41  AND  42  CF  CONE  ATCH  32  ARE  EQUAL, 


‘CDHPLTJrUDAWDf — *5  A*D  16‘'CF  CORE  ATrp  12  ARE  EQUAL, 


QRCUP  NUp.etH 
1  - 


•  1C~“ 


12 


13 


15 


AND  29  CF 

CORE 

ATCM  14 

AND  6  CF 

CORE 

ATCM  3 

-CORE  A  TCP 

and-- 

VtGANDS 

SUBGROUP 

h 

C 

eoRt 

0 

1 

LIGAND 

75 

47 

L  I  (»AND 

1 

0 

- LI  LAND- - 

— 1._ 

-—•0 . 

L  I  Ci-*  ND 

1 

0 

CDHfc 

Q 

1 

UGANP 

73 

*6 

ligand 

3 

1 

ligand 

1 

' 

0 

L  I  (j  AND 

.  -  r  - 

CORE: 

0 

1 

ligand 

5 

2 

-  LIGAND 

3 

1 

LIGAND 

67 

43 

LIGAND 

3 

1 

CORfc  ~ 

. o  - 

.....  x  ... 

ligand 

79 

47 

LIGAND 

1 ' 

0 

LfliAND 

1 

0 

ligand 

1 

0 

core 

0 

1 

VI G AND 

11 

5  ~ 

LIGAND 

69 

42 

ligand 

1 

C 

LIGAND 

1 

c 

CQRfc 

0 

1 

LIGAND 

7b 

47 

- LIGAND - 

...  -1- 

—  c . 

LIGAND 

1 

0 

ligand 

1 

0 

core 

0 

1 

ligand 

1 3 

6 

LIGAND 

65 

3  V 

LIGAND 

0 

1 

ligand 

0 

1 

CORE 

0 

1 

ligand 

13 

V 

LIGAND 

55 

34 

LtGANP 

9 

4 

* - LIGAND - 

- 1- 

- D  ”  ‘ 

CURE 

0 

1 

LIGAND 

23 

14 

LIGAND  - 

13 

7 

LIGAND 

41 

29 

LIGAND 

1 

1 

- -CORE  — 

o- 

—  J . 

ligand 

69 

44 

LIGAND 

7 

3 

-  LIGAND  - 

1 

0 

LIGAND 

1 

0 

CORE 

0 

1 

L  I GAND 

- 69 

—  41 

ligand 

11 

6 

LIGAND 

1 

C 

“  LIGAND  - 

1 

0 

.  core. 

n 

1 

LIGAND 

3V 

?3 

LlliANT 

-  3? 

?3 

ligand 

3 

1 

LIGAND 

1 

0 

core 

0 

1 

LIGAND 

77 

47 

ligand 

1 

C 

- core - 

- 1"'  — 

LIGAND 

71 

45 

LIGAND 

3 

1 

-  LIGAND 

3 

1 

ligand 

1 

0 

CORE 

0 

1 

tl  GAND 

67 

LIGAND 

3 

1 

LIGAND 

5 

3 

ligand 

3 

1 

I  n  each  group 


c 

0 

* 

0 

o 

o 

0 

4 

0 

c 

'O'  - 

0 

0 

0 

4 

c 

0  ■ 

4 

c 

c 

c 

0 

tr  - 

4 

c 

0 

c 

4 

Z  ' 
0 

0 

0 

G 

4 

c 

G  1 

c 

c 

4 

c 

c 

G 

G 

0 

2 

2 

C  - 

4 

0 

c 

0 

0 

4 

0 

0 

0 

0 

k 

-‘i  — 

c 

0 

1 

2 

1 

-C-- 

4 

s 

0 

Q 

0 

-“4— 

0 

0 

0 


In 

0 

4 

0 

■  ■  o  — 
0 
0 

4 

0 

0 

"0  ~ 

o 

G 

0 

4 

0 

C  - 

4 

0 

0 

D 

G 

'  V  " 

4 

0 

G 

G. 

4 

-  c  • 

0 

G 

C 

0 

2 

1 

1 

0 

2 

2 

0 

-o  - 

0 

2 

0 

2 

0. 

0  ■ 
4 
0 
0 
0 
0 


0 

0 

0 

0 

2 

2 - 

0 

0 

0 

4 

G 

"C  . — 

4 

0 

0 

c 

0 

4 - 

0 

0 

c  ; 


60 


16 


CORE  0  1  0  0 

LIGAND  41  25  2  2 

VtGAND - 1 - —  1 - 2 - 0 - 

LIGAND  23  14  0  2 

LIfiAND  13  7  0  0 

- *7~  - - CO«E  - 0-  j - 0  j m  0  „  . . .  . . . 

LIGAND  75  47  4  4 

LIGAND  10  0  0 

- tmt^D - 1 - 0 - 6 - 0 - 

LIGAND  1000 

1*  COflt  0010 

- - - - LIGAND - 77 - 4« - 3  -  4 -  - 

LIGAND  1  0  0  0 

19  CORE  0100 

- fcrWANB - W - 47 - 4 - 4 - 

LIGAND  1000 

ligand  1  0  0  0 

- - - LIGAND - 1 - 0 .  0  - 0  -  -  - - 

20  CORE  0  l  0  0 

LIGAND  73  45  4  4 

- fciCAftfl - 3 - 1 - 0 - 0- - 

LIGAND  2100 

21  CORE  0100 

- - LIGAND - 7> - 47 4  -  ~  4  - - - -  - 

LIGAND  1000 

LIGAND  1000 

. ‘  •  ~  •  LIGAND  J-  •  0  0  0  . 

22  CORE  0  1  0  0 

LIGAND  75  47  4  4 

- : - 1 1  liAND - 1 - 0 - 0 - 0 - 

LIGAND  1000 

ligand  1  0  0  0 

- 2* - CORE - o - 1 - 0 - 0 - - 

LIGAND  76  47  4  4 

LIGAND  1000 

.  ligand  l  '  0  0  0 

24  COHE  o  1  0  C 

LIGAND  75  47  4  4 

- LIGAND - 1 - 0 . 0  0  -  - - - - - 

LIGAND  1000 

LIGAND  1  0  0  0 

25  CORE  0  1  0  0 

LIGAND  75  47  4  4 

LIGAND  1  0  0  0 

- L1&AN0 - i—  o - o  —  o - •  . . 

LIGAND  1000 

26  CORE  0  1  1  C 

. .  -  LIGAND  -77  47  .  2  4 

ligand  1  01  0 

27  CORE  0  1  0  0 

- LTD  AND - 55 - 34“"~4 - 2 - - - 

LIGAND  13  9  0  2 

LIGAND  «400 

-  - - - LIGANO  10  0  0 

26  CORE  0100 

LIGAND  65  41  4  4 

- LIGAND - 11 - 6 - -0 - 0 - 

LIGAND  .1000 
LIGAND  1000 

— - 2f - ~  -  -  CORE  "  0  0  10 

LIGAND  77  46  i  4 

LIGAND  10.00  ... 

- frfl - CORE - 0 - 1 - 6 - C  — - - — “ 

LIGAND  65  39  4  2 

LIGAND  13  6  3  0 

■: -  ■  -  LIGAND.  0  1-  -  0  1 

LIGAND  o  1  o  1 

31  CORE  0100 

- L-fGANn - 36 - -44 - 4 - 4 - 

LIGAND  7  3  0  0  -  • 

LIGAND  1  0  0  0 

-  .....  - LIGAND  - . 1-0-0  0 

32  CORE  0  1  0  0 

LIGAND  67  42  4  4 

- tTGANtl - 5 - 3 - 1 - 0 - - - - - 

LIGAND  3  l  0  0 

LIGAND  3  1  o  0 
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33 


34 


.  CO«fc 
UliiANt) 
LIGAND 
-  tTG*  NO¬ 
LI  GA-ND 
CORE 
LtGAND  - 
LIGAND 
L [ RAND 
'TTGAPflT- 


66 

11 

1 

0 

-71 

3 

3 

0 

67 

-  3 
3 
*> 

0 

75 

1 

—  1 
1 
D 

75“ 


1 

42 

5 

0 

1 

<3 

1 

1 

~0~ 

1 

43 
1 
1 
2 
1 

47 

0 

'  O  ’ 

0 

1 

-47- 

0 

0 

0 

1 

46 


0 

0 

-  1  - 

47 

0 

““0“ 

0 

1 

4  5 

1 

1 

-  1- 

47 

0 

0 

0 

1 

-47  - 
0 
0 
0 
1 
47 

-■  O' 

0 

o 

i 

47 

G 

-0  ' 

1 

47 

0 

0 

0 

~1 

47 

0 

•  -o 

0 


35 


37 


3fl 


CORE 
LIGAND 

- “~ 

LIGAND 

LIGAND 

COKE 

LIGAND 

LIGAND 

“UtGAND - 

LIGAND 

CQHfc 

-tTtGAND - 

LIGAND 

ligand 
-ligand - 

CQRfc 
LIGAND 
-LIGAND- 


“37— 


-41— 


42 


43 


43 


LIGAND 

LIGAND 

-  CORE  - 

ligand 

ligand 

— LTGAND— 
LIGAND 
CORE 

-  LIGAND 
LIGAND 
LIGAND 

— CORE — 
LIGAND 
UGASO 

-  LIGAND 
LIGAND 

CO«E 

~t  F  GAM5 — 
LIGAND 
ligand 

-  LIGAND 

CORE 
LIGAND 
— trlUAND— 

ligand 
ligand 
CORE  ■ 
LIGAND 
LIGAND 

-ligand- 

core 

LIGAND 

LIGAND 

ligand 
ligand 
— CORE- 
LIGAND 
LIGAND 
LIGAND 
LIGAND 


1 

1 

*1- 

0 

73 

-3- 

.  1 
1 

-0 

75 

1 

— 1- 

1 

0 

73  • 

3 

2 

— ih- 

75 
1 
1 
1 
0 

•—?  S— 

1 

1 

0 

76 

1 

1 

0 

76 

1 

~1- 

D 

76 

1 

1 

1 

-0- 


73 

1 

1 

l 


~4— 

0 

0 

0 

0 

4 

•  9' 

0 

0 

0 

4 

G 

O'-  - 

0 

0 

4 

0 

0 

~o — 

4 

0 

0 

0 

0 

-4 - 

0 

0 

0 

0 

4 

o  — 
0 
0 
o 

4 

0 

-o- 

0 

4 

0 

0 

0 

-O'" 

4 

0 

0 

0 


— o 
0 
4 
0 
0 
0 
0 
4 
0 
o 
c 
0 

— 4— 

0 

0 

-  or 

0 

4 

— O' 

0 

0 

0 

4 

G 

0“ 

0 

0 

4 

0 

0 

~o- 

4 

0 

0 

0 

0 

-  4  -  ■ 

0 

0 

0 

0 

4 

— 0~ 

0 

0 

o 

4 

0 


9 

0 

0 

4 

0 

0 

0 


NUMRER  or  CORE  aTOHS  IN  LONGEST  CHAIN  a  13 


GHOLP  NLHBER  Or  CORE  CONSTITUENTS  OF  LONGEST  CHAIN  ■  1,  2’,  3,  5,  7.  6,  9,  12,  14,  27, 

30,  33,  35,  36,  39, 
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HfiHttn  OttP — tNYER-AtTTtON-  AND — — RflTAY  W#Atr-iYHH«TRt-  -  COnYRIRUTIONI - - - . 


TYPE 

GAUCHE  ALKANE 
GAUCHE  ALKANE 
-HHHiCHfc-*LHAHE — 
GAUCHE  ALKANE 
GAUCHE  ALKANE 
-OAUCHE-*t:HA*€ — 
GAUCHE  ALKANE 
GAUCHE  ALKANE 

-  gauche-alkane — 

GAUCHE  ALKANE 
GAUCHE  ALKANE 
OAUCMfc  AtrKANf! — 
GAUCHE  ALKANE 
GAUCHE  ALKAl-E 

-  GAUCHE— ALKANt" " 
GAUCHE  ALKANE 
GAUCHE  ALKANE 


INTERNAL  ROTATins 
INTERNAL  ROTATION 
-J  K  T  e«  H  A  tr~RO  T  A  T  totr- 
jktehnal  rotation 
INTERNAL  ROTATION 


- INTERNAL  RPTATtOfr- 

Internal  rotation 

INTERNAL  ROTATION 

fnttn fwm L  f^UT  "  *un 

INTERNAL  ROTATION 
INTERNAL  ROTATION 


ring  NUMBER  or 
QR0UP  NUMBER 

2  3 

3  5 


7  a 

8  9 

—a — ttf- 

9  11 

9  12 

“10  —  1R  “ 

11  19 

12  16 


AMOUNT 

2 

2 

- 2 - 


nfct A T  runn»f 

*  KCAL 
1,600 
1.600 


-1C — 27- 
16  26 
27  30 

-27  ,11 

26  32 


30  33 


*-lT606- 
3,200 
3,200 
— ,-fiOff- 
1.600 
3,200 
~  ,800 
1,600 
3,200 
~3,2C  (T 


ENTROPY  HSAT  CAPACITY  CQEFF 1C IENT8 

CAL/MI  n  ■  CAL/CIt'R  CAL/010  R**l  CAL/DEQ  K«*|  GAl'SIQ  MM 


•2.103 
•2,163 
•2,103 
•2,183 
“2,103 
-*2,103“ 
•2,103 
•2,103 
•2,103 
-2,183 
-2,103 
*?;  103"' 
•2,103 
-2,103 
•2,103 
•2,103 
•2,103 


- - INTERNAL*  ~ROTAT.!OhAL  -SYMH6TRY  C.CNTRIeUT  I  OK  - 

SOURCE  SYMMETRY  NUMBER  ENTROPY 

- CgWTg-r - 1 - — CyO  0000“ 


OPTICAL  1SQMER  CONTRIBUTION 


COHPUTEC 


12 


-rEso— — Ttrnft-crTtCTtc-!  sohers- 

4  16 


ENTROPY^ 

0,90969 


GROUP  NUMBER  OF  ASYMMETRIC  CARBON  ATOMIS)  ■ 


16,  27, 


-OWR-tfrnnER-or-pSeUDCASYHrETfllC-CAHBON  *T0*  B  12 


f 


therkcchemical  properties  or  molecule 

HfcAT  CAPACITY  COEFFICIENTS 

CAL/DEG»MOLE  CAL/CEG* •2*MULE  CAL/0tG**3»MCLE  C AL/CEO* M-MOLE 

•2,564i2399E#oi  l ,  c22Q0455£*CQ  .0 , 53203824 E*Q4  1 , 6068o3406*07 


T 

-nEO-K- 
298,0 
300 , 0 


«*0,c  • 
9CC,0 
600,0 
— 7CO.-C-  ' 
0CO.C 
900,0 


1000,0" 

1200,0 

1400,0 


C(P) 

- CAL/L'fcC-HOLE  - 

225,1*4519 

226,933717 

- 2A6T?frfeC99— 

342,200876 

3P7.195968 

-  424,9lfe09l~ 

456,333764 

482,407904 


1600,0- 

1600,0 

2000,0 


534 , 1C5B3T“ 
936,234314 
566,441797 
450704- 
639,983699 
692,763266 


CAL/DFC-MOLfl 

525,954945 

527,465711 

-6*1  i  417773 - 

671,779735 
736 ,2*5249 
00  <*,07169? 
659,727219 
915,026692 


— 907.0126, . 
1062,087460 
1147,201622 
1224,735565 
1297.136695 

.1*4  ,BC„1 


KCAL/MDLB 

3212957?6 

32,747495 

-901603328- 

90t234474 

126i76092O 

167i*3llb2 

211x542310 

298i5l9B64 

307,0 - 

412,275070 

522,777006 

630i97i873- 

761i992224 

894,515125 


•  (QI T )«H  f 290 )  )/T 
CALm-OKCLE 
525,954945 
525 ,959901 

- 535,646093-- 

955,902338 

580,61000’ 

* “607, 671116 
635,669051 
663,666572 

1.-43C300 - 

745,436060 
790,857800 
045,963025  - 
091,751739 
935,790090 


H|T)-h|298) 

KCAL/MOLE 
0,000000 
.  ,451719 

-26,307552 
57,936^98 
94,473144 
135,115406 
179,246534 
226,224100 
-275,582241  — 
379,979294 
490,481231 
606,676097 

740  ifliiiA 


DELHr<290)*KT)-H(290) 

KCAL/MOLE 

•326,703000 

•326,248201 

- - *300,39244?  - - 

•266,761302 

•232,226856 

•191,504594 

•147,453406 

•100,475892 

- •51,117759  . — 

53,279294 

103,781231 

279,970097 

402,990449 

535,519390 
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Table  4-3  (Page  1 

LIST  OF  INPUT  DATA  FOR  SAMPLE  CALCULATIONS 


1  0  0  0 
H  H 


0  0  0  0 


HHHHHHHH 


H-OC-C-C-C-C-C-C-O-H 


0  0  0  0 
H  H  H  H 


0  0  0  0 
H  H  H  H  H  H 


H — C-C - 0 - C-C — H 


Table  4-3  (Page  2  of  6) 

LIST  OF  INPUT  DATA  FOR  SAMPLE  CALCULATIONS 


H  —  C  C  —  H 

~H - B  H - R — R — H — U  - - •- - 


H — C — C — C — C--C — C — C- — C  C — H 


Table  4-3  (Page  4  of  6) 

LIST  OF  INPUT  DATA  FOR  SAMPLE  CALCULATIONS 
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Table  4-3  (Page  5  of  6) 

LIST  OF  INPUT  DATA  FOR  SAMPLE  CALCULATIONS 
0  0  0  0 

C  ....  • 


Table  4-3  (Page  6  of  6) 

LIST  OF  INPUT  DATA  FOR  SAMPLE  CALCULATIONS 


0  0  0  '  0  -  . . hypothetical  molecule 

H  H  H  H  H  H 

-  -  -  -  . _H._r.il_ . H  -  H  H  H 

c  H  —  -  H  -  - 

-  -  -  C  —  H  H  H  C  -  C  C - C  —  H 

H  C  -  -  ~  .  .  H— C - 

H  C  H— C  C  H  H  C  H 


-  H — C — -C  4 — C - C — H 

H  “  “  '  “ 

H-C-H  H  H 

H — C 

—  H  —  C-C - C - C-H 

H  -  -  C  —  H  C  H 

H  H  — 

H  H  H 
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Section  5 

LOGIC  FRAMEWORK  OF  PROGRAM 

The  TGAP  program  is  written  in  the  proposed  American  Standard  version  of 
the  FORTRAN  IV  interpretive  language.  No  machine  language  routines  nor 
nongeneralized  FORTRAN  statements  that  are  unique  to  specific  compilers 
are  employed.  The  intent,  of  course,  is  to  be  able  to  run  the  program  on 
essentially  any  computer  system  without  the  need  of  modification  or  revision. 

The  program  is  composed  of  55  routines.  These  are  listed  in  Appendix  B.  A 
complete  glossary  of  the  program  variables  is  presented  in  Appendix  C.  TGAP 
is  also  the  name  of  the  main  routine  in  the  program  which  calls  the  principal 
subprograms  into  execution.  Figure  5 -  1  gives  the  logic  flow  chart  foi  TGAP. 

There  are  three  main  divisions  in  the  logic  structure  of  the  program.  Each 
of  these  is  discussed  in  some  detail  below,  Described  herein  are  also  the 
basic  rules  and  decisions  which  govern  the  calculation  of  the  various  nonring 
interaction  corrections  and  symmetry  contributions  of  the  molecule. 

5.  1  PROCEDURE  FOR  GROUP  IDENTIFICATION 

Section  1  of  the  program  is  concerned  with  the  identification  of  the  groups  and 
group  components  as  well  as  the  group  weights,  bond  vectors,  and  bond  types. 

Subroutine  STAND  is  the  control  subprogram  for  Section  1  .  There  are  a  total 
of  1 2  routines  employed  in  this  section.  The  cross  reference  map  is  dis¬ 
played  in  Figure  5-2.  Several  routines  such  as  IDENT,  ASSIGN,  BOND, 
NUMBER,  SUMATM  perform  different  basic  functions  and  are  employed  more 
than  once. 

Subroutine  STAND  controls  the  scanning  of  the  graphic  formula  input  array 
GRID  (i.j)  of  the  molecule.  The  scan  order  is  first  rows  then  columns,  start¬ 
ing  with  the  left  uppermost  location  (1,  1)  of  the  array.  Unless  an  error  was 
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STAND 


Figure  5-2.  Subroutine  Cross  Reference  Map  for  Section  1 


committed  on  input,  the  first  nonblank  character  is  a  chemical  symbol. 
Subroutine  IDENT  determines  if  the  symbol  is  valid  and  if  it  is  composed  of  a 
single,  char acte r  or  more  than  one  character,  such  as  the  radicals  CO,  CN, 
NO  and  NO£*  Subroutine  FIND  ascertains  if  the  firs.t  atom  is  a  core  or  non- 
,  core  atom.  If  it  is  a  noncore  atom,  it  finds  the  core  atom  of  this  particular 

group.  The  identification  of  all  the  other  group  components  in  the  molecule 
is  performed  by  SCAN,  The  scan  operation  proceeds  from  group  to  group 
I  until  a  terminal  core  atom  is  reached.  If  branch  structures  are  present,  the 

scan'  operation  is  then  shifted  to  a  new  branch  which  has  not  yet  been 
examined. 

i  •  * 

The  scan  directions  are  set  by  Subroutine  ASSIGN.  Subroutine  BOND  checks 
the  validity  of  symbols  used  as  chemical  bonds  and  scans  same  until  a  dif¬ 
ferent  character  is  encountered.  In  Section  1,  Subroutine  SUMATM  computes 
the  atomic  composition  of  the  molecule  and  later  in  Section  3  the  atomic  com¬ 
position  'of  feach  ligand.  i 


t 
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The  identification  weight  for  the  core  atom  and  all  ligand  atoms  of  a  group  is 
computed  by  Subroutine  DELTA  1.  The  group  weight  is  modified  by  Subrou¬ 
tine  DELTA2  when  one  or  more  of  the  ligand  atoms  is  also  bonded  to  a  sur¬ 
rounding  carbon  atom  by  means  of  a  double  or  triple  bond,  as  for  example  in 
the  groups 


H 

\  ■ 

C  =  C  and  H  -  C  -  C  =  C 

\ 

C 

The  group  weight  is  altered  once  more  in  Section  3  of  the  program  if  aro¬ 
matic  carbon  atoms  are  present  in  the  molecule.  In  order  to  execute  this 
operation,  however,  all  ring  structures  must  already  have  been  classified. 


H 

i 

H  -  C  - 

i 

C 


The  presence  of  a  ring  structure  is  manifested  when  the  identity  flag  of  a  par¬ 
ticular  core  atom  about  to  be  scanned  registers  1000,  thus  indicating  the  atom 
has  already  been  processed  in  a  previous  cycle.  This  will  activate  the  ring 
identification  operations  in  Sections  2  and  3. 


Several  important  molecular  property  arrays  and  variables  are  generated  in 
this  section.  These  include  array  IX  (i,  j,  k)  which  contains  the  basic  struc¬ 
tural  data  described  in  3.  2.  3,  the  connectivity  array  KGN(i),  the  core  atom 
arrays  NC(i)  and  NW(i),  the  group- weight  array  MOLWT(i),  the  branch  atom 
array  IDBR(i),  the  total  number  of  branch  atoms  NOBR  and  the  total  number 
of  core  atoms  KCC. 


5.  2  PROCEDURES  FOR  CHAIN  AND  RING  IDENTIFICATION 
Section  2  of  the  program  identifies  all- the  chain  segments  and  the  components 
of  each  unique  ring  in  the  molecule.  Subroutine  CHAINM- is  the  control  sub¬ 
program  for  Section  2.  A  total  of  8  routines  are  employed  in  this  section. 
The  cross  reference  map  for  these  subroutines  is  presented  in  Figure  5-3. 

Subroutine  CHAINM  monitors  the  calculation  of  the  chain  properties.  For 
this  purpose  the  first  element  in  the  first  chain  is  set  equal  to  the  first 
terminal  core  atom  in  the  molecule.  If  no  such  atom  is  present,  each  core 
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Figure  5-3.  Subroutine  Cross  Reference  Map  for  Section  2 


atom  in  the  molecule  must  be  a  component  of  a  ring  structure.  In  this  case, 
the  first  element  is  set  equal  to  the  core  atom  of  group  number  one.  The 
core  atom  bonded  to  element  one  becomes  the  second  element  in  the  chain, 
and  so  on,  until  a  ring  closure  or  a  terminal  atom  is  reached.  All  linkage 
information  is  derived  from  the  IX  molecular  property  array. 

If  a  terminal  atom  is  reached  and  there  are  branch  junctions  in  the  molecule, 
Subroutine  NEWCOL  constructs  the  new  chain  by  first  locating  the  nearest 
branch  atom  with  unused  core  ligands.  Subroutine  NEWKC  finds  the  group 
number  of  the  core  ligand  which  is  to  replace  the  previous  chain  segment  of 
the  branch.  The  procedure  is  continued  until  all  the  core  ligands  of  every 
branch  atom  have  been  processed.  In  this  cycle  the  group  numbers  of  the 
components  of  each  chain  are  stored  in  array  NBC(i,  j),  the  locations  of  the 
branch  atoms  in  each  chain  in  array  NBX  (i,  j),  and  the  total  number  of 
components  and  branch  atoms  in  each  chain  in  array  NBS(i,  j). 


\ 
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In  the  event  the  presence  of  a  ring  was  detected  in  Section  1,  the  program  now 
checks  each  new  element  in  the  chain  for  ring  closure.  When  this  occurs, 
the  particular  chain  is  closed  and  a  new  chain  is  started  by  Subroutine  NEW- 
COL  providing,  of  course,  there  are  still  new  chain  segments  to  be  formed. 
One  advantage  of  this  procedure  is  that  the  ring  components  are  defined 
simultaneously  with  the  chain  elements.  Subroutine  RING  copies  the  ring 
components  and  ring  size  in  array  IRING(i.  j).  When  polycyclic  fused  ring 
structures  are  present,  the  same  ring  may  be  identified  more  than  once  with 
the  ring  elements  stored  in  a  different  order.  Subroutine  RING  also  checks 
for  situations  such  as  this  and  eliminates  the  duplicate  ring  structures. 

When  all  the  chain  segments  have  been  constructed,  three  data  checks  are 
performed  that  may  lead  to  a  redefinition  of  the  number  and  size  of  the  ring 
structures  and  the  length  of  the  chains: 

A.  If  more  than  one  ring  is  present,  each  ring  is  checked  by  Subroutine 
LESSEN  to  determine  whether  it  is  a  basic  cyclic  unit  and  not  a 
composite  of  two  or  more  ring  units.  All  composite  ring  structures 
are  discarded. 

B.  If  the  first  atom  in  the  first  chain  is  a  ring  atom  and  two  or  more 
rings  are  present  separated  by  a  chain  segment,  if  is  likely  the 
lengths  of  all  the  chains  can  be  increased  the  same  amount  simply 
by  shifting  the  initial  scan  location  of  the  ring  next  to  the  branch 
point.  This  is  accomplished  by  Subroutine  RESETR.  Case  1  5  in 
Section  4  is  an  example  of  this  type  of  situation. 

C.  Subroutine  CHANGE  will  also  increase  the  lengths  of  all  the  chains 
an  equal  amount  if  there  is  a  chain  that  contains  only  one  branch 
atom  and  more  "non- similar  "  atoms  than  "similar"  atoms.  The 
"similar"  atoms  are  those  atoms  listed  first  in  the  chain  and  com¬ 
mon  to  all  the  chains  in  the  molecule.  This  type  of  situation  is 
illustrated  below  with  an  example  of  a  chain  and  ring  identification 
procedure.  The  core  atom  group  numbers  that  appear  in  the  skele¬ 
ton  drawing  of  the  molecule  demonstrate  the  order  in  which  the 
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graphic  formula  is  scanned  in  Section  1.  Each  column  contains  the  group 
numbers  of  the  elements  in  the  chain.  Prior  to  the  execution  of  CHANGE 
the  "similar’1  atoms  consisted  of  atoms  1,  2,  and  3. 
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5.  3  PROCEDURES  FOR  DETERMINATION  OF  NON  GROUP  INTERACTIONS 
AND  SYMMETRY 

Section  3  of  the  program  is  a  multi-functional  logic  structure  that  computes 
the  various  nonbonded,  nongroup  molecular  interaction  corrections,  the 
contributions  for  internal  and  external  rotational  symmetry  and  the  entropy  of 
mixing  contribution  associated  with  the  asymmetry  or  optical  activity  of  the 
molecule . 

Subroutine  CORCIG  is  the  main  or  control  subprogram  of  Section  3.  A  total 
of  33  subroutines  compose  the  logic  structure  of  this  section.  The  subrou¬ 
tine  cross  reference  map  is  depicted  in  Figure  5-4.  The  logic  and  computa¬ 
tional  procedures  associated  with  the  principal  functions  of  this  section  are 
described  below. 

5.  3.  1  Gauche  Interactions 

The  program  identifies  the  gauche  alkane,  alkene  and  ether  interactions  of 
the  molecule  and  computes  the  total  contribution  of  each  to  the  heat  of 
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Figure  5-4.  Subroutine  Cross  Reference  Map  for  Section  3 


formation.  The  individual  gauche  contributions  are  small  as  shown  in 
Appendix  A.  In  highly  branched  molecules,  however,  the  sum  of  these  con¬ 
tributions  can  amount  to  several  kilocalories. 

Table  5-1  depicts  the  various  types  of  gauche  interactions  and  gives  the  rules 
used  to  establish  the  presence  or  absence  of  gauche  configurations.  Actually 
the  two  central  atoms  and  M2  in  the  gauche  alkane  and  gauche  alkene  type 
interactions  can  be  either  (a)  two  non-ring  atoms  or  (b)  one  ring  and  one  non¬ 
ring  atom0  An  example  of  the  latter  is  found  in  the  compound  cis-1- 
methyl,  2-ethyl  cyclopentane 0  In  the  drawing  below  the  two  central 
carbon  atoms  involved  in  the  interaction  are  labelled  1  and  2,  Also  refer  to 
Case  12  in  Section  4  for  the  thermochemical  property  calculations  of  this 
molecule. 

H  H 

V  / 


H  H 


Table  5-2  gives  the  number  of  gauche  interactions  associated  with  the  most 
stable  configuration  of  each  type  which  are  not  already  included  in  the  group 
value  of  M-  or  M  .  For  example,  in  the  linkage  (C)0  C- C(  C)(- C)  the  gauche 
interaction  between  the  double  bonded  carbon  ligand  and  the  ligands  of  the 
first  carbon  is  taken  into  account  by  the  group  C(C)  (C^),  in  which  repre¬ 
sents  a  carbon  atom  double  bonded  to  another  carbon  atom,  Hence,  this 
interaction  is  not  included  in  the  table*  Furthermore,  in  accord  with 
Reference  2,  the  gauche  interaction  for  (C )(H)  C- C( C)( C^)  is  ignored  due  to 
the  small  magnitude  of  this  correction*  The  central  atoms  represented  in 
the  table  are  of  the  carbon- carbon,  carbon-nit r ogen,  and  ether  types*  The 
central  atoms  are  not  included  in  the  ligand  count*  The  columns  entitled 
Type  2  give  the  number  of  ligands  bonded  to  atoms  one  and  two  which  are  of 
the  type  designated  by  Rule  2  (cf.  Table  5- IK  Those  configurations  with  a 
total  of  6  ligands  (3  and  3)  are  gauche  alkane  interactions  whereas  the  2-3 
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combination  represent  gauche  alkene  types.  In  the  ethers,  the  maximum 
number  of  interactions  is  less  since  oxygen  has  only  the  one  CO  ligand  in 
addition  to  the  CO  ligand  that  constitutes  the  central  bond. 

Table  5-1 

RULES  FOR  DETERMINATION  OF  GAUCHE  INTERACTIONS 


W 


1.  Atoms  1  and  2  must  be  a  carbon  and  oxygen  atom  of  the  indicated 
connectivities. 

2.  Ligand  U,  V,  and  W  must  be  a  heavy  (non-hydrogen)  atom.  Atom 
3  must  be  a  carbon  atom. 


79 


The  gauche  interaction  configurations  are  identified  by  Subroutine  GAUCHE. 
Each  pair  of  carbon  atoms  in  the  molecule  is  checked  fpr  the  presence  of 
gauche  formations  according  to  the  rules  of  Table  5-1.  .  Use  is  also  made  of 
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Table  5-2  to  compute  the  total  contributions  to  the  heat  of  formation. 
Subroutines  CORCIG  and  SORNGI  monitor  the  selection  of  the  non- cyclic 
atom  pair  and  the  cyclic/non-cyclic  atom  pair,  respectively. 

5.3.2  Ditertiary  Ether  Contributions 

In  addition  to  the  four  gauche  ether  contributions,  ditertiary-butyl  ethers 
require  an  additional  correction  to  the  heat  by  formation  of  8.4  kcal/mole. 

If  gauche  ether  interactions  are  present,  Subroutine  DITERE  determines 
whether  the  carbon  atom  in  every  two  adjacent  CO  ether  linkages  has  a  ter¬ 
tiary  butyl  structure  and  applies  the  appropriate  correction  to  the  heat  of 
formation. 

5.3.3  Cis  Contributions 

In  some  respects  the  identification  and  correction  of  cis.  effects  is.  similar  to 
that  of  gauche  interactions .  Each  atom  pair  is  examined  individually,  non¬ 
ring  atoms  and  non-aromatic  ring  atoms  may  be  involved,  and  the  overall 
operation  is  monitored  by  Subroutines  CORCIG  and  SORNGI,  as  before. 

There  are  however,  some  differences  in  the  application  of  these  two  types  of 
interactions  in  addition  to  the  obvious  dissimilarities  in  structure.-  In  the 
case  of  cis  configurations,  the  thermochemical  corrections  apply  to  the 
entropy  (in  specific  cases)  and  the  heat  capacity  as  well  as  to  the  heat  of 
formation.  Furthermore,  the  corrections  vary  for  certain  structural  groups. 
Lastly  the  rules  for  identifying  the  cis  interactions  are  different.  These 
rules  and  modifications  are  listed  in  Tables  5-3  and  5-4,  respectively.,  The 
standard  cis  thermochemical  corrections  are  listed  in  Appendix  A.  Table  5-4 
lists  only  those  corrections  that  are  non-standard. 


The  calculation  of  the  cis  contributions  is  performed  in  two  consecutive 
interphased  steps.  The  first  step,  which  is  executed  by  Subroutine  CIS,  tests 
for  the  presence  of  the  central  carbon  linkage  and  identifies  the  atom  pairs  on 
either  side  of  the  central  bond.  Subroutine  CISCOR  determines  whether  the 
composition  and  bonding  properties  of  each  ligand  pair  correspond  to  rule  2. 

If  the  criteria  are  satisfied,  a  search  is  initiated  for  the  special  structural 
configurations  listed  in  Table  5-4.  On  the  basis  of  these  analyses,  the 
appropriate  cis  corrections  are  selected  and  assigned  to  the  pertinent 
functions . 
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Table  5-3  ! 

RULES  FOR  DETERMINATION  OF  CIS  INTERACTIONS 


Y  Y' 


\ _ / 

(a)  Nonring  '  (b)  Ring 

1.  0^  and  must  be  carbon  atoms  with  a  connectivity  of  3(a) 
or  4(b). 

2.  X  and  X'  or  Y  and  Y1  must  also  be  carbon  atoms  with  a  connec¬ 
tivity  >3.  If  all  four  ligands  are  carbon  atoms  with  connectivi¬ 
ties  >  3,  there  are  two  cis  corrections  for  the  atom  pair. 


The  molecule  cis- 1 -methyl  -  2-ethyl  cyclopentane  also  provides  an  example 
of  a  cis  interaction  in  a  ring  structure.  *  The  central  atoms  (1  and  3)  are 
indicated  in  the  graphic  formula  of  the  molecule  presented  in  Section  5.  3.  1. 

5.3.4  Ortho  Corrections 

In  addition  to  the  three  aforementioned  types  of  nongroup  interactions,  the 
program  also  corrects  for  still  another  type  of  nonbonded  next-nearest 
neighbor  interaction  exhibited  by  aromatic  structures  -  the  ortho  interaction. 
In  the  case  of  pyridine,  there  is  a  correction  for  para  as  well  as  ortho  sub¬ 
stitution  in  the  ring. 

The  assignment  of  these  corrections  as  well  as  those  for  the  ring  structures 
(see  Sections  5.  3.  5)  is  preceded  by  a  series  of  logistic  tests  and  determina¬ 
tions  that  establish  the  identity  of  the  ring  structure  itself.  For  this  purpose, 
the  program  has  to  distinguish  between  simple  and  fused  ring  structures, 
aromatic  and  nonaromatic  cyclic  rings,  and  benzene,  pyridine,  and  higher 
nitrogen  -  containing  aromatics.  In  addition,  it  must  establish  the  size 
composition,  intra-ring  bonding  properties,  and  the  types  of  nonring  atoms 
bonded  to  the  ring.  A  few  of  these  properties,  such  as  the  size  and. 
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Table  5-4 


t-  Bu 

\ 

C 

X 

H 


t-  Bu 

\ 

C 

/ 

t-Bu 


t-Bu 

\ 

C 

c 


t-Bu 

\ 

C 

/ 

t-Bu 


Two  cis 


NON-STANDARD  CIS  CORRECTIONS  FOR 
SPECIFIC  C ONFIGURA T IONS 


f 

^Hf,  298 

| 

(kcal/mole) 

(kcal/mole) 

c 

/ 

-  c 

\ 

H 

4.  0 

t-  Bu 
\ 

/c 

c 

C 

C\ 

6.0 

C 

/ 

=C\ 

c 

8.  0 

t-Bu 

/ 

H 

t  -Bu 

\ 

H 

,10.0 

t-Bu 

/ 

— C 

.> 

12.  0 

t-Bu 

\c_ 

/ 

t-Bu 

t  -  Bu 

<c 

14.  0 

S 

(cal/0  K-mole) 

t-Bu 

/ 

“C 

\ 

t-Bu 

20.  0 

H 

C 

X 

c 

1 .  2 

interactions 

3.  0 

C-C-.C  = 

/ 

H 

C-C-C-C 

\ 

H 

-0.  6 

83 


components  of  the  individual  ring  units,  are  computed  in  Section  2  of  the 
program,  as  described  earlier.  Subroutine  FUSION  initiates  the  entire  ring 
calculation  procedure  of  Section  3.  FUSION  is  a  dual  function  routine  that 
identifies  the  set(s)  of  fused  ring  system(s)  present  in  the  molecule  and 
determines  certain  symmetry  characteristics  of  these  polycyclic  structures 
such  as  the  number  and  type  of  atoms  that  are  common  to  any  and  to  all  ring 
pairs  in  the  set.  This  information  is  also  used  later  in  the  symmetry 
calculations. 

In  order  to  maintain  the  nomenclature  of  the  structural  input  as  close  to 
standard  chemical  notation  as  possible,  no  special  symbols  are  assigned  to 
aromatic  carbon  and  nitrogen  atoms.  Instead  the  identification  of  these 
aromatic  groups  is  performed  directly  by  the  program.  Subroutine  CYCORR 
executes  the  identification  procedures.  Figure  5-5  illustrates  the  various 
different  types  of  aromatic  groups  in  question.  Once  identification  is 
achieved,  Subroutine  HEXGON  alters  the  identification  weights  of  all  the 
conjugated  aromatic  groups  in  the  molecule.  The  resultant  weights 
correspond  to  the  true  identification  weights  assigned  to  the  various 
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Figure  5-5.  Carbon  and  Nitrogen-Containing  Ring  Groups  Requiring  Special  Processing 
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aromatic  groups.  These  weight  modifications  are  performed  on  both 
monocyclic  and  polycyclic  fused  ring  systems  as  well  as  on  nonring  groups 
bonded  to  these  ring  systems. 

Subroutine  HEXGON  also  determines  the  number  of  ortho  corrections  in 
benzene  and  pyridine-like  rings  as  well  as  the  number  of  para  corrections 
in  the  latter.  Both  monocyclic  and  fused  polycyclic  ring  systems  are 
included.  In  order  to  qualify,  the  ortho  and  para  substituted  ligands  can  be 
any  atom  or  group  other  than  hydrogen  (See  Section  4,  case  14). 

5.  3.  5  Ring  Contributions 

The  ring  corrections  to  the  thermochemical  properties  of  the  molecule  are 
associated  with  the  ring  structure  itself.  In  most  cases  the  ring  correction 
for  the  heat  of  formation  equals  the  strain  energy  of  the  ring.  These  correc¬ 
tions  are  applied  to  all  ring  structures  listed  in  Appendix  A. 

Subroutine  CYCORR  assigns  the  appropriate  correction  to  single  ring  units. 
Subroutine  CRINGS  performs  this  function  for  spiropentane  and  the  six 
bicyclo  compounds  and  Subroutine  NRINGS  does  it  for  the  fused  tricyclic 
nitrogen-containing  ring.  In  each  case,  the  procedure  entails,  a  check  of 
the  size,  composition  and  bonding  characteristics  of  the  ring.  In  the  case  of 
fused  ring  systems,  the  number  of  rings,  the  properties  of  each,  and  the 
types  and  number  of  atoms  in  common  have  to  be  determined  also. 

A  few  of  the  sample  calculations  in  Section  4  involve  ring  corrections. 
Aromatic  systems  appear  to  be  exceptions.  This  is  because  in  aromatic 
systems  the  ring  correction  is  already  included  in  the  group  thermochemical 
contributions. 

5.3.6  Symmetry  Logic  Structure:  Part  I 

The  logic  schemes  for  the  determination  of  the  symmetry  of  the  molecule  and 
of  the  individual  groups  thereof  are  the  most  extensive  and  complex  of  the 
logic  structures  of  the  program.  These  symmetry  operations  are  divided  into 
two  parts.  Part  I,  which  is  discussed  in  this  section,  performs  the  basic 
internal  group  by  group  symmetry  evaluations  of  the  molecule.  In  particular, 
the  ligands  of  each  group  are  compared  one  versus  another  and  the  similarity 
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or  dissimilarity  of  their  structures  is  determined.  Here  as  before,  the  term 
ligand  denotes  all  the  atoms,  one  or  more,  which  form  the  bond  chain  or 
"branch”  of  a  particular  core  atom.  A  tetra,hedral  carbon  atom,  for  example, 
has  four  ligands  each  composed  of  one,  two  or  any  number  of  atoms. 

The  principal  steps  in  the  logic  framework  of  Part  I  are  outlined  in  Figure  5-6. 
The  figure  lists  only  the  control  routines  for  the  indicated  steps.  The  names 
of  the  auxiliary  subprograms  are  given  in  Figure  5-4. 

The  symmetry  calculations  of  Part  I  are  initiated  and  controlled  by  Subrou¬ 
tine  SYMTRY  The  atomic  composition  calculations  performed  by  SYMTRY 
and  its  auxiliary  routines  SCANCH  and  SCANBR  serve  as  a  screening  medium. 
In  this  manner,  complex  ligand  pairs  with  dissimilar  atomic  compositions  are 
from  the  start  automatically  eliminated  from  the  more  time  consuming 
structural  identity  tests.  Groups  in  which  the  core  atom  of  the  group  is  also 
a  ring  atom  are  likewise  exempt  from  undergoing  additional  symmetry  anal¬ 
ysis  in  Part  I, 

The  structural  similarity  of  two  or  more  ligands  which  contain  no  core  atoms, 
such  as  H,  CN  etc.  ,  is  determined  at  once  by  directly  comparing  the  atomic 
composition  of  the  ligands.  If  the  ligands  are  complex,  that  is  they  contain 
core  atoms,  the  structural  similarities  of  each  ligand  pair  must  be  deter¬ 
mined  bond  by  bond  and  branch  by  branch.  A  pair  of  ligands  is  not  consid¬ 
ered  dissimilar  until  all  possible  branch  pair  combinations  of  the  two  ligands 
have  been  tested.  This  operation  is  performed  by  Subroutine  EQUAL,  If 
ring  structures  are  encountered  during  the  structural  analysis  of  a  ligand 
pair,  the  equality  of  the  ring  structures  is  determined  by  Subroutine  EQUALR 
and  by  auxiliary  routines  FIRSTR  and  DELETE.  The  program  can  actually 
handle  several  ring  formations  per  ligand.  The  limits  imposed  on  the  various 
variables  are  those  for  ring  formations  in  general,  namely,  30  components 
per  ring  and  a  total  of  40  ring  units  per  molecule. 

Once  the  structures  of  a  ligand  pair  are  found  to  be  similar,  several  blocks 
of  symmetry  identif ication  data  pertinent  to  the  particular  ligands  and  groups 
are  generated  and  stored  by  Subroutine  SAME.  These  comprise:  (a)  the 
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Figure  5-6.  Main  Logic  Flow  Chart  for  Symmetry  Calculations:  Part  I 


group  numbers  of  the  atoms  in  the  ligand  pair(s)  which  are  bonded  directly 
to  the  core  atom*  (KCSAME(i,  j));  (b)  the  number  of  elements  in  (a),  the 
number  of  different  ligand  pairs,  and  the  code  type  assigned  to  each  similar 
ligand  (NOSAME(i,  j));  and  (c)  the  total  number  of  similar  ligands  of  each 
type  in  the  group  (NTOTAL(i,  j)). 


5*307  Symmetry  Logic  Structure:  Part  II 

The  symmetry  data  calculated  in  Part  I  are  utilized  in.  Part  II  to  compute 
the  various  molecular  properties  listed  in  Figure  5-7.  Each  of  these 
properties  is  discussed  below. 


5.  3.  7.  1  Identification  of  Asymmetric  Atoms 

The  first  function  in  the  final  sequence  of  symmetry  calculations  consists  of 
determining  the  total  number  of  asymmetric  atoms  in  the  molecule  and 
the  group  numbers  of  these  atoms.  The  following  are  the  criteria  employed 
by  the  program  for  the  selection  of  asymmetric  atoms. 

A.  The  core,  atom  must  be  a  carbon  atom 

B.  It  must  be  a  nonring  atom 

C.  It  must  have  a  connectivity  equal  to  4,  i.  e.  ,  tetrahedral. 

D.  The  four  ligands  of  the  core  atom  must  be  dissimilar 

These  tests  are  conducted  by  Subroutine  ASYMC.  If  all  the  aforementioned 
algorithms  are  satisfied  by  a  particular  atom,  the  asymmetry  counter 
NASYMC  is  incremented  by  one  and  the  group  number  of  the  atom  is  stored 
in  array  IOPATM. 

The  data  are  utilized  later  to  compute  the  entropy  effects  arising  from  opti¬ 
cal  isomerism.  These  data  are  also  used  in  the  computation  of  the  longest 
chain  which  is  described  next. 

5,  3.  7.  2  Determination  of  the  Longest  Chain 

The  calculation  of  the  longest  chain  is  a  prerequisite  to  the  computation  of 
the  external  symmetry  number.  If  asymmetric  atoms  are  not  present,  the 
longest  chain  is  simply  taken  as  the  longest  of  the  chains  in  the  molecular 
chain  array  which  was  identified  and  constructed  in  Section  2  of  the  program. 
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When  there  is  more  than  one  chain  of  maximum  length,  the  one  located  first 
in  the  chain  array  is  selected  as  the  longest  chain  of  the  molecule. 

In  the  event  asymmetric  atoms  are  present  and  there  is  more  than  one  chain 
of  maximum  length  available,  the  program^computes  the  number  of  asym¬ 
metric  atoms  present  in  the  first  chain  of  maximum  length  identified.  If  all 
the  asymmetric  atoms  are  contained  in  this  chain,  no  redefinition  is  neces¬ 
sary.  If  only  some  of  these  atoms  are  contained  therein,  all  other  chains 
equal  to  the  longest  chain  in  size  are  examined  for  their  content  of  asym¬ 
metric  atoms.  That  chain  of  maximum  length  containing  the  greatest  number 
of  asymmetric  atoms  is  designated  the  longest  chain  in  the  molecule. 

The  above  tests  are  executed  by  Subroutine  MAXCHN.  The  number  of  the 
components  and  the  chain  number  of  the  longest  chain  are  stored  in  array 
MBS,  and  the  group  numbers  of  the  chain  components  in  array  MBC.  These 
data  are  also  output  by  the  routine. 

5.  3.  7.  3  Determination  of  the  Symmetry  Number 

The  computation  of  the  symmetry  number,  a  property  that  characterizes  the 
external  rotational  symmetry  of  the  molecule,  is  a  multi-phase  logic- 
controlled  operation  executed  by  Subroutine  EXTROT.  The  principal  steps 
involve  finding  the  central  atom  or  atoms  in  the  molecule  and  the  number  of 
structurally  identical  ligands  bonded  to  this  atom  or  atoms. 

For  the  purpose  of  this  analysis,  the  central  atom  is  defined  as  the  atom  in 
the  center  of  the  longest  chain.  If  the  chain  contains  an  odd  number  of 
atoms,  the  molecule  has  one  central  atom.  If  the  number  of  chain  constitu¬ 
ents  is  even,  there  are  two  central  atoms  in  the  molecule  whose  ligands  have 
to  be  analyzed  for  similarity  of  structure. 

Once  the  even  or  odd  multiplicity  df  the  chain  is  established  and  with  it  the 
number  of  central  atoms,  additional  classifications  are  made  to  determine 
the  symmetry  number  of  the  molecule.  These  new  classifications  include 
the  type  and  connectivity  of  the  central  atom(s)  and  the  symmetry  charac¬ 
teristics  of  the  ligands.  The  principal  algorithms  of  the  classification 
scheme  are  listed  below.  The  symmetry  number  of  configurations  not 
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listed  therein  is  set  equal  to  one.  The  symbol  Cp  refers  to  the  axial 
rotational  symmetry  of  the  ligand. 

Rules  for  the  Assignment  of  the  Symmetry  Number 
A.  Odd  Number  of  Atoms  in  the  Longest  Chain 
1.  Central  Atom  is  a  Carbon  Atom 
a.  Two  Ligands 
(1)  X  -  C  =  Y 

(a)  X  and  Y  are  both  linear,  <r  =  1 

(b)  X  is  linear  and  Y  is  symmetrical  (or  vice  versa), 
<r  =  symmetry  of  Y 

(c)  X  and  Y  both  have  the  same  axial  rotational 
symmetry  C^, 

cr  =  C 


(2) 


X  =  C  =  Y 
(a)  X  =  Y 


<y  =  2  x  C 


(b)  X  ^  Y,  same  as(l) 
Three  Ligands 


/\ 


c . 


X  =  Y  and  Z  is  linear  or  has  C^  =  C2>  <r  -  2 
Four  Ligands 


(1) 

CX2YZ, 

<T  —  1 

(2) 

cx2y2, 

cr  =  2 

(3) 

CX  Y, 

cr  =  3 

(4) 

CX„  and  X  is  linear  or  has  C  =  Cv  <r 
4  P  ° 

CX4  and  X  not  linear  or  C^  4  Cy  cr  =  6 
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2.  Central  Atom  is  an  Oxygen  Atom 


/°\ 

X  Y 

X  =  Y  a=2 


3.  Central  Atom  is  a  Nitrogen  Atom 

a.  Two  Ligands 

X  -  N  =  Y  <7  =  1 

b.  Three  Ligands 
NXYZ 

X  =  Y  =  Z,  cr  =  3 

Even  Number  of  Atoms  in  the  Longest  Chain 

1.  Two  Central  Atoms  are  Unique:  <7=  i 

2.  Two  Central  Atoms  are  Carbon  Atoms 

a.  Two  Ligands 

x  =  c  -  C  =  Y  or  X  =  C  =  C  =  Y  Same  as  1.  2. 

b.  Three  Ligands 

Y\, 


(1) 


(2) 


X 


C  2  z 


cr  =  1 


\ 

/C  ~  C  =  Z  X  ~  Y  and  Z  is  linear  or  has 

C  =  CX,  <r  =  2 
P  2 


X 


F our  Ligands 

(i) 


*x 


(2) 


W  =  C  -  C  — Y 

^Z 

X  =  Y  =  Z  and  W  is  linear  or  has  C  =  C_,  <x  =  3 

P  3 


X> 


W 


,C  =  C* 


-Y  >  . 
'Z 


(a)  W  =  X  =  Y  =  Z  <r  =  4 

(b)  W  =  Z  and  X  =  Y  =  2 

(c)  W  =  X  and  Y  =  Z  <7  =2 
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3. 


4. 


(3>  X^c_c^Y 


W' 


X  =  Y  and  W  =  Z  (r  =  2 


d.  Five  Ligands 


- C'—Y 

V*^  z 


o-  =  1 


e.  Six  Ligands 


V — C  — C— Y 

(1)  U  =  V  =  W  =  X  =  Y  =  Z  0-  =  6 

(2)  U  =  V  =  Y  =  Z  and  W  =  X  cr  =  2 

(3)  U  =  V  =  W  and  X  =  Y  =  Z  *  =  1 

(4)  W  =  X,  V  =  Y,  U  =  Z 

Two  symmetry  configurations  are  possible 
(a)  Dissymmetric,  nonasymmetric  form 


(b)  Nondissymmetric  form 
W 

C 

// 

Vu/ 

Two  Central  Atoms  are  Oxygen  Atoms 


X  -  O  -  O  -  Y  (assuming  nonplanar  structure) 

X  =  Y  o-  =  2 

Two  Central  Atoms  are  Nitrogen  Atoms 
a.  Two  Ligands 
X  -  N  =  N  -  Y 

X  =  Y  <r  =  2 
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As  indicated  in  Section  2.2,  the  program  does  not  determine  the  symmetry- 
number  of  a  molecule  whose  central  atom(s)  is(are)  part  of  a  ring  structure 
unless  the  ring  is  a  monocyclic  benzene-type  ring.  In  the  latter  case,  the 
symmetry  number  assignments  are  made  according  to  the  following  rules: 


U 


Condition  & 

U  =  V,  X  =  Y,  W  =  Z  _  ~ 

U  =  W,  X  =  Z,  C  =  C,  for  V,  Y  2 

'  *  p  2 

U  =  W,  X  =  Z,  V  =  Y,  C  =  C,  for  V,  Y  2 

p  2 

U  =  X,  V  =  Y,  W  =  Z  2 

U  =  V  =  W,  X  =  Z,  C  =  C  for  V,  Y  2 

U  =  V  =  W,  X  =  Y  =  Z,  Cp  =  C2  for  all  ligands  2 

U=W=Y,  V=X=Z  3 

U  =  W  =  Y,  V  =  X  =  Z,  Cp  =  C2  for  all  ligands  6 

U  =  V  =  W  =  X,  Y=Z  2 

U  =  V  =  W  =  Y,  X  =  Z,  Cp  =  C2  for  V  2 

U  =  V  =  X  =  Y,  C  =  C2  for  W,  Z  2 

U  =  V  =  X  =  Y,  W=Z  2 

U  =  V  =  X  =  Y,  W  =  Z,  Cp  =  C2  for  W,  Z  4 

U  =  V  =  W  =  X  =  Y,  Cp  =  C2  for  all  ligands  2 

U=V  =  W  =  X  =  Y=  Z  6 

U=V=W=X=Y  =  Z,  C  =  C_  for  all  ligands  12 

_ P  _ _ . _ _ _ 
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5.  3.  7.  4  Calculation  of  Internal  Rotation  Contributions 

The  program  computes  the  contributions  to  the  entropy  associated  with  the 
internal  rotation  of  one  or  more  groups  of  atoms  in  the  molecule.  The  rotation 
about  the  bond  axis  can  be  either  of  C^  or  C£  symmetry.  Specifically,  the 
program  applies  this  correction  to  tetrahedral  carbon  atoms  of  C^  symmetry 
aromatic  rings  (benzene  and  pyridine  types)  of  symmetry,  and  -NC^  groups. 

In  the  case  of  tetrahedral  atoms,  the  rotation  can  involve  a  simple  structure 
like  -CH^  or  groups  in  which  the  structures  of  the  identical  ligands  are  highly 
complex.  (For  example,  carbon  group  12  of  case  calculation  17,  Section  4.  ) 
Subroutine  INTROT  looks  for  core  atoms  which  are  carbon  atoms  with  a 
connectivity  of  four  and  three  structurally  identical  ligands.  Each  of  these 
groups  is  assigned  a  threefold  rotational  axis  and  an  entropy,  correction 
equal  to  -R  fn3.  If  the  particular  internal  rotation  has  already  been  included 
in  the  external  rotational  symmetry  of  the  molecule,  it  is  eliminated  from 
this  calculation. 

The  auxiliary  routine  C  TWO  assigns  a  twofold  axis  of  rotation  and  an  entropy 
correction  of  -R  in  2  to  NO^  groups  and  to  aromatic  rings  of  C ^  symmetry 
providing  each  of  these  structures  is  bonded  to  a  nonlinear  ligand.  The 
linearity  of  the  ligand  is  determined  by  Subroutine  LINEAR.  The  twofold 
rotational  symmetry  of  a  ring  is  established  by  Subroutine  EQUAL  and  its 
auxiliary  subprograms.  This  calculation  is  feasible  even  when  there  are 
complex  ligand  subroutines  present  in  the  structure  of  the  ring. 

5.  3.  7.  5  Calculation  of  Optical  Isomer  Contributions 

In  order  for  a  molecule  to  be  optically  active,  it  must  in  general  be  dis¬ 
symmetric,  that  is,  it  must  lack  an  alternating  (Sn)  axis  of  symmetry  (5). 

From  the  intramolecular  point  of  view,  the  most  frequently  encountered 
albeit  not  necessary,  condition  for  optical  activity  is  the  presence  of  an 
asymmetric  carbon  atom.  The  criteria  used  by  the  program  to  identify  the 
presence  of  asymmetric  carbon  atoms  were  presented  in  Section  5.  3.  7.  1. 

Subroutine  ENTSYM  computes  the  optical  properties  (number  of  enantiomers 
and  me  so  structures)  of  a  molecule  in  which  the  optical  activity  is  due  to  the 
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presence  nonring  asymmetric  carbon  atoms.  The  number  of  active 
conformations  or  enantiomers  (a)  and  the  number  of  meso  structures  (m)  is 
computed  directly  from  equations  A,  B,  or  C  (6). 

A.  All  asymmetric  carbon  atoms  are  different: 

a  =  2^  m  =  o 

B.  Like  asymmetric  carbon  atoms  are  present  and  V  is  even: 


C.  Like  asymmetric  carbon  atoms  are  present  and  V  is  add: 


where  V  is  the  total  number  of  asymmetric  carbon  atoms  including  the 
pseudoasymmetric  carbon  atom  (if  any).  The  latter  is  defined  as  follows; 

If  two  like  asymmetric  carbon  atoms  are  joined  to  a  carbon  atom  attached 
to  two  different  nonasymmetric  groups,  two  active  forms  and  two  meso 
forms  are  possible.  The  central  atom  in  such  compounds  is  called  a 
pseudoasymmetric  carbon  atom.  The  pseudoasymmetric  carbon  atom  and 
meso  structures  are  identified  by  Subroutine  EXTROT.  Only  that  asymmetric 
carbon  atoms  contained  in  a  single  chain  of  the  molecule  can  be  included  in 
the  above  equations.  The  program  takes  into  account  all  asymmetric  carbon 
atoms  contained  in  the  longest  chain  where,  if  necessary,  the  latter  is  rede¬ 
fined  in  the  manner  described  in  Section  5.  3.  7.  2. 

Subroutine  ENTSYM  also  computes  the  enthropy  contribution  due  to  external 
rotational  symmetry 
S  =  -R  lu  (<r) 
and  optical  isomerism 
S  =  R  In  (n) 

where  n  =  a  +  m.  The  assumption  is  made  here  that  all  isomers  are  present 
in  equal  amounts. 
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ABSTRACT 


Group  additivity  thermodynamic  methods  are  applied  to  the 
calculation  df  Arrhenius  parameters.  Studies  of  several  classes 
of  unimolecular  reactions  are  presented  together  with  preliminary 
methods  developed  for  the  calculation  of  kinetic  rate  constants. 

The  problems  associated  with  bimolecular  reactions  are  reviewed 
and  some  preliminary  solutions  suggested.  Where  required  to 
support  this  work,  thermodynamic  functions  of  radical-containing 
groups,  single  bond  strengths,  and  pi  bond  strengths  for  various 
compounds  are  calculated.  A  computer  scheme  to  calculate  the 
rate  constants  for  simple  unimolecular  fission  is  presented. 
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Section  1 
INTRODUCTION 

Although  the  advantages  of  applying  chemical  kinetics  to  the  solution  of  many 
physico-chemical  problems  are  well  recognized,  the  measurement  of  basic 
kinetic  data  is  tedious  and  time-consuming.  In  addition,  the  sample  size  and 
sample  purity  required  for  such  measurements  are  often  prohibitive.  This 
makes  it  especially  difficult  to  evaluate  new  and  potentially  useful  chemical 
compounds.  A  reliable  method  for  estimating  the  kinetic  parameters  of  as 
yet  un-studied  compounds  would  be  particularly  valuable. 

Systematic  exploration  of  current  data  on  gas  phase  reactions  by  Benson^^ 
has  revealed  that  for  some  reaction  categories,  the  Arrhenius  parameters 
can  be  fitted  with  the  aid  of  some  simple  rules  into  the  frame  work  of  transi¬ 
tion  state  theories.  These  rules  involve  in  part  the  extension  of  the  methods 
of  group  additivity  to  the  calculation  of  kinetic  parameters.  The  purpose  of 
the  present  work  was  to  investigate  more  thoroughly  and  improve  upon  those 
reaction  categories  to  which  these  methods  have  been  applied  and  to  extend 
the  group  additivity  methods  to  other  reaction  categories. 

The  first  task  of  the  kinetic  investigation  was  to  examine  critically  empirical 
methods  of  estimating  activation  energies  for  elementary  chemical  reactions. 
This  has  proved  to  be  a  fairly  involved  task  because  the  problems  associated 
with  different  reaction  categories  are  very  specific  for  each  category.  Before 
a  systematic  set  of  rules  for  the  computation  of  kinetic  parameters  can  be 
formulated  and  ultimately  adopted  to  computer  programming  methods,  it  is 
necessary  to  understand  fully  the  problems  and  techniques  involved  in  calcu¬ 
lating  activation  energies.  Therefore,  the  majority  of  work  on  this  project  was 
devoted  to  the  first  task  of  estimating  activation  energies. 

Several  specific  classes  of  reactions  were  investigated  in  depth.  This  approach 
was  used  rather  than  a  general  survey  involving  less  detailed  analysis  of  a 
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large  number  of  reaction  classes  in  order  to  understand  more  fully  the  problems 
associated  with  empirical  methods  of  estimating  activation  energies. 


The  ultimate  goal  of  this  project  is  to  develop  uniform  techniques  to 
calculate  rate  constants  of  elementary  chemical  reactions  rapidly. 
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Section  2 

METHODS  FOR  ESTIMATING  THE  ACTIVATION  ENERGY  AND 
FREQUENCY  FACTOR  OF  GAS  PHASE  CHEMICAL  REACTIONS 


where  E  is  the  activation  energy,  R  is  the  gas  constant,  T  is  the  reaction 
temperature  in  °K,  and  A  is  the  frequency  factor.  Transition  state  theory 
postulates  the  participation  of  an  ’/activated  complex11  or  excited  intermediate 
in  chemical  reactions.  For  the  simple  unimolecular  reaction  this  is  repre¬ 
sented  by 


A  —  —  B 

where  A^  is  the  excited  intermediate  of  the  transition  state.  The  frequency 
factor  A,  as  given  by  the  transition  state  theory  is 

A  =  exp  (AS^/R) 

where  Ast  is  the  entropy  of  activation,  that  is,  the  change  in  entropy  in  going 
from  reactant  A  to  the  transition  state  A^.  The  other  terms  have  their  usual 
significance.  Similarly  the  activation  energy  E  is  equal  to 

E  =  AH*  +  RT 

where  AH^  is  the  enthalpy  of  activation,  namely,  the  difference  in  enthalpy 
between  the  transition  state  and  reactant.  These  equations  can  easily  be 
generalized  to  more  complex  elementary  reactions  with  more  than  one  reactant 
or  product.  In  any  case,  the  values  of  AS^  and  AH*^  are  required  to  calculate 
the  rate  constant. 
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(la) 

Benson  has  examined  the  experimentally  determined  Arrhenius  parameters 
for  the  most  common  reaction  categories  and  has  derived  fairly  simple  rules 
for  estimating  entropies  of  activation  through  consideration  of  the  changes  in 
molecular  structure  in  going  from  reactants  to  transition  state.  The  rules 
for  computing  the  A  factors  are  uniquely  defined  for  each  particular  transition 
state  model,  and  are  based  on  changes  in  bond  lengths,  bond  bending  and 
stretching  frequencies,  and  barrier  heights  to  internal  rotation.  The  method 
used  by  Benson  is  fairly  general  and  can  be  applied  with  success  to  a  wide 
variety  of  reaction  categories. 

In  contrast,  the  problems  encountered  with  the  calculation  of  AH*f  for  different 
reaction  classes  has  proved  to  be  specific,  or  nearly  specific,  for  each 
reaction  class.  In  practice,  only  those  classes  of  reaction  were  chosen  for 
the  study  of  methods  for  calculating  AH^,  for  which  methods  have  been  well 
developed  for  calculating  AS^.  There  are  several  reasons  for  this  approach. 
First,  the  transition  state  model  adopted  for  these  reaction  classes  has  been 
shown  to  be  accurate  due  to  agreement  of  calculated  and  experimental  A 
factors.  Second,  for  these  classes,  in  general,  a  large  amount  of  experimental 
data  is  available  for  comparison  of  results.  Finally,  the  methods  for  calcu¬ 
lating  AH*f  can  be  combined  with  those  for  AS*f  to  generate  reaction  rate 
constants  for  those  classes. 

The  methods  applied  to  several  classes  of  unimolecular  and  bimolecular 
reactions  will  be  discussed  in  this  section.  In  addition,  the  work  performed 
to  determine  some  bond  energies  and  heats  of  formation  required  for  this 
work  will  also  be  discussed. 

2.  1  CALCULATION  METHOD  FOR  AHf°  OF  RADICAL.  GROUPS,  BOND 
STRENGTHS  IN  SINGLE  BONDS  AND  Pi  BONDS 

After  the  onset  of  this  work,  it  was  evident  that  group  additivity  thermodynamic 

functions  were  required  for  radical  containing  groups.  The  need  for  these  values 

values  arose  several  times  during  the  course  of  this  work.  Benson  and 

O’Neal^  have  given  a  partial  tabulation,  but  additional  data  were  required. 

In  addition,  their  tabulation  does  not  take  into  account  delocalization  in  radical 

groups  which  possess  multiple-bonded  ligands.  That  is,  Benson  and  O’Neal's 

formulation  would  give  the  same  group  value  for  the  heat  of  formation  of  the 
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group  [•  C-IC^MH^)]  as  for  [’C-tCHt^)]  .  In  their  calculation,  a  second  step 
is  required  to  account  for  the  delocalization  of  the  unpaired  electron  by  the 
double  bonded  carbon.  The  present  calculation  incorporates  the  effect  of 
delocalization  into  the  group  values.  This  is  compatible  with  group  values  for 
non-radical  groups  and  removes  an  additional  step  from  the  calculations.  This 
is  advantageous  for  machine  calculations. 

A  second  requirement  for  these  calculations  was  a  tabulation  of  bond  strengths, 
Benson^k)  has  given  a  fairly  complete  tabulation,  but  additional  data  were 
required.  For  example,  previously  in  group  additivity  calculations,  the  same 
value  of  bond  energy  was  assigned  to  a  carbon-carbon  single  bond  in  a  hydro¬ 
carbon  chain  as  to  a  C-C  bond  alpha  to  a  carbon-oxygen  single  bond.  It  was 
found  to  be  necessary  to  differentiate  between  two  cases  such  as  this  in  order 
to  accurately  account  for  differences  in  activation  energy.  In  addition,  pi  bond 
energies  of  C==C  and  C=0  double  bonds  were  required.  It  was  possible  to 
calculate  all  required  bond  strengths  using  group  additivity  thermodynamic 
methods  and  group  values  of  radical  groups. 

2.  1.  1  Method  of  Calculation  of  AHf°  of  Radical  Groups 
The  formula 

AH°f  [•  R]  =  DH°  [C-H]  +  AHf°[HR]  -  AHf °  | •  H ] 

is  used  to  calculate  the  heat  of  formation  AH^°  of  the  radical  group  -R. 

DH° [C-H]  represents  the  bond  strength  of  the  carbon-hydrogen  bond  of  the 
saturated  group  HR  which  on  dissociation  yields  the  radical  *R.  AHf°  [ *  H 1 
is  the  heat  of  formation  of  a  hydrogen  atom,  namely,  52.  1  kcal/mole,  All 
heats  of  formation  are  for  the  standard  state  at  298°K. 

It  is  necessary  to  use  the  value  of  DH°  [C-H]  which  best  reflects  the  ligand 
environment  of  the  group  HR.  For  example,  the  difference  in  the  C-H  bond 
strength  for  the  CH^  group  of  C^Hg  and  CH ^=CH^CH  is  10.  5  kcal/mole. 
Therefore,  the  approximation  will  be  used  that  the  bond  strengths  given  by 
Benson  for  saturated  hydrocarbons  can  be  corrected  for  the  presence  of  a  C^ 
ligand  in  the  HR  group  by  subtracting  10.  5  kcal/mole.  Similarly  the  C-H  bond 
strength  for  the  CH^  group  of  propane  is  98  kcal/mole  compared  to  the  value 
for  acetone  of  98.  3  kcal/mole.  A  value  of  0.  3  kcal/mole  is  thus  added  to  the 
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C-H  bond  strengths  in  saturated  hydrocarbons  to  correct  for  the  replacement 
of  a  carbon  ligand  by  a  carbonyl  ligand.  Finally,  a  value  for  the  C-H  bond  of 
the  group  of  propyne  of  91.  5  kcal/mole  and  of  methanol  of  95.  9  kcal/mole 

gives  correction  factors  of  -6.  5  kcal/mole  for  ligands  and  -2.  1  kcal/mole 
for  single  bond  oxygen  ligands  respectively. 

Table  2-1  contains  the  AHf°  for  a  series  of  radical  groups  *R.  Included  in 
the  table  are  the  values  of  AH^0  [HRl  used  as  obtained  from  Benson'1, 
and  the  C-H  bond  strength  used.  Also  included  is  the  model  compound  used 
for  the  bond  strength  along  with  the  resultant  correction  factors.  Where  two 
or  more  ligands  are  present  which  require  a  correction  factor,  the  correction 
factors  are  considered  additive. 


Table  2-1 

HEAT  OF  FORMATION  OF  RADICAL  GROUPS 


•R 

AHf°  t-R] 
kcal/ mole 

AHf°  [HR] 
kcal/ mole 

DH°  [C-H] 

kcal/mole 

Model 

Correction 

kcal/mole 

•C-(C)(H2) 

35.  8 

-10.08 

98 

H-CH2CH3 

- 

•C-(Cd)(H2) 

25.  5 

-9. 95* 

87.  5 

H-CH2CHCH2 

- 

•C-(CO)(H2) 

36.  1 

-10. 08 

98.  3 

H-CH2CH3 

+0.  3 

•C-(C2)(H) 

37.  5 

-4.  95 

94.  5 

H-CH(CH3)2 

- 

•C-(Cd)(C)(H) 

27.  1 

-4.  76 

84.  0 

H-CH(CH3)2 

-10.  5 

•  C-(  CO)  (C)  (H) 

37.  7 

-5.  0 

94.  8 

H-CH(CH3)2 

+  0.  3 

•C-(C3) 

37.  0 

-1.  90 

91 

h-c(ch3)3 

.  - 

•C-(Cd)(C2) 

26.  9 

-1. 48 

80.  5 

H-C(CH3)3 

-10.  5 

-C-(H2)(0) 

33.  7 

-10. 08 

95.  9 

H-CH2CH3 

-2.  1 

*C -(H)  (0(0) 

31.  8 

-8.  5 

92.  4 

H-CH(CH3)2 

-2.  1 

•C-(0)2(C) 

17.  5 

-17.  2 

86.  8 

H-C(CH3)3 

M  a  \ 

-4.  2 

-  /  i  \ 

*The  value  of  Hr  [  ]  was  not  available  from  Benson  .  This 

i  a  J 


value  was  calculated  from  the  AH^°  =  4.  9  kcal/mole  for  CH^CHCH^. 
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2.  1.  2  Method  of  Calculation  of  Single  Bond  Strengths 

Various  C-C,  C-O,  C-H,  and  O-H  bond  strengths  were  required  to  compute 
the  energies  of  activation.  The  general  values  of  these  bond  strengths  were 
not  accurate  enough  since  they  do  not  take  into  account  the  effect  of  nearest 
neighbor  environments  in  the  compounds  where  these  groups  occur.  Some  of 
these  bond  values  were  already  available  but  most  values  had  to  be  calculated 
using  group  additivity  methods.  For  example,  for  a  C-C  bond  of  the  type 
C-CH— CH2~0(where  the  longer  bond  indicates  the  one  to  be  broken)  no  value 
for  the  bond  strength  could  be  obtained.  The  model  compound  CH^Ct^CH^OH 
was  chosen  and  the  bond  energy  calculated  from  the  reaction 

AE  +  CH3CH2CH2OH— CH3CH2-  +  •  CH2OH 

Group  additivity  thermodynamic  calculations  are  applicable  to  this  problem, 
and  the  AH^0  of  the  radical  groups  previously  derived  can  be  used  for  this 
purpose.  Two  examples  are  given  in  Table  2-2. 


Table  2-2 

SINGLE  BOND  STRENGTHS 


Bond  Type 

Bond  Environment 

Model  Compound 

Bond  Strength 
kcal/ mole 

C-O 

C-O - 

-  ch2.o 

CH3OCH2OH 

88.  0 

C-C 

c-ch2- 

ch2-o 

CH3CH2CH2OH 

83.  0 

2.  1.  3  Method  for  Calculation  of  Pi  Bond  Strengths 

Values  of  the  pi  bond  strength  for  various  C=0  and  C=C  bonds  are  required. 
For  symmetric  compounds  such  as  ethylene,  the  pi  bond  strength  is  defined 
as  the  difference  of  the  C-H  bond  strength  in  the  cleavage  of  ethane. 

i.  e.  E7r[CH2CH2l  =  DH  [CH3CH2-H]  -  D  [CH2CH2-H] 
E7r[CH2CH2]  =  98  -  38.  5  =  57.7  kcal/mole 

Experimental  values  are  available  for  this  example,  but  group  additivity 
thermodynamic  methods  along  with  radical  group  values  have  been  used  to 
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calculate  the  pi  bond  energy  of  those  compounds  which  lack  such  experimental 
data 

H  H 

e.g.,  E^[CH3CH  =  CHCH3J  =  DH  [CHj  CHCH2CH3]  -  DH  [CH3CHCHCH31 

Experimental  values  are  available  for  the  first  bond  energy  but  the  second  bond 
energy  is  calculated  by  group  additivity. 

The  same  principle  is  applied  to  asymmetric  compounds.  The  pi  bond  energy 
for  propene  is  defined  by  the  scheme  : 

CH„  CHCH„  +  H- 


ch3ch2ch3 


where 

E7rtCH3CHCH2]  =  DH.  -  DH2  =  DH,  -  DH2 

In  Table  2-3  are  the  pi  bond  energies  along  with  the  bond  environment  and 
model  compound  used  in  the  calculations. 
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Table  2-3 

PI  BOND  STRENGTHS 


Bond  Type 

Bond  Environment 

Model  Compound 

Bond  Strength 
kcal/ mole' 

C  =  O 

c_c/ 

/° 

CH,C 

82.  6 

C  =  O 

\> 

y° 

H - C 

\* 

CH,0 

75.  0 

o 

il 

u 

x° 

c  —  C 

2 

/° 

CH-C 

74.  5 

c  =  c 

X 

CH2  =  CH  -  O 

3  \ 

H 

CH2CHOH 

52.  8 

c  =  c 

CH2  =  CH2 

ch2ch2 

59.  5 

o 

ii 

o 

C-CH  =  CH-C  CH3CHCHCHCH3 

57. 7 

c  ■•=  c 

C-CH  =  CH2 

CH3CHCH2 

58.  4 

2.  2  UNIMOLECULAR  REACTIONS 

The  majority  of  work  on  individual  reaction  classes  has  centered  on  unimolec- 
ular  reactions.  The  emphasis  on  unimolecular  reactions  was  for  two  reasons. 
First,  these  reactions  constitute  the  simplest  reaction  category,  and  if  a 
trial  method  to  calculate  activation  energies  fails  here,  it  is  not  likely  to  work 
for  more  complex  classes  of  reactions.  Second,  many  bimolecular  reactions 
can  be  treated  as  the  reverse  of  unimolecular  reactions.  Kinetic  information 
calculated  for  the  latter  case  can  then  be  used,  by  means  of  detailed  balance, 
to  compute  comparable  data  for  the  more  complex  bimolecular  reactions. 
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2.  2.  1  Simple  Fission  Reactions 

Simple  fission  of  a  molecule  into  an  atom  and  molecule  or  two  molecules  is 
the  simplest  case  of  unimolecular  reactions.  The  procedure  for  calculating 
the  A  factor  for  these  types  of  reactions  is  well  documented  and  straight  for¬ 
ward.  In  addition,  with  the  assumption  that  the  back  reaction  for  the  associa¬ 
tion  of  two  radicals  or  a  radical  and  atom  has  zero  activation  energy,  the 
activation  energy  for  simple  fission  reactions  is 

E  =  AHg  +  RT 
a  R 

where  AH.^  is  the  endothermicity  of  the  reaction  defined  as 

AH^  =  AH°  [products]  -  AH°  [reactant] 

these  AH^  are  easily  calculated  by  group  additivity  thermodynamic  methods 
using  the  radical  group  values  given  above.  Two  examples  are  : 

1.  C2H6—  2  [.CH  1 

Observed  Ea  =  91.7  kcal/mole  at  T  =  800°K 
Calculated  from  group  additivity: 

AH°  =  2AH°,  800oK  [*CH3i  -  AHf°,  80QoK  [C^l  =  90.  6  kcal/mole 
Ea  =  AH°  +  RT  =  92.  2  kcal/mole 

2.  c2h3ono2-^c2h5o  +  no2 

Observed. Ea  =  41.  2  kcal/mole  at  T  =  420°K  • 

Calculated  from  group  additivity 

AHR  =  AHf°’  298°K  tC2H50.]  +  Hj;  42QoK  [NO^ 

-  AH°  42Qok  [C2H50N021  .  =  41.  6  kcal/mole 
Ea  =  AH°  +  RT  =  42.  4  kcal/mole 

£Y 

The  agreement  in  both  cases  is  good,  and  even  with  this  method  only  applied  to 
a  few  examples,  it  seems  fairly  certain  that  the  method  is  applicable  to  a  wide 
variety  of  simple  fission  cases. 
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Some  possible  problems  can  be  anticipated.  In  many  cases  there  is  a  lack 
of  group  heat  capacity  values  available  to  allow  extrapolation  of  AH°  at  298°K 
to  the  appropriate  reaction  temperature.  In  the  examples  given  above,  this 
extrapolation  has  been  made  for  the  cases  where  heat  capacity  exists,  but  an 
important  deficiency  in  the  group  value  tables  is  this  lack  of  heat  capacity  data. 
This  may  limit  the  effectiveness  of  higher  temperature  calculations. 

One  advantage  of  this  present  method  is  that  the  use  of  the  radical  group 
values  given  above  already  accounts  for  resonance  stabilization  of  the  unshared 
electron  in  the  free  radical  molecule.  This  is  an  improvement  over  earlier 
methods  which  had  to  account  for  resonance  stabilization  by  a  general  corr 
rection  factor  in  a  separate  calculation  step.  But  some  problems  may  be 
encountered  when  conjugated  chains  or  rings  are  available  to  stabilize  the  free 
radical.. 

In  any  case,  since  the  calculation  for  the  A  factors  for  these  types  of  reactions 
is  well  established,  the  inclusion  of  this  method  to  calculate  the  activation 
energy  should  provide  a  successful  method  to  calculate  rate  constants  for 
simple  fission  reactions  without  any  great  difficulty. 

2,  2.  2  Three  Membered  Ring  Reactions 

This  class  of  reactions  represents  a  fairly  involved  reaction  system.  O'Neal 
and  Benson^  have  postulated  a  biradical  mechanism  for  these  reactions,  and 
they  have  developed  a  system  to  calculate  AS^  for  individual  reactions.  In 
addition,  they  have  done  extensive  work  on  a  method  to  estimate  activation 
energies  for  this  class  of  reactions.  This  class  of  reactions  was  studied  in 
this  project  so  as  to  better  understand  the  application  of  group  additivity 
methods  to  the  calculation  of  Arrhenius  parameters  for  more  complex  classes 
of  reaction  and  also  to  see  what  improvements  could  be  offered  on  the  present 
method. 
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The  simplest  reaction  of  this  class  is  the  pyrolysis  of  cyclopropane  to  give 


ch2  =  ch  -  ch3 


Using  a  biradical  mechanism,  the  reaction  can  be  divided  into  several  distinct 
steps. 

1.  Ring  Opening.  For  ring  openings  to  the  corresponding  biradical, 

AH^  can  be  estimated  directly  using  group  additivity  thermodynamic 
methods.  Again,  the  group  values  of  the  radical  groups  are  required 
for  this  step. 

2.  Temperature  Correction.  The  AH°  given  in  the  tables  (lc)  is  for 
298°K.  Heat  capacity  data  has  to  be  used  to  correct  this  to  the 
reaction  temperature.  There  is  an  important  lack  of  heat  capacity 
data  for  radical  groups,  and  this  may  prove  to  limit  the  effectiveness 
of  some  calculations. 

3.  Transition  States.  Modification  o£  the  bi-radical  structure  to  the 

structure  of  the  transition  state  must  be  made  and  the  corresponding 
1* 

AH’  deduced  for  this  step.  This  is  the  important  step  since  the  value 
^  t 

for  AH^  for  this  step  can  only  be  obtained  by  empirical  methods.  In 
addition,  for  a  given  cyclopropane  ring  reaction  several  transition 
states  are  possible,  and  all  of  these  have  to  be  accounted  for  in  this 
step.  As  an  example  the  cyclopropane  decomposition  is  considered 
in  detail. 

Step  1.  Ring  opening  via  a  biradical  mechanism 


CH. 


CH 


/  \ 


CH. 


I 


I 

i 
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Ah\  =  2AH°  be  -  (H2)(C)1  -  2AH°  [C  -  (H2)(C2)]  -  AH°  lAl 

=  54,  0  kcal/mole 

where  AH°  ( A )  is  the  three  member  ring  strain  energy  as  given  by 
Benson(2). 

Step  2.  For  a  reaction  temperature  of  770°K  a  correction  of  0.  8  kcal/ 

t 

mole  must  be  added  to  298°K 

Step  3.  Transition  State,  The  decomposition  of  cyclopropane  yields 
propylene.  The  transition  state  involves  a  hydrogen  atom  migration: 


The  measured  activation  energy  is  65.  6  kcal/mole  at  770°K.  As 
calculated 

Ea  =  AH^  ^  AH?,  ^  ^  =  65.  6  kcal/mole 

f 

All  the  factors  in  this  equation  are  known  except  AH^,  that  is  the 
correction  factor  for  hydrogen  atom  migration  in  the  transition  state. 
Its  value  can  be  deduced  to  be 


=  9.  3  kcal/mole 


The  numerical  values  differ  somewhat  from  those  of  O’Neal  and 
Benson^,  due  mainly  to  different  values  for  the  radical  groups. 


This  deduced  value  for  the  energy  associated  with  hydrogen  atom 
migration  in  a  1,  3  biradical  will  be  used  for  other  3-membered  ring 
calculations.  The  assumption  is  made  that  this  value  can  be  used  at 
all  temperatures. 
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A  second  possibility  is  available  in  step  3  above,  that  is  closure  of 
the  three  membered  ring.  The  reaction  of  the  cis-trans  isomeriza¬ 
tion  of  1,  2-dideutero-cyclopropane  has  been  studied  and  an  experi¬ 
mental  activation  energy  obtained  of  64.  2  kcal/mole  at  720°K. 
Therefore 

AhL  =  8.0  kcal/mole 
U3 

Again,  the  assumption  is  made  to  use  this  value  at  all  temperatures. 

The  biradical  mechanism  is  demonstrated  for  the  case  of  cis-1,  2  - 
dimethylcyclopropane. 

The  reactions  are: 

ch3  ch3 

(1) 

CH -  CH  - 

/ 

CH2 

(2) 

_ ..  CH3CHCHCH2CH3 

(3)  CH3CH2CCH2 
CH3 

(4)  CH  CHCCH 

- -  J  I 


CH. 


CH 


CH 


cis-trans 

isomerization 


\/ 


CH. 


CH. 
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An  important  point  in  three  membered  ring  reactions  is  that  several 
sites  are  available  for  the  biradical,  for  example. 


ch3  c 


CH 


\/ 

CH, 


CH3  .  CH 

CH  -CH 

\/ 

CH, 


CH3 

CH 

CH  — 

-CH 

\/ 

CH2 


CH.  CH- 

i  3  |  3 


•  c: 


H 


c: 


H 


/ 

■  CH2 


The  calculated  and  observed  activation  energies  for  the  above  reactions  are 
shown  in  Table  2-4. 


Table  2-4. 

CALCULATED  AND  OBSERVED  ACTIVATION  ENERGIES 


Product 

Reaction  Temperature 
°K 

E-  (Calculated) 
kcal/mole 

E  (Observed) 
\cal/mole 

1 

689 

61.  4 

59.  4 

2 

689 

62.7 

61.  3* 

3 

689 

64.  1 

61.  9 

4 

689 

64.  1 

62.  3 

^Average 

of  61.  2  and  61.  4  for  cis  and  trans  form  form  of  pent 

-2-ene. 

The  agreement  is  good  for  the  activation  energies.  This  methods  appears  to 
be  generally  applicable  to  these  small  ring  compounds.  O'Neal  and  Benson^ 
have  also  investigated  cyclobutane  reactions  and  these  methods  also  are 
applicable  in  that  case.  Nevertheless,  several  problems  are  evident.  First, 
if  a  multiple  bond  is  available  alpha  to  an  unshared  pair  of  electrons,  delocal¬ 
ization  will  occur.  Presumably,  since  the  group  values  for  radical  groups 
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already  account  for  this  delocalization,  this  can  be  handled  in  a  straightforward 
manner.  The  problem  of  conjugated  chains  or  rings  being  available  for  further 
stabilization  of  the  free  radical  has  not  been  studied.  More  important  though, 
double  bonds  can  allow  for  migration  of  the  unshared  pair.  This  can  lead  to 
new  mechanisms  with  different  transition  states.  For  example,  the  major 
product  of  vinyl  cyclopropane  decomposition  is  cyclopentene. 


CH, 


iCH 


CH. 


CH 


CH  - 


\  / 

CH. 


■CH. 


•  CH 


Vh/ 


■CH. 


CH  CH, 

.  |  2 

CH 

_ _ _ .  *1 

CH, 

•  La 

•CH  >CH2 

.chn 

'CH  ^ 

CH2 

CH  - 

II 

- CH 

1 

CH 

CH 

\  / 
CH2 


This  five  membered  ring  closure  mechanism  can  also  be  treated  in  a  manner 
similar  to  the  three  membered  ring  closure.  This  results  in 

AH*t,  =  8.  3  kcal/mole 

S 

Applying  this  value  of  AH^,  to  the  reaction 


CH,=  CH  CH  _ 

1  1 

CH3 

1 

CH 

- ►  /  \ 

1  1 

CH - CH 

/ 

CH, 

\ 

CH, 

\  / 

1  2 

1  2 

ch2 

CH  = 

=  CH 
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yields  an  activation  energy  of  47,  8  kcal/mole.  The  observed  Ea  is  48.  6  kcal/ 
mole.  Again  the  agreement  is  good,  but  this  simple  example  does  give  some 
indication  of  the  problems  associated  with  complex  side  chains.  They  greatly 
increase  the  total  number  of  possible  products  and  provide  new  transition  states. 
Extensive  experimental  data  will  be  required  to  provide  model  systems  for 
many  possible  transition  states. 

2.  2.  3  Six  Centered  Elimination  Reactions 

A  large  fraction  of  unimolecular  gas  phase  reactions  appear  to  take  place  by 
way  of  four  centered  and  six  centered  cyclic  activated  complexes.  The  vast 
majority  of  the  four  centered  reactions  are  the  elimination  of  hydrogen  halides 
from  the  alkyl  halides  to  produce  olefins.  Losses  of  water,  ammonia  and  H^S 
from  tertiary  alcohols,  amines,  and  mercaptans  have  been  observed  in  only  a 
few  cases.  Since  this  program  is  not  concerned  with  halogen  containing 
compounds  we  will  not  investigate  four  centered  reactions  at  this  time.  Instead 
we  will  investigate  six  centered  reactions  which  involve  esters,  alcohols  and 
ketones.  These  reactions  are  more  applicable  to  the  interests  of  this  program, 
and  it  is  believed  that  the  techniques  developed  for  six  centered  reactions  will 
be  equally  applicable  to  four  centered  reactions. 

Methods  have  been  developed  for  calculating  AS^  for  different  classes  of  six 

centered  elimination  reactions  but  little  work  has  been  done  on  estimating  AH^ 

(3) 

for  these  types  of  reactions 

The  information  and  experience  gained  in  the  study  of  the  two  previous  cate¬ 
gories  of  unimolecular  reactions  were  applied  to  several  six-centered  elimina¬ 
tion  reactions. 

Below  are  given  model  reactions  for  the  various  categories  of  six  centered 
reactions  which  were  studied.  In  all  cases,  the  transition  state  involves  a 
six  membered  ring  containing  bonds  of  fractional  order.  This  results  in  an 
overall  polarization  of  the  molecule  between  the  two  reacting  ends. 
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1.  Esters 


ch3c. 


ch3c 


nOCH2CH3 
Ethyl  acetate 


6"  fit 

/O  .  H 

X 

v  CH, 

\  '/ 

No  CH£ 


C2+H4 

CH3COOH 


2.  a  -  Unsaturated  Alcohols 


^CH2 


,CH2  •  H 


C  H, 
3  6 


ch2-ch2 

But-3-ene-l-ol 


'CH,  .  ch: 

C  i 


H2co 


Allylic  Rearrangement 


^CD2  CD2< 


CH 

\ 


ch2-ch2 


'CD2  *  CDZ, 


‘CH2  .  CH/ 


.CD-,  —CD, 

X  2\ 

CH  CH 

^*CH2  CH2^ 


6,  6  -  Tetradeutero  -1,  5  -  hexadii 


4.  Allylic  Ethers 

XCHz 

CH 

—  ch: 


6"  6+ 
•  CH2  •  H 


O  •  CH  - 


C2H4 


ch3cho 


Ethyl  vinyl  ether 
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5.  Diesters 


/ 

X - CH, 


Q  CH_ 

V  ^  3 


o  •  'c'6+ 


o 

I 

CH. 


"O  —  c 


O  CH2  q 


\ 


CH. 


/ 

O-^-C 

\ 


CH, 


Methylene,  ddacetate 


6.  H  -  atom  transfer 


1,  3  -  hexadiene  (cis) 


.CH-/ 


HC 

Jc 


‘CH 


h6+ 

CH 

/  \ 


HC 

II 

HC 


CH. 


CH 

'CH  '  CH . 


\  . /  \ 


The  assumption  was  made  that  all  bonds  formed  or  broken  in  the  transition  state 
of  the  model  compound  are  of  order  one-half.  The  energy  involved  is  equated  to 
one-half  the  appropriate  bond  strength.  The  various  bond  strengths  required 
were  calculated  by  group  additivity  thermodynamic  methods  as  outlined  above. 

The  six-membered  ring  of  partial  bonds  produces  a  polarization  of  the  reacting 
ends  of  the  six  membered  chain.  The  polarization  energy  is  the  difference 
involved  in  the  energy  required  to  produce  the  partial  bonds  and  the  observed 
activation  energy. 


The  various  model  compounds  listed  above  were  chosen  as  the  simplest 

examples  for  which  experimental  data  were  available.  The  experimentally 

observed  activation  energies  are  listed  in  Table  2-5  along  with  the  correction 

for  end  interaction  AHt  deduced  by  this  method. 

P 
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Table  2-5 

ACTIVATION  ENERGIES  AND  END  INTERACTION  CORRECTION 


Reaction  Category 

Model  Compound 

E 

a 

kcal/mole 

T°K 

AHt 

P 

kcal/mole 

Esters 

Ethyl  Acetate 

48.  0 

725-876 

13.  8 

a  -  unsaturated 
alcohols 

But-3-ene-  l-ol 

41.  0 

643-685 

13.  2 

Allylic 

Rearrangement 

1 ,  1 , 6, 6  -  tetra  - 
deutero  -1,5- 
hexadiene 

35.  5 

530 

4.  4 

Allylic  Ethers 

ethyl  vinyl 
ethe  rs 

44.  4 

640-721 

3.  5 

Diesters 

methylene 

diacetate 

3  6.  4 

493-578 

17.  2 

H-atom  transfers 

1,3-  hexadiene 
(cis) 

32.  5 

474-518 

14.  1 

E  was  calculated  from  the  summation  of  one-half  the  bond  energies  of  those 

cl 

bonds  being  broken  in  the  transition  state,  and  the  subtraction  of  (1)  one-half 
the  bond,  energies  of  those  bonds  being  formed,,  and  (2)  AH^,  to  account  for 
the  stabilization  of  the  transition  state  due  to  the  end  polarization. 

One  important  consideration  for  these  reactions  is  the  effect  of  substituent 
groups  attached  to  the  six  membered  ring.  These  groups  do  not  enter  directly 
into  the  transition  state,  but  have  the  effect  of  adding  or  withdrawing  electron 
density  from  the  ring,  thereby  increasing  or  decreasing  the  magnitude  of  the 
end  interaction.  The  class  of  ester  reactions  was  chosen  to  investigate  sub¬ 
stituent  effects  since  they  encompass  a  large  amount  of  experimental  data. 

The  model  compound  for  ester  reactions  is  the  decomposition  of  ethyl  acetate 
into  ecetic  acid  and  ethylene.  In  studying  compounds  with  hydrocarbon  sub¬ 
stitution  at  the  a  or  P  position  of  the  transition  state  ring,  account  must  be 
taken  of  the  changing  bond  environment  as  well  as  the  change  in  electron 
density  of  the  ring. 
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The  two  cases  of  methyl  substitution  at  the  (3  position  of  the  ester  along  with 
the  nonsub stituted  model  compound,  ethyl  acetate  are  illustrated  below. 


(I)  CH3C 


‘OCH2CH3 


(II)  CH,C 


'OCH(CH3)2 


o 


(III)  CH.C 

•3  N 


'o-c(ch3)3 


CH3COOH 

+ 

CH2CH2 

CH3COCH 

+ 

CH2CHCH3 

CH3COOH 

+ 

CH2C(CH3>Z 


The  calculated  and  observed  activation  energies  are  given  in  Table  2-6  along 
with  the  correction  per  methyl  group. 


Table  2-6 

ACTIVATION  ENERGIES  AND  CORRECTIONS 
FOR  p  METHYL  SUBSTITUENTS 


Reaction 

Observed 

Calculated 

Number 

Correction  per 

Number 

AE-kcal 

AE-kcal 

CH^  groups 

CH^  group-kcal 

i 

48.  0 

48.  0 

0 

1 

ii 

45.  0 

48.  8 

1 

-3.  8 

. 

hi 

40.  0 

49.  4 

2 

-4.  7 

Therefore  an  average  correction  of  -4.  3  kcal  per  CH3  group  will  be  applied. 

In  a  similar  manner,  the  following  are  examples  of  ethyl  substituents  at  the 
P  position;  and  the  corresponding  calculated  corrections  per  ethyl  group.  The 
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activation  energies  and  the  calculated  corrections  per  ethyl  group  are 
presented  in  Table  2-7. 


(IV)  ch3c 


O 


'OCHCH. 


CH2CH3 


CH3COOH 

+ 

CH2CHCH2CH3 


(V) 


o 


CH  C 

~>  V 


'O-CCH. 


(CH2CH3)2 


CH  COOH 
+ 

ch2c(ch2ch3)2 


Table  2-7 

ACTIVATION  ENERGIES  AND  CORRECTIONS 
FOR  (3  ETHYL  SUBSTITUENTS 


Reaction 

Number 

Observed 

AE-kcal 

Calculated 

AE-kcal 

Number  of 
Ethyl  Groups 

Correction  per 
Ethyl  group-kcal 

i 

48.  0 

48.  0 

0 

IV 

45.  2 

48.  8  . 

1 

-3.  6 

V 

38.  0 

49.  4 

2 

-5.  7 

The  average  correction  is  -4.  7  kcal  per  ethyl  group. 


Methyl  substituents  at  the  a  position  are  shown  in  (VI)  and  (VII): 


(VI)  CH,C 

J  N 


sOCH2CH2CH3 


CH3COOH 

+ 

CH2CHCH3 
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*o  ch3cooh 

(VII)  ch3c  - -  + 

X>.<:H2CH(CH3)a  CH2C(CH3)2 


The  data  for  these  substituents  are  shown  in  Table  2-8. 


Table  2-8 

ACTIVATION  ENERGIES  AND  CORRECTIONS 


FOR 

a  METHYL  SUBSTITUENTS 

Reaction 

Number 

Observed 

AE-kcal 

1 

Calculated 

AE-kcal 

Number  of 
Methyl  Groups 

Correction  per 
Methyl  group-kcal 

I 

48.  0 

48.  0 

0 

VI 

47.  7 

46.  8 

1 

+0.  9 

VII 

47,  3 

46.  1 

2 

+0 .  6 

The  average 

correction  is 

therefore  +0. 

7  kcal  per  CH.  group. 

Examples  of  ethyl  substitution  at  the  a  position  are: 


(VHI)  CH3C 

ch3cooh 

r  1 

no(ch2)3ch3 

CH2CHCH2CH3 

ch3cooh 

(IX)  CH3C 

— -  + 

OCH2CH(CH2CH3)2 

CH2C(CH2CH3) 

Data  for  the  ethyl  substituents  at  the  a  position  are  shown  in  Table  2-9. 
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Table  2-9 

ACTIVATION  ENERGIES  AND  CORRECTIONS 
FOR  a  ETHYL  SUBSTITUENTS 


Reaction 

Number 

Observed 

AE-kcal 

Calculated 

AE-kcal 

Number  of 
Ethyl  groups 

Correction  per 
Ethyl  group-kcal 

I 

48.  0 

48.  0 

0 

— 

VIII 

46.  0 

46.  8 

1 

-0.  8 

IX 

45.  8 

46.  1 

2 

-0.  2 

This  results  in  an  average  correction  of  -0.  5  kcal  per  ethyl  group  at  the 
a  position. 

This  method  for  hydrocarbon  substituents  can  be  applied  to  other  classifica¬ 
tions  of  groups  as  well.  The  limiting  factor  is  the  lack  of  sufficient  experi¬ 
mental  data  for  model  compounds. 

Table  2-10  lists  a  variety  of  ester  reactions  to  which  our  six-centered 
reaction  method  was  applied  to  calculate  the  activation  energy.  The  observed 
values  of  AE  are  from  a  tabulation  of  experimental  values  by  Bensonv  The 
values  in  parentheses  are  activation  energies  calculated  from  observed  rate 
constants  and.  Benson's  tabulation  of  A  factors.  In  most  cases  where  a  wide 
discrepancy  exists  between  experimental  values  and  Benson's  "calculated." 
activation  energies,  Benson's  values  are  probably  more  accurate.  The 
values  of  AE  calculated  by  the  present  method  are  in  good  agreement  with 
the  observed  values. 


Table  2-10 

ACTIVATION  ENERGY  FOR  SIX  CENTERED  REACTIONS 


AE-kcal 

Observed 

AE-kcal 

Calculated 

i.  c2h5cooc2  — 

48.  5 

48.  0 

C,H.  +  C0H  COOH 

2  4  2  5 

2.  Isopropyl  acetate  — - 

45.  0 

44.  5 

C0H,  +  CH  COOH 

JO  3 

46.  3 
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Table  2-10 


ACTIVATION  ENERGY  FOR  SIX  CENTERED  REACTIONS  (Continued) 


AE-kcal 

Observed 

AE-kcal 

Calculated 

3. 

t  -  Butyl  acetate  — ► 

(ch3)2c  =  ch2  +  ch3cooh 

40.  5 

40.  0 

40.  8 

4. 

sec  -  Butyl  acetate  — *■ 
cis-trans-2butene  +  CH,COOH 

1-butene  +  CH3COOH 

46.  6 
(45.  2) 
(45.2) 

43.  8 

44.  1 

5. 

t-Butyl  propionate  — 
isobutene  +  ^H^COOH 

39.  16 

40.  8 

6. 

N-Propyl  acetate  — ► 

C.H.  +CH.COOH 

3  D  3 

47.  7 

48.  3 

47.  5 

7. 

N-Butyl  acetate  — *■ 

1-butene  +  CHjCOOH 

46.  0 

46.  3 

8. 

N-Pentyl  acetate  — ► 

1-pentene  +  CH^COOH 

46.  4 

46.  3 

9. 

3-Methylbutyl  acetate  — - 

ch3cooh  +  ch2  =  chch(ch3)2 

47.  9 

46.  3 

10. 

2-Methylpropyl  acetate  — ^ 
isobutene  +  CH3COOH 

48.  8 

47.  5 

11. 

2-Methylbutyl  acetate  — ► 

ch3cooh  +  ch2  +  c(ch3)c2h5 

45.  7 

46.  3 

12. 

2-Ethylbutyl  acetate  — ■ ► 

ch3cooh  +  ch2  =  c(c2h5)2 

49.  7 
(45.  8) 

45.  1 

13. 

2-Pentyl  acetate  — ► 

CH^COOH  4-  2-pentenes  (cis,trans) 

43.  7 
(43.  8) 

42.  6 

CH^COOH  +  1-pentene 

(43.  8) 

44.  1 

14. 

3-Pentyl  acetate  — ► 

CH^GOOH  +  2-pentenes  (cis,trans) 

44.  7 

43.  4 

15. 

2-Heptyl  acetate  — ► 

CH^COOH  +  2 -heptenes  ( cis, trans) 

45.  3 
(43.  6) 

42.  6 

CH3COOH  4-  1-heptenes 

(43.  6) 

44.  1 
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Table  2-10 

ACTIVATION  ENERGY  FOR  SIX  CENTERED  REACTIONS  (Continued) 


AE-kcal 

Observed 

AE-kcal 

C  alculated 

16. 

3-Heptyl  acetate  — ► 

CH^COOH  +  2-heptenes  (cis, trans) 

46.  9 
(42. 8) 

43,  4 

CH^COOH  +  3-heptenes  (cis,  trans) 

(42.  8) 

42.  2 

17. 

4-Heptyl  acetate  — *- 

CH^COOH  +  3-heptenes  (cis,  trans) 

44.  9 

42.  7 

42.  2 

18. 

3-Methyl  -  2-pentyl  acetate  — *- 
3-Methylpent  -  1-ene  +  CH^COOH 

43.  2 
(43.  5) 

44.  1 

3-Methylpent  -  2-enes  (cis, trans) 

+  CH3COOH 

(43.  5) 

43.  3 

19. 

2,  4  -  Dimethyl  -  3-pentyl  acetate  — ► 

(CH3)2  C  =  CH(i  =  C3H7) 

44.  7 
(43.  0) 

44.  1 

20. 

2,3-  Dimethyl  -  2 -butyl  acetate  — 

CH^COOH  +  tetra  methylethylene 

41.  3 
(37.9) 

41.  7 

ch3cooh  +  ch2  =  c(ch3)ch(ch3)2 

(37.  9) 

40.  4 

21. 

3 -Methyl  -  3 -pentyl  acetate  — ► 

CH^COOH  +  3-methylpent  -  2-ene  (cis,  trans) 

42.  4 
(38.  0) 

39.  7 

CH^COOH  +  2-ethylbut  -  1-ene 

(38.  0) 

40.  0 

22. 

2-Methyl  -  2-pentyl  acetate  — ► 

CH^COOH  +  2-methylpent  -  1-ene 

40.  6 

(38.  3) 

40.  4 

CH3COOH  +  1-methylpent  -  2-ene 

(38.3) 

38.  9 

The  above  results  demonstrate  the  systematic  method  as  applied  to  substit¬ 
uent  effects  for  ester  reactions.  A  similar  method  can  be  used  to  gain 
correction  factors  for  oxygen,  nitrogen,  and  halogen  containing  substituents. 
Similarly,  it  is  expected  that  this  method  can  be  applied  to  the  other  classes 
of  six  centered  reactions. 

2.  3  BIOMOLECULAR  REACTIONS 

In  general,  biomolecular  reactions  can  be  separated  into  three  classes  — 
metathesis  or  atom  transfer,  association  reactions,  and  displacement 
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reactions.  Some  of  the  inherent  problems  and  tentative  solution  methods  for 
each  of  these  reaction  categories  are  described  below. 


2.  3.  1  Metathesis  Reactions 

An  example  of  hydrogen  atom  transfer  is: 

H  +  C2H6  —  H2  +  C2H5 

The  transition  state  for  this  type  of  reaction  involves  some  structure  such 
as: 

H 

I 

H-H-*-  C - CH, 

A- 


The  method  of  partial  bond  orders  used,  in  the  six  centered  elimination 
reactions  may  be  applied  here  to  calculate  activation  energies.  Problems 
are  sure  to  be  encountered  in  these  atom  transfer  systems  due  to  unequal 
sharing  of  the  atom  in  question  in  the  transition  state.  But  a  large  set  of 
data  exist  for  many  classes  of  metathesis  reactions  so  empirical  correction 
factors  may  be  deduced. 

A  first  step  in  calculating  an  entropy  of  activation  for  atom  transfer  reactions 
may  be  to  assume  a  constant  A  factor  for  any  one  class  of  atom  transfer.  The 
A  factors  appear  to  be  fairly  constant  within  any  one  class. 

Since  most  atom  transfer  reactions  involve  only  two  bonds  being  simultane¬ 
ously  made  or  broken,  this  category  of  reactions  may  prove  to  be  success¬ 
fully  handled  by  group  additivity  methods.  Most  certainly  some  empirical 
correction  factors  will  be  needed  nevertheless. 

2.  3.  2  Association  Reactions 

Such  association  reactions  as  radical- radical  recombination  may  be  best 
handled  by  calculating  the  back  rate  for  unimolecular  decomposition  and 
using  the  principle  of  detailed  balance  to  calculate  the  desired  rate.  But 


27 


more  complicated  proces ses  of  association  will  require  careful  considera¬ 
tion  of  the  vibrationally  excited  intermediate  produced.  For  the  general 
type  of  reaction: 

A  +  B  *  AB*  *  AB 

the  intermediate  AB*  is  usually  vibrationally  excited  and  unless  this  energy 
can  be  transferred  to  other  degrees  of  freedom,  the  intermediate  will 
decompose  back  to  reactants.  The  rate  of  deactivation  of  this  type  of  excited 
intermediate  will  be  an  important  step  in  solving  reaction  rate  calculations 
for  these  bimolecular  reactions.  A  potential  approach  for  estimating  the  rate 
of  decomposition  of  larger  intermediates  is  in  the  Rice-Ramsperger-Kassel 
model  which  treats  the  decomposition  rate  of  an  activated  molecule  in  terms 
of  transition  probability  and  the  density  of  vibrational  quantum  states.  These, 
in  turn,  are  a  function  of  certain  estimated  or  assigned  variables  such  as  the 
critical  energy  or  minimum  quanta  of  vibrational  energy  required  for 
decomposition.  The  low-pressure  limit  for  deactivation  can  also  be 
calculated. 

For  small  (di-  or  triatomic)  intermediates,  it  may  be  possible  and  indeed 
desirable  to  utilize  the  more  accurate  Rice-Ramsperger-Kas sel-Marcus 
model. 

2.  3.  3  Displacement  Reactions 

This  type  of  reaction  appears  frequently  in  the  chemistry  of  unsaturated 
compounds.  For  example: 


It  may  be  possible  to  calculate  activation  energies  for  these  reactions  by 
assuming  a  stable  structure  for  the  transition  state  and  using  group  addi¬ 
tivity  methods  to  calculate  the  AH  of  this  species.  Free  radical  group 
additivity  data,  such  as  that  estimated  during  the  course  of  this  study,  would 
be  utilized  to  compute  the  heat  of  formation.  In  order  to  account  for  differ¬ 
ences  in  bond  strengths,  chemical  environments,  etc.  /  correction  factors 
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will  most  probably  have  to  be  determined  for  various  classes  of  compounds, 
as  in  our  treatment  of  unimolecular  reactions.  These  calculations  may  be 
hindered  by  the  fact  that  little  quantitative  data  seem  to  be  available  for  this 
type  of  reaction. 
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Section  3 


SCHEME  FOR  COMPUTER  CALCULATION  OF  THE 
MOST  PROBABLE  REACTION  RATE 

For  any  elementary  reaction  to  which  there  are  several  possible  reaction 
paths,  and  therefore,  s everal  pos sible  products ,  the  most  straightforward 
way  to  compute  the  most  probable  path  would  be  to  calculate  the  rate  con¬ 
stant  for  all  paths,  and  eliminate  any  paths  which  have  a  rate  of  less  than 
10  percent  of  the  fastest  path.  This  appears  to  be  a  reasonable  criterion  by 
which  to  predict  the  major  products  of  a  reaction.  To  save  computer  cal¬ 
culation  time  a  preliminary  screening  method  may  be  possible.  Since  many 
reactions  of  a  given  molecule  along  similar  paths  should  have  similar  AS^, 
it  may  be  possible  to  calculate  AH^  and  eliminate  all  but  the  most  important 
paths  having  the  smallest  activation  energies.  Then  the  AS^  could  be  cal¬ 
culated  for  only  these  most  important  paths. 
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Section  4 

METHODS  OF  HANDLING  COMPLEX  REACTION  SYSTEMS 

Complex  reaction  systems  involve  multiple,  parallel,  or  successive  elemen¬ 
tary  reaction  steps.  For  this  study,  complex  reaction  systems  are  divided 
into  two  general  categories:  Additive  Reaction  Systems  and  Successive 
Reaction  Systems. 

4.  1  SUCCESSIVE  REACTION  SYSTEMS 

These  systems  are  defined  as  reaction  systems  which  require  one  or  more 
successive  elementary  steps  to  go  from  reactant  to  product.  This  is  repre¬ 
sented.  by: 

A — -A,  — ►  A  — ►  A  — *  A1  (1) 

12  n 

This  reaction  group  can  be  very  complicated  due  to  the  fact  that  at  each 
elementary  step  different  products  may  be  possible.  This  will  produce 
branching  reaction  systems  leading  to  different  products.  Parallel  processes 
are  also  possible.  In  addition,  the  different  products  at  each  individual  step 
may  influence  other  steps.. 

Fortunately  there  are  some  basic  considerations  which  make  it  possible  to 
sort  out  those  pathways  which  are  important.  In  any  chain  reaction  of  the 
sort  described  in  (1)  it  is  necessary  that  the  intermediate  or  chain  propaga¬ 
tion  steps  be  very  fast.  This  is  necessary  if  the  chain  reaction  is  to  compete 
favorably  with  the  concerted  mechanism.  For  example  in  the  decomposition 
of  H^  and  Br^,  the  following  conce rted  m echanism  is  possible: 

H  -  -  -  B  r  H  —  B  r 


i 

i 


H - Br 


(2) 
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The  chain  reaction  would  be 


H2 

+  Br  — ►  HBr 

+  H 

(3) 

H 

+  Br,  — -  HBr 

+  Br 

(4) 

In  order  that  several  successive  steps  can  compete  with  a  concerted  mecha¬ 
nism,  the  chain  propagation  steps  must  have  low  activation  energies.  In 
addition,  the  intermediates  (A,.,  A~,  --  A  ,  in  reaction  (1)  )  must  be  highly 
reactive  intermediates,  i.  e.  ,  atoms  and  free  radicals.  This  now  limits 
somewhat  the  possible  pathways. 

4.  1.  1  Criterion  for  Favorable  Chain  Propagation  Steps 

Following  the  recommendations  of  Benson  we  will  limit  the  possible 

chain  propagation  steps  to  three  types  of  elementary  reactions:  addition, 
atom  abstration,  and  unimolecular  fission.  In  addition,  Benson  has  proposed 
a  criterion  for  the  activation  energies  of  chain  propagation  reactions.  For 
bimolecular  steps  AE/0  <  5.  5  and  unimolecular  steps  AE/9  <14  where  0  = 

2.  303  RT  kcal/mole.  This  is  an  estimation  of  the  activation  energy  required 
so  that  the  chain  propagation  steps  will  be  rapid  enough  to  compete  with  a 
concerted  mechanism. 

4.  1.  2  Method  for  Determining  Chain  Propagators 

In  order  for  a  chain  reaction  to  proceed,  some  catalytic  agent  must  be  pro- 
duced  and  then  used  to  propagate  the  chain,  the  overall  concentration  of  the 
propagator  not  increasing  with  time.  Benson  has  proposed  a  general  method 
to  find  the  chain  propagator.  We  have  tested  this  for  the  decomposition  of 
ethyl  chloride  and  found  that  it  succes sfully  predicts  the  correct  reaction 
mechanism. 

The  procedure  is  the  following:  Some  radical  R  is  chosen.  It  can  only 
abstract  an  atom: 

—  .  CH2CH2C1  +  RH 

R-  +  C-jH-Cl  —  •  CH  'CHC1  +  RH 
2  5  3  -  * 

— -  CH3CH2  +  RC1  (5) 
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There  are  three  possible  propagation  steps.  The  exothermicity  of  the 
reactions  is  given  in  kcal/mole. 


•  CH2CH2C1 

—  CH2CH2 

+  Cl* 

21.  3 

(6) 

CH3CHC1 

— r  CH2CHC1 

+  H* 

43.  1 

(7) 

ch3ch2 

—  CH2CH2 

+  H* 

39.  0 

(8) 

Clearly,  considering  that  activation  energy  E  >  AH,  reactions  (7)  and  (8) 
are  too  slow  to  be  considered,  only  reaction  (6)  will  be  the  propagation  step. 
Reaction  (6)  produces  the  radical  Cl,  so  in  order  for  this  to  be  a  chain 
process,  the  unspecified  chain  carrier  R  in  reaction  (5)  must  be  Cl,  The 
possible  reactions  with  their  AH  are  given  below. 


Cl* 

+ 

C2H5C1  — 

•  CH2CH2C1  +  HC1 

-5.  0 

(9) 

Cl* 

+ 

C2H5C1  — 

CH3CHC1  +  HC1 

-8.  4 

(10) 

Cl* 

+ 

c2h5ci  — 

ch3ch2  +  Cl2 

+21.  3 

(ID 

Again  only  reactions  (9)  and  (10)  can  compete,  and  as  shown  above  only  the 
decomposition  of  •  CH^CH^Cl  can  compete.  The  final  reaction  system  is 
given. 

Cl-  +  CH3CH2C1  —  •  CH2CH2C1  +  HC1  (12) 

•CH2CH2C1  —  CH2CH2  +  Cl*  (13) 

Therefore,  the  only  products  of  CH3CH2C1  decomposition  are  CH2CH2  and 
HC1  ,  and  we  have  deduced  the  important  elementary  steps  needed  to  calcu¬ 
late  the  overall  rate.  Unfortunately  it  has  been  found  difficult  to  apply  this 
method  to  more  complex  systems.  A  fair  amount  of  chemical  knowledge 
seems  necessary  to  make  a  preliminary  choice  of  important  reactions; 
otherwise  the  number  of  possible  reactions  becomes  difficult  to  handle. 

It  is  necessary  to  include  the  chemistry  of  initiation  and  termination 
reactions  to  develop  a  realistic  list  of  possible  reactions. 
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4.  1.  3  Initiation  Reactions 


and  disappear  two  at 


For  homogeneous  systems,  monoradicals  must  appear 
a  time.  It  will  be  necessary  to  test  all  possible  sources  of  radicals  in  a  given 
system,  the  most  rapid  formation  being  the  one  observed.  It  is  possible  to 
develop  an  appropriate  criterion  to  test  some  of  the  possibilities. 


Shown  in  Table  4-1  are  a  factors  and  activation  energies  for  a  variety  of 
unimolecular  reactions  involving  the  fission  of  a  molecule  into  two  radicals. 
While  there  is  considerable  range  for  the  activation  energies,  logA  =  16 
±1  sec'  for  many  systems.  In  contrast,  bimolecular  reactions  of  two 
molecules  to  two  radicals  (the  reverse  of  radical  disproportionation)  have 
logA  =  9  ±  1  liter /mole. sec. 


Table  4-1 

FISSION  OF  MOLECULES  INTO  TWO  RADICALS 


Reaction 

Log  A 

E(kcal/mole) 

T(°K) 

C,H,  — *  2CH. 

Zb  j 

16.  0 

86.  0 

850 

C4H10  *  2C2H5 

17.  4 

82 

400 

CH3OOCH3  - r  2CH30 

15.  4 

36.  1 

410 

C2H5ONO  — -  C2H5°  +  NO 

14.  2 

37.  7 

430 

C,H  ONO,  — ►  C,HcO  +  NO, 

16.  9 

41.  2 

420 

CH3N  =  NCH3  — ►  CH3  +  NNCH3 

17.  2 

55.  5 

600 

C  ,  H  -  CH  C  ,H  — ►.  C,H_  -  CH,  +  C,Hc 
65  225  65  2  25 

14.  9 

68.  6 

950 

c(no2)4  —  c(no2)3  +  no2 

17.  5 

40.  9 

470 

(le) 

Considering  these  values  of  pre-exponential  factors  Benson  has  shown 
that  for  the  rates  of  unimolecular  and  bimolecular  fission  to  be  equal,  E^  - 
E^  =  (9.  4  ±  2)  9,  where  9  =  2.  303  RT  kcal/mole,  and  Ej  is  the  activation 
energy  for  the  unimolecular  step  and  E^  is  the  activation  energy  of  the 
bimolecular  process.  At  2000°K  E^  -  E^  =  88  ±  18  kcal.  This  is  so  high  that 
bimolecular  initiations  will  almost  never  happen.  At  1000°K  E^  -  E^  =  44  ± 

9  kcal.  It  will  be  necessary  to  test  possible  bimolecular  and  unimolecular 
processes  with  this  criterion. 
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In  most  case,  the  bimoiecular  initiation  process  will  never  occur,  except 

( le) 

with  few  exceptions.  One  example  given  by  Benson  is  the  pyrolysis  of 
olefins. 

C3H6  —  C3Hs  +  H  AHj  =  87  kcal/mole 
2C3H6  —  C3Hg  +  C3H?  AH2  =  48  kcal/mole 

In  this  case  near  1000°K  the  bimoiecular  process  must  be  included  in  the 
kinetics. 

4.  1.  4  Termination  Reactions 

The  rates  for  several  representative  radical  recombination  reactions  are 
given  in  Table  4-2. 


Table  4-2 

RADICAL -RADICAL  ASSOCIATION  REACTIONS 


Reaction 

Log  k 

(Liter/mole  sec) 

T(°K) 

2CH3  —  C2H6 

10.  5 

400 

2C2H5  —  C4H10 

10.  4 

400 

CH3  +  C2H5  -  C3H8 

10.  7 

400 

CH3  +  n  -  C3H7  —  n  -  C^ 

10.  7 

400 

2  i  -  C3H7  —  C6H14 

10.  8 

400 

2  t  -  C4H9  —  C8Hlg 

9.  5 

373 

2  t  -  C4H90  —  CgH1902 

00 

00 

420 

CH3  +  NO  — *-  CH3NO 

00 

00 

300 

Most  of  these  reactions  fall  in  the  range  of  Log  k  =  9.  8  ±  0.  5  liters /mole  sec. 
Disproportionation  reactions  for  large  radicals  also  have  the  same  rate  range. 
To  a  first  approximation  all  radical- radical  recombinations  or  disproportion¬ 
ations  will  be  considered  to  have  nearly  the  same  rate.  As  the  radical  chain 
carriers  build  up  in  concentration,  the  termination  reactions  will  be  governed 
by  this  rate. 
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The  only  exception  to  this  case  is  the  recombination  of  two  small  radicals 
(CH^,  NO,  Cl).  In  certain  pressure  regions  these  reactions  will  be  third 
order.  These  molecules  will  have  to  be  treated  as  special  cases.  Possibly 
the  number  of  small  radical  reactions  involved  in  most  reactions  can  be 
minimized,  and  specific  experimental  data  used  for  these  reactions. 

4.  2  ADDITIVE  REACTION  SYSTEMS 

These  systems  are  defined  by  a  molecule  or  molecules  in  which  two  or  more 
distinct  elementary  steps  are  possible  with  the  initial  reactants.  This  can  be 
represented  by  the  following:  If  reactive  center  A  can  undergo  an  elementary 
step  to  produce  A1,  and  similarly  B  to  B  '  ,  then  an  additive  reaction  system 
would  contain  both  A  and  B,  that  is: 

A  —  B  — ►A'  —  B' 

This  represents  the  additive  reaction  class  in  its  simplest  form.  If  both 
reactive  center  A  and  B  are  present  in  a  molecule,  then  the  product  of  the 
reaction  is  just  the  sum  of  the  reactions  of  the  two  reactive  centers.  An 
additional  mechanism  may  be  possible.  We  will  term  this  the  interference 
mechanism.  That  is  to  say  with  both  A  and  B  present  in  the  same  molecule, 
the  possible  elementary  steps  which  lead  to  products  will  be  something 
different  than  A  to  A'  and  B  to  B'.  In  other  words,  the  simultaneous  occur¬ 
rence  of  both  A  and  B  may  interfere  with  the  "usual"  elementary  reactions. 
For  instance  B  to  B'  may  be  blocked  by  the  presence  of  A. 

A  —  B  — -  A1  —  B 

Or  another  possibility  is  that  A  and  B  together  may  produce  some  different 
elementary  step. 

A  —  B  — ►  A'  —  C 

It  is  clear  that  in  such  cases  as  two  carbon-carbon  double  bonds  widely 
separated  in  a  hydrocarbon  chain  so  that  there  is  no  conjugation,  the 
appropriate  chemical  entities  will  undergo  their  usual  reaction,  unhindered 
by  the  presence  of  each  other.  This  would  be  a  simple  additive  case.  As  long 
as  we  are  able  to  calculate  the  individual  elementary  steps,  the  simple  addi¬ 
tive  case  will  just  be  the  sum  of  the  possible  elementary  steps.  But  it  will 
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be  necessary  to  investigate  more  thoroughly  the  possible  "interference" 
mechanism  possible  for  different  combinations  of  A  and  B. 

For  example  the  molecule 


H2C' 


H2C 


CH-CH2-C 


0-CH2CH3 


would  be  expected  to  undergo  reactions  specific  to  both  3-membered  rings 
and  to  esters.  The  presence  of  the  ring  will  not  influence  to  any  great  extent 
the  six  centered  elimination  pathway  of  the  ester;  nor  will  the  ester  group 
interfere  with  the  biradical  mechanism  of  the  small  ring.  The  final  prod¬ 
uct  distribution  depends  only  on  the  relative  rates  of  the  two  independent 
pathways  and  the  reaction  conditions.  This  is  a  simple  example  of  the  addi¬ 
tive  reaction  class. 


A  more  complex  situation  arises  when  "interference"  between  two  distinct 
reaction  centers  can  occur.  An  example  of  this  type  are  the  reactions  of 
bicylopropyl.  The  usual  biradical  mechanism  is  present  for  each  ring 
separately  and  4  mono-  and  9  di-ene  products  are  found.  Some  examples 
are  given  below. 
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All  these  products  arise  due  to  ‘'addition11  of  the  pathways  possible  from  the 
two  distinct  reactive  centers.  But  several  pathways  are  possible  only  if  both 
reactive  centers  are  present.  These  reactive  pathways  are  not  due  to  simple 
’'addition"  of  reactive  pathways.  For  example 


This  pathway  in  fact  produces  the  major  product  for  the  bicyclopropyl  mol¬ 
ecule.  Another  example  of  an  interference  process  is  when  a  double  bond  is 
adjacent  to  a  3-membered  ring.  The  usual  reactive  pathways  of  the  ring  and 
double  bond  are  available,  but  the  following  reaction  produces  one  of  the 
major  products. 


This  five-membered  ring  closure  mechanism  is  not  possible  unless  both  the 
3-membered  ring  and  double  bond  are  present  together. 

We  have  not  investigated  the  possible  interference  mechanisms  for  a  wide 
variety  of  systems.  But  these  few  examples  indicate  that  they  are  of  major 
importance.  Some  complex  molecules  which  contain  more  than  one  reactive 
center  will  be  able  to  be  treated  by  simple  addition  of  our  present  methods. 
But  additional  study  will  be  required,  to  locate  important  classes  of  inter¬ 
ference  mechanisms.  Due  to  all  the  possible  combinations  between  reaction 
categories,  there  may  be  many  special  cases  which  require  special 
treatment. 
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Section  5 


POTENTIAL  PROGRAMMING  METHODS 
FOR  COMPUTER  APPLICATION 

Two  examples  have  been  chosen  to  demonstrate  the  application  of  computer 
methods  to  the  reaction  kinetics  technique.  These  are  simple  fission  reactions 
and  reactions  of  simple  fission  with  resonance  stabilization. 

5.  1  SIMPLE  FISSION  REACTIONS 

The  conditions  for  this  reaction  category  are  that  the  products  of  a  unimo- 
lecular  reaction  are  radicals  and/or  atoms:  and  no  concerted  reaction  path 
with  a  much  lower  activation  energy  is  available.  The  activation  energy  at 
temperature  T,  is  defined  by  E  =  ^  +  RT,  where  ^  is  the  bond  dis¬ 

sociation  energy  of  the  reaction  coordinate.  For  the  reaction 

A-  B  — ►  A*  +  •  B 

g  =  -AH^(A-B)  +  AH^(A-  )  4-  AH^(B*  ).  The  assumption  here  is  that  the  back 
reaction  for  the  recombination  of  two  radicals  has  zero  activation  energy. 

The  calculation  for  the  entropy  of  activation  is  straightforward.  The  various 

contributions  to  the  entropy  of  activation  are  (a)  reaction  path  degeneracy, 

1/3 

R  ln(g)  (b)  external  rotation  which  can  be  shown  to  be  R  ln(6D^  g/RT) 

except  for  the  case  of  H  or  D  atom  cleavage  where  it  is  zero,  (c)  hindered 

internal  rotation  to  free  rotation  in  the  transition  state,  AS^  =  -  S^  (these 

values  are  tabulated  by  Benson  ),  (d)  reduction  in  frequency  of  bending 

modes,  =  5  e.u.  for  the  case  of  atomic  and  free  radical  products,  = 

10  e.u.  when  the  products  are  two  radicals,  (e)  entropy  associated  with 

reaction  coordinate,  AS  =  S  where  S  is  the  entropy  of  stretching  of 

rc  st  st 

a  bond  undergoing  cleavage. 
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The  entropy  of  activation,  therefore  is 

ASt  =  RLng  +  RLn  (6DA_B/RT)1 /3  +  ASrc  +  Sbf  +  ASfr 

and  the  calculated  rate  constant  is 

k  =  (~^“ )  exp  (AS^/R)  exp  (-E/RT) 

This  method  for  the  calculation  of  the  rate  constant  of  a  simple  fission 
reaction  can  be  adapted  to  computer  programming.  This  is  represented 
below  in  Figure  5-1. 

5.  2  REACTIONS  OF  SIMPLE  FISSION  WITH  RESONANCE  STABILIZATION 
This  case  is  very  similar  to  the  simple  fission  case  with  some  exceptions. 
The  conditions  are  the  same  with  the  addition  of  the  fact  that  a  C-C  pi  bond, 
is  one  C-C  single  bond  removed  from  the  reaction  site.  This  means  that 
cleavage  results  in  an  atom  and  olefin  or  radical  and  olefin.  The  carbon- 
carbon  double  bond  stabilizes  the  radical  product  by  resonance  delocaliza¬ 
tion.  This  mechanism  produces  changes  in  both  activation  energy  and  entropy 
from  that  in  the  simple  fission  case. 

The  activation  energy  is  calculated  in  the  same  way  as  in  the  simple  fission 
case,  but  it  is  necessary  to  use  the  group  values  for  the  radical  groups 
developed  earlier  in  this  program.  These  group  values  already  account  for 
resonance  stabilization  of  the  radical. 

The  entropy  calculation  is  similar,  but  the  participation  of  the  ir  bond  in  the 
transition  state  must  be  accounted  for.  This  bond  contributes  mainly  to  the 
internal  rotation  where  a  1/2  pi  bond  torsion  is  assigned  -3.  5  e.u.  Also  the 
entropy  decrease  for  the  tightening  of  the  C-C  single  bond  due  to  resonance 
is  obtained  by  increasing  its  normal  vibrational  frequency  by  about  1.  3.  The 
resulting  computer  flow  diagram  is  given  below  in  Figure  5-2. 

It  would  be  a  simple  matter  to  integrate  these  two  schemes  together  to  treat 
simple  fission  cases  with  and  without  resonance  stabilization. 
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GO  TO  SIMPLE 
FISSION  CASE 
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